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Wavelengths and estimated intensities of approximately 2,000 spectral lines emitted by 


singly ionized rhenium atoms were available for this analysis. 
for 220 Re 1 lines facilitated finding and interpreting the first energy levels. 
19 low-even, and 85 high-odd levels were found. 


continued until 


Observed Zeeman patterns 
The analysis 
Transitions between these 


groups of levels account for 1,014 Re 1 lines including about 70 percent of the total observed 


intensity. 
1. Introduction 


The term analysis of the first spectrum of rhenium 
(Re 1) progressed intermittently since it began In 
1931; it culminated in a recent report [1]* on the 
resent state of its analysis and interpretation. 
F wo decades elapsed before the term analysis of the 
second spectrum of rhenium (Re u) could begin. 
The delay is explained by the lack of basic observa- 
tions on this spectrum. No useful description of 
the Re u spectrum was available until 1952 when 
Meggers [2] published mM description of the are and 
spark spectra of rhenium including wavelength and 
intensity data for approximately 6,000 lines, about 
4200 of which were assigned to Ret, and nearly 
1.800 to Rew (and possibly Rew). The wave- 
lengths of spark-enhanced rhenium lines reported 
in that paper ranged from 2000.47 A to 6026.60 A. 
The spark spectrum of rhenium, like that of any 
heavy metal, is rich in the ultraviolet but gradually 
fades in the visible region, and vanishes in the 
orange. Because this spectrum is strong in the 
ultraviolet it would be expected to extend far 
beyond the violet limit of observation with large 
spectrographs at atmospheric pressure. In order 
to record the spectrum in the shorter-wave range 
Catalan made are and spark spectrograms with a 
vacuum spectrograph at Princeton University; 
these spectrograms were measured in Madrid by 
Sales who thus extended the wavelength and in- 
tensity data of rhenium spectra to about 1500 A. 

The splitting of spectral lines by magnetic fields 
(Zeeman effect) is extremely useful in analyses of 
jcomplex spectra because it reveals (by the number 
lof components) the inner-quantum numbers (J) of 
both levels, and usually identifies (by the splitting 
factors, g) the term types (/) and the multiplicity (7). 
In 1940, Zeeman spectrograms of rhenium were 
made [3] with the Bitter magnet and grating spec- 
trographs in the Spectroscopy Laboratory of the 
Massachusetts Institute of Technology, but these 
spectrograms never yielded any information about 
Reu. By courtesy of G. R. Harrison in 1949 
Meggers made new Zeeman spectrograms with the 
M.I. T. equipment and measured them at the 
National Bureau of Standards. Wavelengths of 
' Deceased. 
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The ground state of Re* is described by the spectral term 5d5(a®S)6s a78y. 


the Zeeman patterns were computed, and from these 
the inner-quantum numbers and the magnetic 
splitting factors of the first Re u energy levels were 
derived in 1950. <A preliminary announcement of 
the Zeeman effect and first spectral terms of Re 
was made [4]in 1952. When 22 low-even and 20 high- 
odd levels were established, the entire mass of 
Zeeman data was transferred to Madrid where 
Catalin and Sales continued the analysis until 1954 
when it included about 830 classified lines resulting 
from transitions between 47 low-even and 69 high- 
odd levels. Most of the Zeeman data were exploited 
in this analysis. A final search among wave- 
numbers, at the Bureau disclosed several additional 
(provisional) levels that combined 7 or more times 
with other levels, so that a total of about 1,000 Re u 
lines could be classified. Details concerning the 
classified lines of Re 11, the energy levels, the Zeeman 
patterns and magnetic splitting factors, are given 
below, followed by a brief discussion of spectral 
terms, and electron configurations. 


2. Basic Data of Re II 


The basic data of Re 1 consist of measured wave- 
lengths and estimated relative intensities of more 
than 2,000 lines emitted by electric ares and sparks 
between metallic electrodes. Nearly 1.800 of these 
have wavelengths between 2000 A and 6000 A which 
have been published by Meggers [2]; more than 1,000 
lines of shorter wavetength (to 1500 A) were observed 
with a vacuum spectrograph by Catalan, but only 
the classified lines will be reported in this paper. 
The wavelength measurements were complicated by 
the hyperfine structure (hfs) exhibited by many 
rhenium lines. This has been discussed in the paper 
Term Analysis of the First Spectrum of Rhenium 
[1] and need not be repeated here. It suffices to 
say that in the cases of resolved hfs we have quoted 
the average wavelength published by Meggers [2]. 

By comparing the intensities of lines emitted by 
ares and sparks, it was usually possible to identify 
the lines emitted by neutral atoms and by ions. 
In general, the intensities of Re 1 lines are greater in 
are than in spark spectra, whereas Re 1 (and pos- 
sibly Re 111) lines appear stronger (or only) in spark 
spectrograms. Unfortunately there are no reliable 
criteria for differentiating Re 1 and Re 11, so all the 
lines enhanced in spark spectra above 2000 A were 
included in the Re 1 list although a few may belong 








to Reus. In the vacuum ultraviolet the are spec- 
trogram showed many more lines than the spark and 
there was insufficient information to assign the lines 
to successive spectra; for this reason only those lines 
that appear to satisfy combinations of Re 1 levels 
are reported. 

The Zeeman-effect spectrograms at 81,000 oersteds, 
previously described [1, 4] yielded information con- 
cerning some 220 Re 11 lines. The Zeeman patterns 
of 85 of these were adequately recorded and resolved 
to yield J- and g-values for both atomic energy 
levels responsible for producing a spectral line. 
These initiated this analysis and interpretation of the 
Re spectrum. The Zeeman patterns are of types 
1, 2, 3, and 7 according to the description given by 
Back and Landé [5]. Type 1 occurs when the J- 
values are unequal and the level with the larger J 
has the smaller g; the strongest p-component is not 
displaced, and the strongest n-components have the 
smallest displacement. Type 2 occurs when J and 
q change in the same direction, in which case the 
strongest p-component not displaced, but the 
strongest n-components have the largest displace- 
ment. Type 3 results from equal J-values but un- 
equal g-values which produce a maximum displace- 
ment of the strongest p-components, and the dis- 
placements of two pairs of strong n-components are 
exactly equal to the g-values of the combining levels. 
Type 7 has only one (undisplaced) p-component, 
and two displaced n-components, either because 
one of the J-values is zero, or the two g-values are 
equal. , 

When Zeeman patterns of types 1, 2, and 3 were 
not resolved, only the wavelengths corresponding 
to the maximum intensity in both polarizations 
were measured. Then if the line is classified and the 
q-value of one level is known, the other may be 
calculated. 


3. Classified Lines of Re II 


The classified lines of Rew are listed in table 1 
where information is given as to estimated relative 
intensity [2], wavelength [2], vacuum wavenumber, 
spectral-term combination, and Zeeman effect. In 
column 1, the intensity numbers of many lines shorter 
than 2000 A are followed by the letter @ which indi- 
cates that the line was observed only in the are 
spectrum. Because it is impossible to assign term 
symbols to many Re 11 energy levels these are repre- 
sented in column 4 by the first 4 digits of their 
numerical value referred to the ground state as 
zero. The daggers attached to 7 wavelengths in 
column 2 mean that the line was also classified as 
Re [1]. An asterisk attached to the Zeeman type 
in column 5 means that the Zeeman effect disagrees 
with the term combination of column 4. Numbers 
in column 6 express the separation of the resolved 
Zeeman components, or Ag, in Lorentz units. In 
columns 7 and 8, instead of giving all the observed 
Zeeman components, only the strongest parallel (p) 
and normal (n) components are shown; this infor- 
mation, combined with the g-values of columns 9 
and 10, and the corresponding J-values in column 4, 
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permit the calculation of the complete 


pattern, h 
columns 9 and 10, the g-values in parentheses » 


mean values from resolved patterns adopted t 
compute the other g-values from unresolved pat 
when the types are indicated in column 5 and 


J-values are known from column 4. After th, 
Zeeman data of Re 1 were fully exploited in est, 
lishing energy levels some further progress was 

in searching for additional levels by means of gop 
cidences between wavenumber differences and know 
level intervals. About 830 Re lines were thy 
explained by the end of 1954, but 900 lines wig 
wavelengths greater than 2000 A still remaing 
unclassified. In 1955 this search was continued » 
the Bureau by Jack Tech who found many prov. 
sional odd levels, nine of which have been Tetaingi 
in this analysis. In the meantime the possibjliy 
of applying electronic digital computers to sysi 
matic searches for atomic energy levels in comple 
atomic spectra had been tested with Tc 1 by Bozmy 
[6] who used the IBM-604, and with Rut ly 
Kessler, Prusch, and Stegun [7] who employed th 
SEAC. Since more than half of the known Rey 
lines (mostly of low intensity) were still unclassifig 
it appeared desirable to enlist the aid of an electron; 
computer in a final search for energy levels. Thi 
was done with an IBM-—704 computer which ds 
placed the machines previously used at the Bureay 
Coding this problem for IBM-—704 was successfull 
carried out by W. R. Bozman and C. Coleman, an 
searches for new levels of Re 1 were made by them 
The first computer search was for possibly undis 
covered high-odd levels between 40000 and 8000 
K (=em~') that could combine with the 47 low. 
even levels (0 to 37000 K.), according to the selection 
rule AJ=0, +1. The second search started with 
intervals between the known odd levels (44000 t 
77000 K) and sought new even levels, both lower 
and higher, that might combine with those @ 
opposite parity and thus account for observed Re! 
lines. These systematic computer searches wen 
both disappointing and illuminating. Several hu- 
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dred hypothetical levels were found, but when ther 
combinations were displayed in a square array, they 
usually appeared to have a random distribution, 


they generally involved only 6 or 7 lines, and the) 


often used the same line in two or more places. In) 


other words, the great majority of predicted levels 


appeared to be based on accidental coincidences 
This conclusion was confirmed by calculation of the 
number of probable coincidences to be expected 
among 900 unclassified Re mu lines. 

S. A. Goudsmit is credited [8] with having devised 
the following approximate formula for the most 
probable number of accidental coincidences, ( 
caused by a random distribution of L spectrum lines 
extending over a range of J wavenumber units 
C=2¢(L?/1) x (I—d)/I if all pairs of lines having 4 
wavenumber difference between (d—e) and (d+e) are 
regarded as coincident. If d is small compared with 
I, the second term may be neglected. Given 900 
Re u lines with wavenumbers between 20000 ané 
50000. K, a tolerance of +0.2 K for coincident 
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yields 6 or more chance coincidences for each value 
of d, depending on the size of d. Hundreds of such 
coincidences were found among the unclassified 
Re m lines but the energy levels calculated from them 
were discarded unless they accounted for lines in a 
nonrandom fashion, or combined with at least 8 
other levels. ‘This resulted in the discovery of only 
9 additional low-even levels, and 3 odd levels. 
The search for high-even levels revealed nothing 
that could not be explained by chance coincidences. 
For example, the fact that such tentative high-even 
levels usually combined only with high-odd levels 
indicated that they were fortuitous because the 
density of unclassified Re 1 lines happens to be 
oreatest below 2500 A. This was disappointing 
hecause it destroyed any hope of finding one or more 
spectral series that might be extrapolated to its 
limit and give an approximate value of the ionization 


TaBLE 1. Classified lines of the second spectrum of rhenium (Re 11) 





potential of Re*. This experience is a warning 
that many false levels may be found in complex 
spectra by searching with wavenumber intervals, 
unless the tolerance on coincidences is reduced 
sufficiently to exclude most accidentals, or unless 
the levels can be verified by Zeeman data. 

When frustration ended the search for energy 
levels it was seen that the accepted levels account 
for 1,014 Re 1 lines listed in table 1. Only 33 of 
these lines are doubly classified. Although only 
half of the rhenium spark lines with wavelengths 
greater than 2000 A are classified in this analysis, 
these account for about 70 percent of the total 
observed intensity, and include almost all lines of 
intensity greater than 30. The combinations 
a’S,—z'P2, 3, 4 (2275.25, 2414.27, 1973.13 A), yield 
the strongest Re u lines, but whether they arise 
from an s—p or a d—p transition is still in doubt. 


Zeeman components g 
Inten- | Wavelength Wave- Term Zee- Separa- __ RAS Ay ” 
sity air number combination man tion 
spark type Strongest Strongest Ist 2d 
Pp n 
] 2 3 j 5 6 7 8 9 10 
A K 
5h 5418. 37 18450. 60 3699 2 SF's 
5h 5286. 68 18910. 21 a 3G4—z *P3 
lh 5066. 60 19731. 60 37313; TDs es 
tou 1909. 71 20362. 12 29072— z 5Pj 
10cu 1904. 33 20384. 46 a 5'P3—z *P3 
iO: 1853. 73 90596. 96 a 5P,—z ?P3$ 
iGcu 1703. 70 21253. 91 a 3Fy—2z ?P 
100¢ 1673. 15 21392. 85 a *F,—2z *P3 
8 1664. O1 21434. 78 3493,.— 2: 5D} 
30 1584. 49 21806. 57 a 5P,—z?P 
hou 1542. 89 22006. 25 3738,—z *P3 
20 1520. 97 22112. 94 3493, iD = 
6 1502. 98 22201. 29 3493.—2z 5D5 
100¢ 1481. 32 22308. 59 a °F z *P3 
y/ 1475. 73 22336. 46 3699, ‘1p 
2h 1473. 60 22347. 09 3731 59663 
10 $452. 68 22452. O8 3022;— 2 ®P3 
th 1429. 76 22568. 25 3287.—z 5F3 ‘ 
9 92 49 99809 6 a *F,—z 'P3 
20 1423. 0 22602. 64 1 2699 ; 52 
3 1409. 56 22671. 63 26762 5Pi 
10cwl 1389. 60 22774. 72 2666,—z ®Pj 3 0. 962 | 0. 962w 1. 462, 2. 424 1. 462 2. 424 
2 1384. 44 22801. 52 3225 7DS ios 
LOcw! 1381. 00 22819. 43 3098,—z ®P3 
65 1356. 29 22948. 86 2972.—z 5P3 
2/ 1335. 23 23060. 34 3398; ‘D3 
i) 1330. 69 23084. 52 30713;—z ®P3 
20h 1311. 68 23186. 29 3225;—z 5F3 l Ow 1. 105 (1. 114) 1. 09 
4h 1248. 27 23532. 37 3606;—z ®F3 - . waa 
tow 4240. 18 23577. 27 3022;—z ®P3 . 
2h 1186. 26 23880. 95 33164—z 7D3 7 ne ' 
5Hw 1149. 46 24092. 73 32954—z 7D3 , Ram ‘ 
2h 4135. 43 24174. 46 3287,.—z 7D3 ; a = 
20h 4091. 96 24431. 28 3101,—z ®F3 3 .340 0. 670 1. 078, 1. 418 1. 418 1.078 | 
| 2h 4088. 69T 24450. 83 3493,—2z °P3 cuem L occeuw | tduane | eenwewewaeeeen) Be : . 
1 3 4085. 14 24472. 06 32954—z °F} : ‘ 








TABLE 1, (¢ lassified lines of the second spectrum of rhenium (Re u Continued 








Zeeman components q 
Inten- | Wavelength Wave- Term Zee-  Separa- . 
sity air number combination man tion In 
spark type Strongest Strongest Ist 2d s 
p Mi sp 
l 2 3 } 5 6 7 Ss ) 10 
A AK 
2h O51. 45 24675. 55 3738,—6 2053 
3 1043. 10 24726. 51 3071 KF 2 0.180 O 1. 624 1. 264 1. O84 
3 1042. 78 24728. 47 3493. —5 966 7a, | Oh 1. 132 1. 23 1.19 
lh 1041. 11 24738. 69 37315—6205 
5h 1032. 35 24792. 42 3995 D 
20¢6 1031. 42 24798. 15 a °*G P 
Dh 1025. 56 24834. 24 3022 D \ 
bh 1024. 040 24843. 63 3220) 3) | 
10h 1020. 860 24863. 27 3606 I 3 0.193h 1. 178w 1. 21 1. 16 30) 
th 1018. O16 24880. RS 3995 D 
LOA 1009. 94 24930. 98 277 P l 616 | O 0 1. 232 1. 848 
3Hw 3988. 63 25064. 18 3699 6205 
10h 3984. 28 25091. 54 a XG P 
3Hw 3971. 63 25171. 50 3225 7 . 
2Hw 3964. 11 25219. 21 322 I . 
) 
3e 3963. 53 95222 95 3220, |: 
5h 3939. 381 25377. 53 2977 DD 7a 2.272 | O 2. 272 0/0 2. 279 
LOA 3915. 382 25533. O08 3738 6291 2 0.166 O 1. 570 1. O72 1. 238 
Qewl 3900, 24 25632. 18 9942 D } 
Dh 3893. 59 25675. 97 3398 966 
7CWSs 3888, 92 25706. 80 2809 P 
5h 3886. 21 25724. 72 3731;—630 | 
Ow 3873. 501 25809. 13 3731 63123 
5h 8858, 528 25909. 26 2676 P 
Qh 3856. 72 25021. 42 3699 G291 { 
1OHu 3847. 72 25982. 05 2007 
10h 3839. 540 26037. 40 3101 DD 0 0. 572 1. 42 1. 14 
20cew 3836. 786 26056. 09 QT74 P 
2h 3832. 92 26082. 38 3738 6346 
SA 3830. 525 26098. 69 a 5P P | 0. 758 0 912 1. 670 2. 42% : 
Dh 3826. 55 26125. 79 S101 1) 
2Hw 3825. 30 26134. 32 3699 63123 | : 
2 3821. 43 26160. 80 3316, D : 
100cw 3800, 905 26301. 75 aq P 2 1.137 | O t. 055 0. 644 1. 781 | ss 
2 3791. 57 26366. 81 2007 I 
Sh 3782. 95 26426. 89 3316, Fs 
Qh 3776. 70 26470. 62 3609 A346 , 
l5cw 3773. O18 26496. 5t 33164—45966 2 Ow 1. 185h 1. 188 1. 189 , 
lh 3758. 51 26598. 73 2977 I) ss 
th 3752. 15 26643. 82 2972 z 5D 10) 
200h 3742. 257 26714. 26 a 5s P 3 0.463 0. 927 1. 714, 2. 177 1.714 2.177 \ 
5Hew  3732.058 | 26787. 24 a 9Fy—z *P; | 
lOHew! 3731. 662 26790. 10 3287 5966 
5He 3730. 89 26795. 64 3606 Ds; | , 
LUA 3714. 51 26913. 80 3022 ) 2 Ow 1. 525 1. 398 1. 33 si 
5h 3710. 22 26944. 91 2942 D \ 
30h 3698. 565 27029. 82 3738 f / 3 260 0. 522 0. S10. 1. O70 1. O70 0. 810 
20h 3697. 932 27034. 45 2676 P 2 555 0 2. 754 1. OSO 1. 644 
10h 3692. 783 27072. 14 3225 3 0. 063 1. 123 1. 114 1. 135 , 
Sh 36585. 803 27123. 41 3220 l Ow 0. 928 1. O8 1. 13 
Ll5cew 3674. 87 27204. 10 3022 F 3 325 0. O75 1. O56, 1. 38] 1. 38] 1. 056 
7h 3668. 496 27251. 37 3234 7 Oh 1. 125h 1. 15 1. 16 , 
3h 3660. 17 27313. 36 2774 Ds 
7h 3656. 834 27338. 27 3225 F; 2 Ow 1. 314 1. 114 1. 15 5 
10hc 3647. 544 27407. 90 3225 5988 
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spectrum of rhenium (Re 1) 


Zeeman components 


Separa- 
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Strongest 
Pp 
6 7 
0 
Ow 
Ow 
Q. 234 
0. SO4 SO3 
Ow 
0.700 O 
45 0. 730 
0 
207 0 
213 0 


660 . 2e7 
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0. 616 l 


Strongest 
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1. 044 
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TABLE 1. Classified lines of the second spectrum of rhenium (Re 11)—Continued 
—=—= 
Zeeman components 9 
Inten- | Wavelength Wave-" Term Zee- | Separa- Inten 
sity air number combination man tion —_ sity 
spark type Strongest Strongest Ist 2d spar! 
p n 
l 2 3 } 5 6 7 8 9 10 1 
A K 
10¢ 3299. 790 30296. 24 a 5P;—2z 7D3 1 386 0 0. 630 1. 403 1. 790 2c 
2h 3298. 535 30307. 76 3606;—66373 5 
5ew 3298. 184 303 10. 98 2907 »—z °P$ 401 
3 3286. 636 30417. 46 3398,—64403 3c 
Th 3268. 278 30588. 34 { 2907-09663 ” 
J235/ 9 65405 
l 
3h 3263. 297 30635. 02 3493.—65573 7a, 2 0 1. 251 1. 23) 1. 24 6 
5h 3261. 424 30652. 62 3220,—2z 5D} 3 0.206 0. 821 1. O85, 1. 291 1. O85 1. 29] 10 
Oh 3254. 936 30713. 71 2977 9—6048} 7a 1.128 | O 1. 128 0/0 1. 128 . 
l 3252. 482 30736. 89 3493.— 65673 3c 
2 3251. 548t 30745. 71 a 5P,;—z 5D§ y 
lo 
lh 3250. 142 30759. 02 2972,—6048 15 
80 3246. 311 30795. 3 a *D,—z7P 2 0.263 0 2. 318 1. 528 1. 79) i 
15 3244. 464 30812. 84 a*G,—z7D 7a, 2 0 1. 101 I. 12 LU 15 
lOcw «3241. 910 30837. 12 3220 ,—630/ ar 
l 3238. 470 30869. 88 3225 63125 
rh 
2H 3229. 804 30952. 70 3398.—649/ | 
2 3227. 76 30972. 30 3738 6835 
1OAl 3226. 37 30985. 64 3606.—6704 37 ISl 0. 881 0. 96, 1. 14 0. 96 l. 14 on 
3 3220. 258 31044. 45 3101.—6205 of 
5Oc 3219. 906 31047. 84 a 3} g 5P I 559 OO 145 1. 263 1. 82] ’ 
5} 3218. 653 31059. 93 2942, 8048 “ 
bh 3213. 168 S1112. 95 316 >'p y 2 
i 3203. 560 31206. 26 5225. —B346 
3 3198. 330 31257. 28 a 3] 2 5P 1Y 
2h 3197. 59 31264. 52 36005-67324 . 
15 3195. 206 31287. 84 2963 ,—z 5] 2 Ow 1. 365 1. 14 1. 18 : 
3192. 852 31310. 91 oS0br—0624 2 
diol BSA2 ~ 
2} 3175. 094 31486. 02 3738 ASSA “ 
2} 3173. 604 31500. 80 2809 ,—z 5 F3 - 
7} 3170. O74 81535. 88 3287 é i41 3 Q. 247 0. 498 0. SL3. 1. O6O 1. O6O 0. 813 
2144 2 7D; s 
8} 3166. 609 31570. 39 2809, 5966 7a, 2 0 1. 276h L. 2 1. 20 2 
LH 3155. 209 31684. 5 3398.—656 8 
| 3151. 498 31721. 8 a 5P z 5F ( 
lh 3143. 974 31797. 7 3316,—6496 3 
2 3142. 018 31817. 5 a §G;—2z 5D; 
6 3138. 796 31850. 2 277 2 SF; 7a, 0 1. 120 1. 23 1. 20 f 
lh 3138. 663 31851. 5 3731;—69173 d 
low 3136. 178 31876. 7 3098,—z 5D3 Af 
20 3131. 952 31919. 8 2774 59663 7.3 0 1. 200 1. 23 1.18 2 
Th 3130. 203 31937. 6 323 27D: | 
2 3127. 521 31965. 0 3731;—69283 10 
2h 3123. 622 32004. 8 2314,;—2z 7D} . | 5 
2 3118. 002 32062. 6 3234;—64403 5 
3c 3117. 652 32066. 2 3287.— 64943 2 0 1. 920w 1. 08 1. 36 8 
2 3116. 720 32075. 7 3220,—z 7D3 
2h 3113.350 | 32110.5 3398;—66103 a 
| 2h | 3110.326 | 32141.7 3071;—z 5Dj . a 
| 2 | 3109.983 | 32145.2 3098,—63123 1¢ 
6 | 3109. 576 32149. 4 3225;—64403 z 7 : 9 
2h | 3109. 156 32153. 8 3225;—64413 ; - ; 


3c 3107. 298 | 32173.0 a *F;—z °P3 -- || 
2h | 3106.948 | 321766 3699;—69173 la 
7h | 3105.082 | 32196.0 3606;—68253 a Bt 
20¢ 3103. 260 | 3221483 } osc sp | 
| 20 | 3103.055 | 32216.96 |f % S22 ‘Pi ~ Bs : 
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1. 790 


S13 


2") 


a”) 











Inten- 
sity 
spark 


2ew 

lo 
cw 
3 

10 


“Ito WwW Or 


50 
2ew 

l0cew 
2h 


10c 
2h 

80cw 
8c 

10c 


TABLE 1. Classified lines of the second spectrum of rhenium (Re 1)—Continued 


Wavelength 
air 


A 
3099. 064 
3096. 312 
3095. 314 
3092. 607 
3092. 288 


3091. 356 
3087. 854 
3080. 660 
3080. 048 
3074 654 


3066. 020 
3063. 762 
3061. 930 
3058. 149 
3057 


3048. 036 
3046 644 
3044, 969 


3039. 804 
3038. 812 
3038. 694 
3038. 368 
3035. 776 


Sf 26 555 
3025. 36 
3022. 194 
3012. 36 
3010. O77 


3009, 362 
3008. O84 
S008, 826 
3007. 62 


3002. 6337 


3000, 828 
2999. 947 
2998. 700 
2996. 842 
2996. 182 


2996. O74 
2992. 746 
2990. 926 
2990. 604 
2990, 433 


2985. 207 
2980. 90 

2977. 148 
2974. 563 
2971. 734 


2967. 415 
2967. O18 
2957. 907 
2955. 730 
2954. 99 


Wave- 


number 


K 
32258. 5 
32287. 1 
32297. 6 
32325. 8 
32329. 2 


32338. 9 
32375. 6 
32451. 2 
32457. 6 
32514. 5 


32606. 1 
32630. 2 
32649. 7 
32690. | 


Ww 


S2798. 5 
32813. 5 
32831. 5 
32887. 3 
32898. | 


32899, 3 
32902. 9 
32931. 0 


33031. 3 
33044. 3 
33079. 0 
33 186. 9 
) 
) 


33202. ‘ 


33220. 0 
33224, 2 
33225. 9 
33239, 2 
33294. 4 
33314. 5 
33324. 2 
33338. 1 
33358. 8 
33366. I 
33367. 3 
33404. 4 
33424. 7 
33428. 3 
33430. 3 
33488. 7 
33537. 1 
33579. 4 
33608. 1 


33640. 6 


= 


33689. 
33694. 
33797. 
33822. 
33831. 


wn~Ime-< 


Term 
combination 


3398 66243 
a 5P;—2z 7D3 
2314,—z 5F3 
3071 63043 
2972, 6 2055 


E727 


3493, 


(xe) 
a 5P z §D3 
3101.—63463 
3731 69774 
a ®F, 7D 


a iF, 2 7D; 
2942,—6 205 
a oP, 2g 5D; 
3022;—62913 
3287 »—655? 


2666,—z °P3 
32204 
3699;—697 2 
3287.—65673 
3606;—688? 


2809,—z ®F% 


2762;—6052 
267625966 
a 3F,—z 5F$ 
3022 63123 
33164 6610 


a "(a4 Z 3 
33 164—6624 

2972 6291 

3316 HAS ri 
296234 f S53 
3606:—69283 
3287 »—6610 
2314,—2z 5D; 
3022 G3465 
a *( Te g i] 3 
3225 65457 


3398 67313 
3398, 67323 
a C4 ra SFY 
3220,—65573 


33164 66532 
29634 —63043 
2nOe yy 

S098, O44U5 


a 3G4—5966 


$ 


a 3F, z7D$ 
2963, 63123 
a 4 ¥3 2 5P3 
32954— 66533 
a 8Gs—z ®Fj 
2907, 62717 
30713;—64403 
3071;—64413 
a 5P,—z 5D3 
2666;,—6048? 
a 5Go 2 5P3 


Zee- 
man 
type 


~ 


447 


Separa- 
tion 


0. 282 


1. OOS 
0. 194 
. 619 
. 150 


187 


0. 710 


. 363 


Zeeman components 


Strongest Strongest 
Pp n 


7 8 


0. 835 1. 138, 1. 420 


0. 148 1. 459w 

0 0. 138 

0 980 

0. 763 1. 284, 1. 903 
$52 1. 246, 1. 396 

Ow 1. 027 

0 0. 921 

0 2. O27h 

0. 168h 1. 104, 1. 291 

Ow 1. 168 

Ow 0. 962 

0 1. 096, 1. 812 

0 1. 221 

0. 725 1. 312, 1. 675 


1. 420 


1. 422 


1. 206 
1. 368 
1. 903 
1. 396 
1. 22) 
1. 40 

1. 65 

1. 104 
(1. 398) 
(1. 12) 


1. 110 
(1. 22) 


1. 138 


1. 496 


2. 274 
. 174 
. 284 
. 246 


1. 18 
1. 28 
1. 27 
1. 291 
1. 51 








TARI gz i. ( lassified nes of the second Spectrum of rhe nium (Re Il Continued 








Zeeman components g 
Inten- Wavelength Wave- Term Zee-  Separa- 
sity air number combination man tion . It 
| spark type Strongest Strongest Ist 2d 
| p n s] 
l 2 3 } 5 6 7 S 9 10 
A K 
Zh 2954. 454 33837. 3 2907.—6291 
3h 2952. 982 33854. 2 3287 6673 
3h 2950. 731 33880. 0 33164—6704 
Sh 2943. 794 33959. 9 3098,;—6494 | 0.158 O 0. 72 1. 197 1. 355 | / 
3h 2941. 141 33990. 5 3225 6624 
50c/ 2940. 974 33992. 4 2314 2 2 0. 19 0 1. 463ew 1. OS 1. 27 
3 2939. 557 34008. 8 a 3G 2 
10 2934. 529 34067. | a 3G z 3P3 l 620 0 0. 900, 1. 519 0. 899 1. 519 
Th 2933. 542 34078. 5 3731 1395 
6h 2932. 657 34088. 8 a 5P z 5} 
lh 2929. 044 34130. 8 3731 714 
5h 2927. 852 34144. 8 3516,—6731 
20h 2926. 661 34158. 7 3316 67323 
10/ 2926. 13 34164. 8 3220 ,— 6637} 7* 0 lL. OF | 
10/ 2925. 703 34169, 9 323 6651 ; 0. 164 1. 120 l. 11 l. 14 
6h 2922. 694 34205. 0 3493 A914 | . 147 0 0. 789 1. 230 1. O83 
lh 2917. 173 34269. 8 5220 G6475 
20 2916. 731 34275. 0 a 34q 2 iF 2 267 60 1. 974 0. 910 1. 177 
2cwl POL5 645 34287. 7 a Cy BO52 
5 2913. 639 S4311. 3 277 6205 7 0 1. 227 
3) 2912. O25 34330. 3 3220 6653 
the 2910. 818 34344. 6 a 3G 5966 
the 2909. 127 S4364. 5 33408 AS3S5 
2 2907. 23 34387. 0 2007 68346 
2h 2906. 11 $4400. 2 3287 6727 
eu 2905, 854 34403. 2 S699 7139 
6/ 2904. OS2 S4413. 5 3731 ‘173 
6 2902. 886 34438. 4 3287 6731 
6 29000. O12 54472. 5 3225 667 33 ? 0.166 O 0. 776 1. 108 1. 274 
1Oc 2897. 242 34505. 5 a 3] z‘D 7? Oh 1. O79h 
3/ 2892. 724 34559. 4 3101 655 
LO 2892. 093 34566. 9 a 5p 2 5f—) 3 273 =O. 826 1. 140, 1. 413 1. 413 1. 14 
LOcu 2880. 832 34594. 0 a 3] 
10k YPRSS. O6 34615. 3 a *D 2 5P ; OOT O09 1. 502. 2. 406 1. 500 2. 407 
LO/ 2886. 774 34630. 6 731 19 
2) 2885. 710 54645. 4 2963 z7]) 
20/ 2884. 244 54661. 0 S10] ASF 2 209 1 2O5 | 
2) 2882. 367 34683. 6 2972,— 644063 
60) 2881. S74 34689. 5 a 3G, | » Ow 1. 390 i. 1. 18 
S/ 2880. 630 34704. 4 3234 APO 
3/ 2877. 528 34741. 9 22 aAGY, 
30 2875. 698 34764. 0 2809, Ds; 3 0. 2] 1. 26w 1. 23 1. 28 
2 287 1 32 S4815. 6 77 13 
2 2871. 102 34814. 6 2800 #29) 
7/ 2870. 000 34833. 0 a 5p Is 0.759 O 0. 445 1. 409 1. 168 
be 2869. 200 34842. 7 ,99() BOS | 
2H 2866. 204 34878. | 5398. B8x¢ | 
3/ 2865, 554 54887. 1 S208 ARS 
10/ 2864. 280 34902. 6 a 5p NAA, 
10¢ 2861. 260 4930. 4 aw > 5] 7? Ow 1. O28 
20/ 2860. 550 S4948, |] 2809, —630/ | Ow 0. SSO l. 22 1. 33 | 
2 2859. 92 34955. 8 3699 19 
lL5cu 2859. 87 34956. 4 1 54 P 2 998 0 5 634 0. 640 1. 638 
2/ 2857. 662 S49083. 4 5234 fi ‘ 
S/ 2857. 575 S4084, 5 2942 6441 
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168 


a 
— 

















Inten- 
sity 


spark 


Rh 


6 
10h 
Th 


HOcwl 


oe 


20h 
10/ 
10/ 
20% 


ST, 
ol 

lh 
20/ 


of 


10/ 


50« 


lh 
3/ 
1/ 
2 


10/ 
L00cws 
6/ 
80¢ 
20d 
30% 
2} 
300¢/ 


beu 


150¢ 


Wavelengtl 


air 


2851. 720 


2850 
2849 
2848 
2847 


2844 
2843 
2842 
2841 
2841 


2840 
283 

28: 
28: 
28: 


2832 
283 

2831 
2SR25 


2824 


2824 
IRI9 
2821 
2819 
2815 


2813 
2809 
2803 
2803 


2802 


2798 
2798 


2797 


ad 


2794 


2789) 


590 
99 

966 
568 


S54 
762 
58S 
669 
> le A 


azivd 


O24 
782 
O5 

857 
540 


245 


596 
341 
110 
779 


31 








O04 
3OR 
562 


911 


TABLE 1. 


A 


35056. 3 


35070 
35087 
35090 

35107 


35140 
35154 
35168 
35180 
35185 


35200 
35265 
35287 
35289 


352903 
35207 
35302 


35303 


35381. 


35388 
35392 
35421 
35436 
35453 
,5509 
35532 
35583 
35653 
35662 
35666 
5) 
35718 
35739 
35772 
} } 


! 


_ 
é 
~ 
‘ 
- 
‘ 
or 
Star 


Hos 


35933 
3590890 
36114 
36146 
36151 


36214 
56242 
56251 
36263 


36275 


36304 

3633 
36345 

36372 


36459 


36559 
36578 
36587 
36598 
36607 


Wave- 
number 


ee ie Boe | 


Nee D> 


“N s3I St to 


“St + 


( lassified lines 


Term 
combination 


3493 
3098, 
2762 
2676 
3398 


2007 
2963, 


a 5P 


SOUS, 
3295, 
SOUR, 
2972 
3493 


3316, 
aes 
a Cis 
aG 


S316, 


of the second spectrum of rhenium (Re 11) 


69283 


6304 

AS8252 
B494 
63123 
27D; 


663; 4 


2 


6651 7: 
65013 
463 7 j 
> 5p 


2 3p 


“ro ue 


6673 
6346 
6972 


5966 


Li) 


65) 4 


6501 
G55?) 
2°D l 
69285 
6637 


7108 


‘ 


G6475 
672 
ASO 


AOISO 
6732 
4610 

139 


‘ 


697 2 


Zee- 
man 
type 


Continued 


Zeeman components 
Separa- 


tion 
Strongest St rongest 
p n 
6 7 8 
10 0. 31 
0 1.513 
0 1. 119 
0. 27 0 0. 901 
Ow 1. 690 
0 lL. i3¢ 
0.570 0. 574 0. 797, 1. 365 
. 045 0 2. 522 
2.452 | 0 2. 452¢ 
Ow 1. 527 
0. 54 0 0. 05 
0. 402 1. 141 
Ow 1. 595w 
0. 436 0. 873 0. 650, 1. O86 


(1. 188) 117 


1. 188) 1. 18 
1. 366 0. 796 
1. 142 1. 487 


0/0 2. 


1. 191 


1. 68 1. 14 


Ck. 7a) 1. 67 


0. 650 1. O86 








Inten- 
sity 
spark 


10h 

10 
2h 

15h 
3el 


3h 


10cw 


“ac 


3h 
L5c 
Sh 
6h 
10h 


1Ocew 
th 
2h 


lh 


5h 
6h 


10/ 
5OA 


Wavelength 
air 


to 


2730. 
2730. 
2729. 
2726. 
2726. ; 


2724. 
2718. 
2717. 
2716. 
2716. 


dry rnwrdy 
S-1-) 


Ss 


2696. 2 


2695. 


2693 
2691 
2691 


2687. 592 


2683 


2682. S88 


2681 


2677. 2: 
2669. 3: 


2667 


2667. ; 


2665 


2664. 556 


2663 


2662. 
2660. 9: 
2659. 5! 
2659. 2: 


2655 


2653. 
2652. 
2651. 


2650. 


2648. 462 


— mm DOS SO 


Wave- 


number 


K 


36610. - 
36618. 
36621. 
36661. < 
36665. 


36696. : 
36772. 
36785. : 
36799. ! 
36806. 


36837. 
36843. 
36853. 5 
36871. | 


36875. 


36884. 
36890. 
36912. : 
36938. 


37002 


37004. 
37023. ‘ 
37025. 
37033. ! 
37076. ‘ 
37089. 


37120. § 
37140. : 


> 
3 


37589. : 
37593. 


37652. 


37675. 
37689. 
37701. 
37721. } 


37746. ! 


9 
9 


; o> My 
Iwo 


ae 


o> “10 0 


See aie | 


— 


Term 


combination 


3071; 
3225; 
3225; 
a *G4 
2774; 
3220, 
2774s; 


2676: 
32954 
a aK; 
2666, 
a * 14 


2963.4 
a*F;, 
3287 
2804, 
29634 


2963, 
2963, 
a 3F 
a 4G; 
a‘D 


67325 
68863 
68874 
z5Dj 
64403 


-68863 


64413 


63463 
697 23 
25D3 
63463 
63043 


66475 
2 3Ps 
697 23 
64963 
66512 


691 / 

63125 
691 7? > 
6928; 


66 733 


74325 
66105 
69285 
z 5Fs 


69 255 


6731 
5966 
6312; 
6048; 


72133 


463 
70133 
B825:% 
6048; 
69794 


49) 2 
6557 
6727} 


69775 


697997 


672? 
AY 7 7 4 
69794 
62913 
z'’D; 


667 33 


67313 
67323 
z'D3 
63043 


> 5Ps 


Zee- 
man 
type 


to 


“J 


to- 


7b 
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TaBLe 1. Classified lines of the second spectrum of rhenium (Re 11) -Continued 


Zeeman components 


Strongest Strongest 
Pp n 
7 Ss 
0 1. 053 
0 1. 227 
Ow 1. 126 
0 0. 738 
0 1. 166 
Ow 1. 062 
0 0. 515 
0 1. 292 





(1. 114) 1. 


(1. 23) 1. 23 


(1. 20) 1. 


1.114 l 
1. 20 
1. 26 
1. O76 





Inte! 
sity 
spar 











1. 30 


19 
31 


637 








Inten- | Wavelength 


sity 
| spark 


| 


air 


to 


SS 


20 

5 
30cw 
3 
30 


Th 
150c¢ 
200cwl 

2h 

Ss 


20¢ 

30? 
8 

100h 
3 


9 


0 
2h 
10h s 
4h 
150cwl 


5h 
Sc 
20h 
2h 


9] 


400 
10 
Sh 
60ed 
l 


} 


7ew 


10cWs 


4h 
L5hAl 


th 
20 

lh 
30 

5h 


300cew 
3h 
2h 
3h 

300 


4hil 
60 
10 
4 


5h 


15 

20 
200 
100cew 
100cew 


2647. 
2644. 
2642. 
2641. 
2641. 


2638. 
2637. 
2635. 
2634. 
2633. 


2632. 
2630. 
2630. 
2628. 
2622. 


2621. 
2620. 
2619. 
2619. 
2616. 7 
2615. ! 
2610. 
2610. £ 


2610. 
2609. 


2608. 
2606. 
2593. 
2588. 
2586. 


2583. 
2582. 
2582. 
2581. 
2580. 


2580. 
2576. 


QnvTK& 


) 
) 
“ 
2975. 
) 
) 


102 
238 
664 
317 
117 


. 810 
. 186 
. 828 
. 192 
. 638 


». 92 
». 366 


35. 582 


2561. 


ORL 


2559. 


2557. 
2557. 
2554. 
2553. 


2550. 


920 
478 


694 
414 
630 
588 
091 











Zeeman components 














Wave- Term Zee- | Separa- 
number combination man tion | 
type Strongest | Strongest Ist 2d 
| Pp | n 
3 4 | & | 6 7 | 8 9 10 
| | 
K | | | | 
37759. 2 a*F;—z*P3 1 0. 263 | 0 | 0. 731 1. 257 1. 520 
37806. 4 a *G;—6312} 2 . 36 | 0 eS ee (0. 90) 1. 26 
37835. 4 a °S,—z D3 3 ee . | Ow So ae (1. 71) 1. 83 
37848. 9 2942,—67 27} | ey ea | stati mriedi etal a eet 
37852. 9 a 5D,—z 5P3 2 326 | 0 2. 148 1. 496 1. 822 
} | 
37894. 1 3098 ,-—6887; 3 . | 0.342 | 1. 22w 1.18 | 12 
37970. 5 a 5G;—z 5F3 7,3 | Ow 1. 136 (1.12) | L15 
37927. 4 a 5S,—z 7D} I . 564 | 0 | 1. 143 1.707 | 2.271 
37942. 3 a ®F,—6048} | cameliiea eae 
37967. 0 a ®Fy—z 5Fy l | Ow 0. 945 | (1. 26) 1. 18 
37980. 3 a *F,— 60523 | fF eee) 
37999. 1 a §G;—z 5D} l 0.199 | 0 | 0. 692 | 1.090 1. 289 
38005. 2 2809,—66 103 si oe oe i 
38036. 7 a 5P;—z 3P3 3 0. 220w | 1. 198w 1. 234 1. 161 
38123. 5 3699;—76113 : hints 
38129. 3 3022;—68353 4 
38142. 9 a *G;—63463 ee 
38159. 4 te ee ee eee eee, eee 
38169. 9 a §*Gy—64403 : ; 
38204. 4 a 5G,—z 7D3 2 .66 | 0 1. 131 0. 65 1. 131 
38220. 9 a 58.—z ®F3 
20900 a 5Gs v4 5F3 ace 
38289. 9 3101,—69303 : 
38294. 8 a *G;—z 7De 
38301. 0 3098,—6928} 
38305. 1 a 5G;—z 7D; 
38324. 8 a 5Dy—z 5P3 3 . 299 | 0. 595 1. 526, 1. 825 1. 526 1. 825 
38354. 5 2774;—66103 3 . 18 . 530 (1. 23) 1. 41 
38553. 1 32253;—70813 3 .116h | 1. 154h 1. 13 1. 17 
38619. 7 a 5G;—z 5Fj : 
38651. 9 3022;—6887 ; 
38700. 9 a 5P;—63463 
38705. 1 a °*G.—64943 c : 
38716. 4 a 5P,—62053 
38730. 6 2774;—66473 ‘ 
38743. 4 2762;—66373 . 
38746. 6 3098 ,—697 23 on es 
38804. 7 26762.—65573 2 Ow 1. 494 (1. 09) 1. 23 
38813. 3 3098,—69793 ; : 
38818. 5 a *G,—z "Di vil ; imaaiins 
38866. 8 a 5G,—z 5D3 ad 
38871. 5 a 5‘D,;—z 5P3 3 0.21 | 0. 64w 1. 56w 1. 46 1. 67 
38896. 1 a 5Ge—z 5 Fs . — ae See 
38901. 5 2972,—68623 . ae ee | meade 
38911. 1 2314,—62053 | . ' toveknnneend T epemebne if sauna 
38919. 5 a 5D,—z ®P3 l 157 | O | 0. 925 1. 479 1. 636 
38945. 5 3022;—69173 " . oa ee eee eo 
38953. 9 2809,—67043 2; | | Ow 1. 550 | (1.22) | 1.29 
38965. 9 a *F,—z 5D; 7b 0 1. 283 1. 28 1. 28 
39021. 5 a *F,—62913 oe oaiiditekeanis | Se ee "7 
ac 3022;—6928} Za, 1 0 1. 220 | (1. 398) 1. 3 
39058. 8 3071;—69773 | - : BE, ett AE Spesteerg San 
39086. 0 a *G;—64403 2 0.313 | 0 2. 146 0. 894 | 1. 207 
39090. 3 a §G;—64413 2 eonncn | Oe 1. 050 (0. 90) | 0.82 
39132. 9 a 5Gy—z 5FY 7b Ow 1. 183 | 1.18 1, 18 
39148. 8 a 5S.—z 5D? 2 eee geet oe a ee) | Pte 
39202. 5 a 5G,—59663 2 : Ow 1. 207 | (1. 18) 1. 17 





2 ee 





























TABLE 1, ( assified nes of the second spechi em of rhenium (Ret Continued 
’ 
] Zeeman components q 
| Inten- Wavelength Wave- Term Zee- Separa- 
sity air number combination man tion Int 
spark type Strongest Strongest SI 
p , Sp 
l 2 3 5 6 7 S 
, A { 
Sh 2548. 386 39228. 7 2963 ,— 6886 | ” 
7 2548. 020 39234. 4 a 3} 6312% 7b 0 1 262 - 
2 2546. 450 39258. 6 S287 (213 ” 
15/ 2545. 258 39277. 0 2907 A835 , 
1/ 2541. GSI 39332. 7 2676 6610 10 
6 2541. 514 39334. 8 aitG B55 ' 
3 2536. 398 39414. 2 O72 6914 > 
SO 2534. 099 39449. 9 a*D 2 5p » Ow 1. 900 i 
5/ 2532 536 $9474 2 R25 ] ; 
Tow 532. 33 39477. 5 2338 2 523 
3 9526. 536 39568. 0 2774 6731 10 
10 2526. 32 39571. 4 2314 6271 ; 0. 43 0. 42 l. 11, 1. 55 
t 525. 653 390581. 8 9774 6732 
; : ( ) 2 5) 
10/ 2521. 691 39644. 0 x | 
a | é ,40) 
2) 2520. 658 39660. 3 SIRT 253 ‘ 
- ™ lf 
LlOhw 720 39690. 8 32295 ~19 
2h o9 39702. 2 a 5P,—6304 
Qh 190 39746. 5 2338 63123 
2 t14 39790. 4 2007 »— 6886 
10 250 39808. 8 5738 719 7 0 1. 156 
20cu 2510. 11 39826. 8 a I) 
2} 2509. 50 SOR836. 6 007 O81] 
SY, 2507. S54 BORG. 7 a *qy 8A10 
10 2506. 999 39876. 3 3225 213 
} 2505. 210 39904. 8 sil 722 
30 2504. 595 39914. 6 aus ID 
200cu 2502. 348 39950. 4 a *G > 5} 2 0. 3910 1. 230w 
2 2501. O05 39971. 2 S287 IS) 
| 24908 566 LOOTO. 9 a °C, Hh? 
3/ 2495. 386 10061. 9 a *G 6052 
7 2495. 141 10065. 8 2666,—6673 
- 3606 13 
2494. 609 10074. 4 oo O15; 
US (C6 
l 2494. 517 10075, 8 3101 108 
sf 2494. 217 10080. 7 3398 107 
2 2493. 597 10090. 6 2963 ,—697 2 
l 2493. 461 10002 2XOT 691 
10 2491. 234 10128. 7 a°*D > 7]) | 0.319 O 0. 838 
| 2490. 91 10133. 9 G,—6637: 
of 2490. 40 10142. | a*G D3 
50 2490, 161 10145. 9 S098, 112 
! 2489. 482 10156. 9 2063 ,—6979 
30 2489, O04 1O164. O 2809, ASSIS | Ow 4 
ti 2488. 104 10179. 1 a >P BAGS 
teu 2486. 46 10205. 7 15S 2 5] 
2 2484 106 1OY38. 4 am Cy He 
Scu 2481. 98 10278. 2 3225 53 
6 2480. 6351 10300, 2 54903 4 
| 2479. 541 10317. 9 2314 A346 
SO 2478. 99 10326. 8 a *I) D 3 0.7774 O. 778 1. 498, 2. 274 
1 0c 2478. 545 10334. 1 S208 32 
10 2476. 332 10370. | a *} 6291 7 0 1. 340 
: — 310 139 
5h 2475. 470 10384, 2 an ‘oe 
51H, R55 
0) 2475. 172 10389. 0 a *] D 2 Ow 1. 443 
| 2472. 78 10428. | a‘ G205 
20 2471. O49 10456. 4 a GQ Fs 7a. 2 0 1. 264 


2d 


10 


20 


274 














Inten- 


silyV 


spark 


10 
10 
15 


10) 


10 


Oe 


b60 





Wavelength 


“it 


2470. GLO 


2469 
2468 
2467 


2467 


24060 
2465 
2462 
2459 


2458 


2455 
2455 
2455 
2450 
2449 


241 

2414 
2412 
2411 
2409 


2408 
2408 
2406 
2405 
2403 


2398 
2398 
23965 
2395 


M95 


356 
ISO 
85 


574 


S10 
SU 
658 
SY 
833 


827 
635 
097 


UOS 


515 

OO7 
148 

60 


392 
201 


S15 
Ob 


905 
QS 
974 


O47 


sad 


QR9 


— 0 ™— to 
) 


DS 
10) 

96 
1] 


TABLE 1 


Wave- 


number 


» 

A 
10463. 6 
10484. | 
10498. 5 
10508. 8 
10513. 4 
10525. 9 
10542. 5 
10594. 2 
10640. 0 
10657. 4 
10707. 2 
10710. 5 
10719. 3 
10798. 6 
tOS11. 9 
10820. 1 
10895. 1 
10901. 8 
10919. 2 
10926. 9 
10936. 3 
f1001. 4 
41014. 4 
11026. 0 
11050. 1 
11061. 7 
11076. 5 
11082. 9 
11091, 2 
$1121. 7 
11127. 1 
flistl. | 
11155. 1 
{1188.1 
11265 ! 
11285. 6 
11337. 2 
11340. 5 
11347. I 
L360. 6 
11373. 6 
11409. 0 
11430. 0 
$1447. 7 
11485. 0 
11498. 6 
11511. 7 
11548. 3 
11586. 9 
11601. 4 
11678. 7 
f16S1. 4 
11719. 7 
11724. 2 
11733. 8 


Term 
combination 


a Cis 
972 
a (iy 
a ) 
2774 
3101 


3398 
204 
a*qG 
2809, 
2514 


Classified lines of the second spectrum of rhenium (Re 11) 


Zeeman components 


Strongest 


Pp 
7 
0 
0 
0 
0. 672 
0 
0 
0. 579 
0 
Ow 
0. 160 
0) 
Oc 
Oc 
Ow 
0 
0 
0 
0 
Ow 


Strongest 


1. 526, 


2. 260 


1. 198, 
1. O95 
1. 14 
l. Zo 
0. 530 
1. 249 
1. 405 
1. 505 
>») > 
— as 
1 300 
0. 750 
1. 231d 


Continued 


~ 








Inten- 
sity 
spark 


30 
5e 
2 
100cw 
l5c 


10 
2 
5e 
3 

20 


10 
10c 


» 


9 
> 


10cw 
3 


10c 


Wavelengt 
air 


i) 


2359. 157 
2358. 022 
2357. 91 
2357. 66 
2355. 423 
2351. 925 
2350. 653 
2349. 94 
2349. 61 
2346. 987 


2345. 51 
2344. 544 
2341. 31 
2340. 344 
2336. 923 
2336. 47 
2333. 75 
2332. 51 
2331. 81 
2330. 878 
2329. 89 
2325. 424 
2324. 45 
2323. 409 
2321. 75 
2317. 889 
2316. 49 


2315. 186 
2314. 209 


2310. 703 


TABLE 1. Classified lines of the second spectrum of rhentum (Re 11)—-Continued 


h Wave- 
number 


K 


41834. 9 
$1857. 7 
41864. 7 
$1882. 7 


£1909. 6 


41967. 4 
41992. 9 
$2004. 1 
$2011. 3 
12019. 8 


2021.8 
$2030. 0 
$2100. 0 
$2119. 4 
12148. 4 
$2151. 2 
$2167. 7 
$2208. 1 
$2221. 4 
$2231. 9 
$2251. 1 
$2272. 0 
12286. 0 
$2314. 4 
$2333. 6 


$2395. 4 
$2397. 5 
$2402. 0 
$2442. 2 
$2505. 3 
$2528. 3 
$2541. 2 
$2547. 2 
12504. & 


$2621. 6 
$2639. 1 
12698. 0 


$2715. 7 


12778. 2 
12786. 4 


12836. 3 
12859. 0 
12871. 9 


12889. 1 
$2907. 2 
12089. 7 
13007. 8 


13027 
13057. 8 


$3129. 4 


$3179. 8 
13198. 0 
13263. 6 


Term 
combination 


a a] 63 ,63 
3220,— 74063 
2314,—65013 
a 5G;—z Dj 
3022 72133 
a 5p 66733 
ai%G 87313 
2972 71733 
2907 .—7 1083 


a ®D,—z 5D} 


a *¢ 14 68 P52 
2774;—69773 
2676:,—68863 
a 5] ) Zz 7] ) = 


2762;-——69773 


a 5P,—65573 
a is Z 7De 
aig 68252 
3493 77223 
a 5D,—z 5D 
a 5P 656) 
2676 6914 
3225 7463 
a C4 2 Dj 
a 3F GB494 
2676 6917 
als 5968? 
2666 69175 
2314;—656? 
s101 ) 355 
1) z 5F 
D 2 *D 
3022 7284 
a *Ciy 68874 
a *D 27D 
2809, 7081 
a 3] 64405 
a°*D z 5) 
3071 R55 
3225 ,i11 
2907 194 
a *ty 68875 
P,—6624 
2338 B6375 
a 5P 667 33 
a 3} 655? 
3101 7407 
a 4k y—6567 
2338;—66533 


a 3 4—67043 


29554 76135 
33164 76386 
a S; 6048} 


| 
| 


Zee- | Separa- 
man tion 
type 

5 6 

l . 231 

») 

2 0. 234 
7a, 2 

. 251 

2 343 

l * 


] 0. 528 


454 


Zeeman components 


Strongest 
Pp 


7 
0 0 
Ow l 
0 l 
Ow l 
1. 053 ] 
0 2. 
Ow ] 
Ow 0. 
0 0 
Ow l 
Ow l. 


Strongest 
n 


8 


. SO9 


. 303 


. 766 


. 301h 


. 126, 1. 485 


508 
O08 


SUYU 


110% 


127 


g 
Ist 2d 
9 10 
1. 27 1. 502 
(1. 26) 1. 25 
1. 532 1, 208 
1. 23) 1. 25 
1. 481 1. 130 
1. 479 l. 136 
1. 24 1. 31 
1. 694 1. 166 
1. 26 1. 22 
1. 20 1. 28 

















1. 298 


166 











Inten- 


sit 


100¢ 


15 
10 
10 
60 


20 
10 


500c 


10 
50 
5 
10 
) 


10% 


) 


10 
Sh 


y 
spark 


w 


a 


2308. 
2308. 
2308. < 


2308. 


2307. 


2306. 
2305. 
2305. 
2303. 
2302. 
2301. 8 
2301. 


2298 


2298. 


2296. 


2295. 2 


2295 


2293. 
2290. 
2286. 


2283 
2282. 
2281. 
2281 

2280 
2279 
IOTR 


) 71) 


ae | 


2269 


2269 
2268 
2266 
22°66 
2265 
2265 
2261 
2260 
2258 
2254 


ae 


oe 


Wavelength 


ir 


1 
972 
658 


5 
. 060 
215 


268 
934 
021 
993 


345 


948 


. 685 


. 012 


tl 


748 
S46 
500 
oe) 

106 


Wave- 
number 


13296. 0 
13301. 89 
13305. 8 
$3313. 1 
13329. 0 
$3346. 8 
413353. 0 
$3370. 2 
13389. 5 
$3420. 6 


om 
w 
—_ 
— 
_ 
_— : 
re—n 


£3500. 9 
£3538. 8 
13555. 4 { 
13558. 4 
13584. 7 
$3650. 9 
13719. 4 
13775. 2 
13788. 9 
13819. 1] 
$3825. 4 
13846. 1 
13848. 5 
13866. 3 
13872. 4 
13903. 6 
13915. 3 
13937. 7 
13962. 2 
13988. 0 
£4021. 2 


£4054. 2 


14069, 2 
14109 4 
14113. 5 
$4134. 5 
14160. 8 
$4198. 0 
$4224. 3 
14262. 5 
14348. O 
14366. 0 
14379. 2 
14388. 1 


14399, 2 


14408. 0 


Term 
combination 


a 4G; 69283 
a *S;—60523 
a *Gs—z 7D 
a *F;—6 4943 
3022;—73553 


307 13;—7 4063 


3071 74073 
2676:—70133 
a 5P;—66733 
a 3Fy—67313 
a *Gs—64403 
a 3B y—67323 
aGq—697 23 
3287 »—76363 
2762;—71123 
a *®Gg—69773 
a °®G3—68873 
a °F ,—66103 
a *G4—69793 
23 14,—66733 
27743;—7 1393 
a 5G Zz 5D; 
a °G;—-62913 
a 3Gs5s—69773 
a *¢ 4 2 7D: 
2972; 73553 
3022 74073 
a3q 6917 5 
a ®P3;—68623 
a5G 62? lj 
a °G3—63043 
3071 74633 
a ‘Ss; Zz P 
a 3() 69283 
a 5G G3123 
5316,—7719 
2676 70813 
5164 22 
a °G 629] 
3225 636 
a 5p 68873 
098,—7511 
3220, A36 
a 4 6304 
d101.—74523 
3097 523 
3287 1722 
a 3#F.— 6825 
a 5*P GH914 
277 7213 


TABLE 1. Classified lines of the second spectru 


Zee- 
man 
type 


455 


m of rhenium (Re 1)—Continued 


Separa- 
tion 


Zeeman components 


Strongest 
p 


7 


Oh 


Ow 


Ow 


Ow 


Strongest 
nm 


8 


1. 626 


2. 00 


Ist 


(1. 19) 
1. 25 








Inten- 
sity 
spark 


50 
10 

y 
10 


Ye 
! 
5 
10 


LOO< 


10 
2 
> 


50 
Sc 
20 


50 


1000e 


10 
At) 
1) 
10 


my 
50 
10 


20) 
100 
7) 


TABLE 1 


Wavelength 


air 


2226. 


. 851 


yor } 
2223 


999] 


») i) 


(1) 
2219 
2916 
2916 


2214 


POG 


2208. 


2207 
2206 


2903 


2 ") > 
2202 
2199 
2197 


2197 


2197 
2197 
2195 
2104 


21903 


2192 
21900 
2188 
2187 


2184 


2185 

»1R9 
2181 
2181 
2180 


2178 
21a 


500 

Ov 
169 
S03 


865 
964 
292 


“07 


104 
048 


x”) 
845 
920 


130 
S04 


iia 

851 

579 
SZ 


He2 


Wave - 


number 


14731 
14769 
14783 
14806 
14847 


14908 
14932 
14966 
$4993 
15012 


15017 
15036 
15101 
15109 
15147 


15404 
15499 
15537 
15555 


155 > 


15504 
15642 
15684 
156908 


157690 


15784 
15807 
L5S19 
15827 
15840 


15882 
15902 
15908 
15912 
15969 


sified line 


of the second spect) im of rhenium 
Zeeman components 
Term Zee- Separa- 
combination man tion 
type Strongest Strongest Ist 
p fi 
5 6 7 8 q 
2005 , L065 
a%D z 3p 
a *} 6610 
a 5G y—6494 
a 5G 6651 
a *P 69283 
a5P 8930 
2676 7139 
a‘*D z 5F; 
a 5D) ,—2z 5F 3 
a *F 6727 
a 3} 6731 
a 5]),—4966 
16, 7797 
a *F 66475 
2762 P53 
5225 719 
ap 6972 
2007 AL 
s022 523 
220, 7 22 
a*D P 
a 3] 6673 
a "Cag, B55 
a’s z7P ; 0. 545 1. 884w 1. 97 
1 3] G914 
a ®G;—6441 
a 5*P ASA? 
D 966 
2676 213 
OOS, as 
aC, 66375 
2PSO9, 354 
2514 HSA? 
(; OS] 
P ~ 277 1? 
rts HO475 
y (5 A657 
POODT 163 
(, / A 
I) AOS 
D,—60 
| 6731 
F,—6732 
2676 4) 
5) Q 
(, AAS 
(y I12 
a ~ ASO 
ry, AUF 
2879 Aa 
1°] B9774 
p> «> Se AUIS 
( (; He 
2 6135 
2800 O85 


Continued 


Re u 


456 


2d 


10 








2d 


10 











-“ 











Inten- 
sity 


spark 


1O¢ 
20 
9 


TABLI 


Wavelength 


aw 
? 
| 
2172. 618 
2172. 094 
2171. 145 
2170. 970 
2170. 806 
2169. 77 
2169. 240 
2166. SSO 
2166 52 
2163. 07 
2161. 914 
2158. 292 
2158. O5 
2156. 103 
2153. 79 
2152. 98 
2152. 79 
2149. 597 
2148. 261 
2148. 15 
2145. 90 
2144. OSS 
2143. 96 
2140. 190 
2139. 45 
2136. 701 
2136. 22 
2134. 773 
2133. 776 
2133. 122 
2128. 775 
2128. 53 
2128 45 
2127. 28 
2126 UH 
2126. 73 
2122. OSI 
2121. 27 
2119. 67 
2118. 671 
2118. 54 
2116. 23 
21 15 >I 
2114. 25 
2111. 854 
2102. 28 
2100 RH 
2100. 58 
2097. 77 
2093. 792 
20901. 920 
POOL. 535 
2090. 20 
2087. 67 
2087. 30 


- 


Classified lines of the second spectrum of rhenium (Re 11) 


Wave- 


number 


AK 
16013 


16024. | 
16044. 2 
16047. 9 
$6051. 4 
16073. 4 
16084. 6 
16134. 8 
16142. 5 
16216. O8 
$6240. 7 
16318. 4 
16323. 5 
16365. 4 


16415. 


16432. 
16436. 
16505. 
16534. 


16537. ! 


16585. 8 
16625. 2 
16627. 9 
16710. 1 
16726. | 
16786. 3 
16796. 8 
16828. 6 
16850. 4 
16864. 8 
16960. 05 
16965. 9 
16069. O 
16993. 5 
17000. 6 
17005. 7 
17108. 6 
17126. 7 
17162. 2 
17184. 4 
17187. 4 
$7238. 9 
17255. 0 
17283. 1 
17336. 7 
$7552. 3 
$7584 7 
17590. 7 
17654. 5 
$7747. 3 
17787. 7 
17796. 5 
17827. | 
17884. 9 
17893. 4 


Term 
combination 


a Ci, 
314, 
a 5D, 
a 5p 
a iy 
a 5G, 
a*F, 
a*D, 


3022 
3097 


a Sy» 
5098, 
2774; 
a 3F, 
a *P 
a 5p, 
2338 
a 5G 
762 
3071 
a°G 
2809, 
a °C, 
a *k 
a%®G 
a 1) 
a G 
2676 
a 5P, 
a Cs 
a °C, 
a Ci, 
a (; 
3022 
a 5P 
SOUR, 


a 3] 
2762 
a °G 
a iT) 


a Ss 
a r 
a (; 
a 7 
a G 


206034 
avGW 
a | 

2338 
aos 
a Cay 
a 34 
a%*G 
a ®D, 


A6 444 
691?) 4 
> 5 F 

70818 
66512 


66532 
656 25 
6048} 
76363 
77193 


63 16 9 
YOO 

tawwg 
74065 
68863 
71123 
70135 


69; 4 
B55? 3 
14063 

»*o9o 
6567 5 


‘ 16. ; ; 


AOL 
6417 
71945 
GO485 


HDD 4 


1308 
70133 
6567S 


6610 
7719 
7173 


6441 
OH8SS6 
AIO 
BS25 


66475 


7719 


7) 


457 


Inten- 
sity 


spark 


20 
10 
6Cc 
Oc 
30¢ 


10 
20 
20 
30 

9 
10 


200 


10e 


10 
30¢ 


Wavelength 


air 


1 


2085. 76 
2085. 14 
2083. 67 
2083. 342 
2082. 234 
2081. 26 
2079. 33 
2079. 05 
2078. 33 
2077. 78 
2076. 30 
2075. 70 
2075. 11 
2072. 88 
2072. 62 
2072. 09 
2071. 175 
2069. 42 
2067. 64 
2066. 91 
2065. 83 
2065. 64 
2064. 87 
2064. 16 
2059. 75 
2059. 33 
2058. 520 
2058. 10 
2056. 59 
2055. 25 
2053. 597 
2053. 10 
2050. 64 
2048. 380 
2047. 955 
2047. 67 
2046. 318 
2042. 64 
2040. 64 
2036. 20 
2035. 51 
2034. 824 
2034. 020 
2032. 59 
2032. 02 
2027. 69 
2027. 45 
2027. 29 
2026. 40 
2024. 05 
2023. 644 
2020. 475 


2020. ° 


2019. 


2019. : 


Continued 


Wave- 
number 


K 


48640. ¢ 


18679. 
18691. : 
18749. 


17928. 8 
$7943. 1 
47976. 8 
47984. 5 
$8010. 0 
48032. 4 
418077. 1 
18083. 6 
48100. 1 
$8113. 1 
18147. 3 
18161. 2 
18174. 8 
418226. 7 
48232. 8 
$8245. 1 
48266. 4 
48307. 2 
18348. 9 
18366. 0 
18391. 2 
$8395. 7 
18413. 7 
18430. 4 
$8534. 0 
18544. 0 
$8563. 0 
$8572. 9 
18608. 6 


18803. 3 


ISS15. ! 


{8S82O0. ¢ 
18852. 
18940. 
18988. 
19095. : 


19111. 
19128. 
$9147. 
19182. 
19196. 


~ 


w = & Oo 


—oIeoun 


19301. 3 
19307. 2 
19311. 1 
49332. 5 
19389. 8 
19399. 8 
49477. 3 
49480. 7 
19489. 7 
19506. 7 


Term 
combination 


a *Ds 
2314, 
a 5D, 
a *Ds 
a 5P, 


a ‘Ds, 
a 3G; 
a 5P; 
a 3F, 
a *Ds 


a 3F, 
a 5D, 
a ®QGs 
a 5P, 
a *Qsz 


a‘*D, 
a?’ F, 
a ®G; 
a 5S, 


a 5D, 


2774; 
a 5P, 
a*Qy, 
a *Gy, 
a *Go 
a ‘Ds 


a iF, 
a‘*D», 
2666, 
a 5P, 
a*D, 


a Gy, 
a‘*D, 
a °*Qs 
2314, 
a@ 5P, 


a *Gy, 
a*F, 
a *)) 
2314, 
2809, 


a *Ds, 
2809, 
a *Qs 
2879; 
a 5P, 


a *P; 
a 5P, 
a °QGsz 
a *Qy, 
2314, 


a *D, 
a‘Ds 
2774; 
a*Ds 
a ®QGs 
a 5S» 


25D; 
71083 
25D; 
—f 2913 
71733 


62913 
74063 
72843 
697 23 
63043 


69773 
63043 
6731 3? 
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62913 
15233 
71943 
O3463 
GI143 
G3043 
4073 
71943 
72533 
69283 
71393 
62713 
72133 
77193 


74063 
74073 
74633 
69793 


YOR@ 
/ 2533 
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66738 
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3a 1927. 912 51869. 6 a *G.—?? 








40 2018. 54 49524. 8 a 5D,—6440% 
: 2017. 08 49560. 6 a 5P;—74323 7 ; 193 
: 20 15. 91 49589. 4 a °F,—72133 5a 1926. 93 51896. 0 ” bs 66243 
2013.02 | 49660. 6 2389,—73553 5 aT 1—76233 
. , JO . 20394 (3003 oa 1926. 566 K100n 7 5 Bese 
10 2011. 82 | 49706. 2 a 5P ~ L1G - 1° ). ++ 51905. 7 a °D,—66733 | 
} 4403 p 926. 096 51918. 5 a 58,—69143 
5a 1925. 05 51946. 7 sass 
100 2009. 92 | 49737. 1 a 5G;—685 91946. 7 a ®5S;—6917 
. 2008.67 | 4976 a. wanes . asa 
- - - | 49768. 1 a *F.—7 1393 oa 1913. 60 52957. 5 5c . | 
6 2004. 14 49880. 6 2809 ~vOrs “vt. 0 a °Gs3 ‘ 1393 
| 677975 12a 1908. 98 523840 ({ @*Ds—6731; 
1 a 5Gy—728s 
a 10a 1908. 48 52397. 7 a *D.—47en 
vac 6 ( b> = I 3 67325 
ba 1907. 385 52427. 8 a 5P.— 77193 
On e 7 oa 907 R942 F em laiagd._ 
25 1999. 50 50012. 5 a 5D; 64943 ] 07. 10 52435. 6 a 30 . 1063 | 
3a 1998. 77 50030.8 |{ % °G2—6887% 10a 1906. 73 52445. 8 a *D.—6732; 
> Yo ~ . “OMOEA 
| a5G;—69173 15 —_ 49 52452. 4 a *D,—6) 278 
— 2 3a 906. 19 52460. 7 5P.—_ 779% 
7 tae 60 50060. 1 a *D,—64943 5a 1905. 689 52474. 5 ; “4, di. 3 
a 996.663 | 50083. 5 a ®D.—64963 12} P Ee a oe 1 — 76363) 
6. 66 { 83. ! L96 zha 1904. 55 5250: 58 9 
10a 1995. 87 50103. 5 a’ Fs 71733 — @ 1S —0918} 
3a a "Or mnt - A whet ‘ , . - } 
on ey «89 90130. r a *G,—76363 2a 1902. 916 52550. 9 a 5G ~139 
“ 1994. 24 50144. 4 a °G;—69283 5a 1901. 40 52592. 8 a Ca—7179s 
10a 1899. 14 52655. 4 a CG. aaa 
se 1993. 480 | 50163. 5 a 5G3—69303 or 1898. 36 52677. 0 a7S;—25P;_— | 
= 1993. 189 50170. 9 a °F ,—74063 oa 1897. 825 52691. 9 a °F .—7592 
l2cwa 1992. 53 50187. 5 a ®*D,—65013 = 
190 1988. 19 | 50297. 0 a °Gs—69143 20 1893. 60 52809. 5 a 5G.— 7194 
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10a 1981. 89 50456. 9 a 5G,—69303 _ 1885. 845 53026. 7 a 5P,—76363 
ba 1981. 217 50474. 0 a 5P;—75233 . oe : 
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= 1974. 65 50641. 9 a °D;—65573 la 1875 118 os a °G1—72133 
150. 1973. 986 50658. 9 a *D»—6501% — S5SoR. © a °F 4—7722 
50 1973. 13 50680. 9 iS, 7 
9 ~) gt VUVSY. * a *S;—z ?P 100ha Q7 9 raeqoy « — ‘ 
vs 1942. 446 50690. 0 a 5D, B557% “4 oee ren - a 2 a *D)3—6835 rl 
5a 1970. 84 50739. 8 a3F La 3 ; S69. 142 O34 19. 6 a‘*D 6727 4 
4 74633 ve 1865. 657 53600. 4 a Cys > 4063 9 
<i a 1862. 219 53699. 4 ( 6862 + 
— 1970. 68 50743. 9 a 5D),—656? 10a 1861. 94 53707. 4 - aa \ 
10a 1966. 55 50850. 5 a5D BERS 'f a ' 7284 , 
4 062. 58 “1022 2 . a ‘ - 7 
». 1 + 62. = D099. 3 a °C, 7193 200 LS58. 65 53802. 5 ais > 5p t] 
~0 1954 2u 51169. 5 a*D 66108 5a L856. 49 53865. 1 ~ aT | 
la 1952. 830 51207. 7 267¢ 97 25 1854. 00 53937. 4 ; D 5886 : 
aOi De 143943 - v PIN 4 a BRB | 
10 1853. 70 53046. 2 a 51), 68873 ; | 
l2a 1952. 45 51217. 7 a °D,—66103 10 1852. 35 53985. 5 a ), 6886 g 
“ : ; : a 3] 2843 , : | 
= sees. +. 51233. 7 a ®Dyp—65013 a 1851. 31 o4001. 2 1°G 284 V 
Bs 1948. 65 51317. 6 a *D,—6624 20 1849. O4 54082. 1 a 3h 9? 
m4 1948. 48 51322. 1 a 5D) 6567 ba 1846, 060 54169. 4 a Cy 43 
12a 1944. 02 51439. 8 a *D;—66373 = 1845. 916 | = 54173. 6 a *S:—7139 7 
L5a 1844. 61 54212. 0 a 51).—6914 ’ 
loa 1942. 19 51488. 3 a D, 663) e e 
a" 1940. 769 51526. 0 a 3F.—7407% ioe 1842. 98 94259. 9 a *Dy—6914 | 
<9 1940. O4 51545. 3 a *Ds—66473 I2a 1840. 05 54346. 3 a °1),—69173 
30a 1938. 224 51593. 6 a‘*D, 66473 18 1839 81 54353. 4 aD 69283 . 
la 1937. 857 51603. 4 a 5P 636 os 4 _ os401. | a ®D4—69285 I 
ee 3a 837. 73 54415. 0 a 5G 355 ; 
20 1936. 79 51631. 8 pS 
936. 7{ 516: a °D,—66513 27 Fe — 
10a 1935. 94 51654. 5 a 5D, $665. 1002 toss. SS 54420. 0 a °D4—6930 ! 
20h 1935. 32 51671. 0 a 081% ~n a 319 o 4516. 1 a®D 6886 
25 1932. 43 51748. 3 a °D,—66103 : 29, 992 pega a °Py—?77973 
1a 1931.372 | 51776.7 a 8P—75113 20. 1829. 992 | 54645. 05 a3F—7719 | | 
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TABLE 1. 
| Intensity | W avelength W a Term comb 
| spark vac number nation 
l 2 3 4 
A = — el 
J B 
5a 1824. 85 54799. 0 a °Ds—69723 
92 92 4 * 9 f a 5D, ~697 23 
12a 1823. 28 54846. 2 1 a 5D 69773 
3da 1821. 681 54894. 3 a ®Dy—69773 
18 1821. 03 54914. 0 a *D4—69793 
50 1816. 19 55060. 3 a *S;—z 7D3 
la 1812. 149 55183. 1 a ®Gs—74323 
10a 1810. 80 55224. 2 a 5G:—74073? 
la 1807. 873 55313. 6 a ®D;—70133 
5a 1807. 37 55329. 0 a ®G;—74463 
5a 1802. 54 55477. 3 a 5G:—74323 
10a 1797. 82 55622. 9 a 8:—72843 
10a 1796. 184 55673. 6 a *Ga—7613; 
50a 1792. 53 55787. 1 a 5Gs—7 4633 
3a 1787. 994 55928. 6 a 5D4—70813 
5a 1777. 11 56271. 1 a 5G:—75113 
2a 1771. 202 56458. 8 a ®D,—70813 
10a 1770. 97 56466. 2 a ®*D;—71395 
3a 1767. 65 56572. 3 a §D,—71393 
5a 1762. 84 56726. 6 a 5Gs—77193? 
5a 1762. 56 56735. 7 a 5D»,—7108}? 
5a 1759. 13 56846. 3 a °8,.—74073 
3a 1753. 83 57018. 1 a §D;—71943 
2a 1753. 020 57044. 4 a‘*D 71398 
18 1752. 85 57049. 9 a 7§,;,—2z 7D3 
200 1750. 14 57138. 3 a *S;—z5°Ds 
4. Energy Levels of Re II 
This term analysis of the second spectrum of 


rhenium yielded more or less information concerning 
49 even and 83 odd energy levels as shown im table 
2. In column 3, CS stands for Catalan and Sales, 
M for Meggers, and T for Tech. Durimg the course 
of the analysis deliberate effort was made to group 
the levels in spectral terms according to the general 
rules governing g- values, and level intervals. In 
the case of Re 1 [1] the interpretation of levels was 
guided by that for Mn 1, and similar clues as to the 
interpretation of Re um were sought by comparison 
with the homologous spectra Mn 11 [9], and Te 1m [10]. 
Indeed the ground states of these three homologous 
spectra are identical, d°(a°S)sa‘S;=0.0, except 
that the principal quantum numbers of d and s 
electrons increase by one unit in successive periods. 
Likewise the companion term d°(a%S)sa°S, is 
comparable in the three spectra except that it 
increases from 9472.9 K in Mn to 12617.2 K in 
Te u, to 17223.5 K in Re u. The analogies end 
here, for the next even level, a 5—D,, in Mn mr starts at 
14325.6 K, in Te m at 3461.3 K, and in Re 1 at 
14882.6 K. In Mn irand Te 1 these levels certainly 
arise from d°® because they are completely inverted. 
If the a °D term of Re u arose from d* it would have 
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Intensity a a Wave- Term combi- 
spark | va | number nation 
1 | 2 3 4 
A B 
3da 1748. 129 57204. 0 a 5D;—72133 | 
3a 1747. 42 57227. 2 a °G3—76363 
5ha 1744. 853 57311. 4 a *Dy—7108; 
5a 1738. 74 57512. 9 a 5Gq—77973 
2a 1738. 458 57522. 2 a ®Gr—76363 sO 
5a 1735. 93 57606. 0 a *D;—72533 
100 1725. 20 57964. 3 a 5Dg—7 2843 
5a 1704. 40 58671. 7 a ®D4—73553 
la 1690. 903 59140. 0 a ®*D3—74073 
40 1686. 87 59281. 4 a 7§;—z "Dj 
60 1685. 49 59329. 9 a 7§;—z D3 
25 1683. 86 59387. 4 a 78;—z °P3 
10a 1680. 67 59500. 1 a 5D, — 74323 
30 1677. 98 59595. 5 a 7S§3—z 5F4 
50 1676. 01 59665. 5 a "S;—59663 
25 1666. 66 60000. 2 a 58s—77223 
l5a 1646. 05 60751. 5 a °S:—77973 
l5a 1632. 53 61254. 6 a §D4—7613; 
20a 1590. 87 62858. 7 a *8;—z °Dj 
3a 1584. 97 63092. 7 a §D4—77973 
2 1584. 10 63127. 3 a *8;—6312; 
30 1552. 65 64406. 0 a 78;—64403 
3 1539. 84 64941. 8 a "S3—64943 
5 1525. 03 65572. 5 a *8;—65573 


very large negative intervals but since its intervals 
are relatively small and predominantly positive it 
must be ascribed to d‘s*. Theoretically, of course, 

d° and d*s’ produce the same number and types of 
spectral terms, viz, *D, *(HGFDP)'(IGFDS)*(FP) 

'(GDS), composed of 34 energy levels. A much 
larger family of even terms consisting of 74 levels 
theoretic ally arises from the d‘s configuration which 
produces 7S, °S, (GD) *(GD) *(IGFDS)' (GFDS) 
FP) (FP) *(HGFDP) '(HGFDP) *D, 'D. Thus 
the total number of levels expected from these two 
configurations of optical electrons is 108, whereas 
only 49 have been found in Re u. Moreover, the 
intervals between levels and terms is so large that 
strong configuration imteraction may account for 
large deviations from Landé g-values, and for many 
interval and intensity anomalies. For these reasons 
in Atomic Energy Levels [10, p. 168] electron con- 
figurations and spectral term designations were 
cautiously given for only 6 low-even terms of Re 11. 
Rather than list the remainder as_ miscellaneous 
levels in table 2, the advice of R. E. Trees was taken 
to attempt to group the levels and assign term 
symbols in a probable but provisional manner. 
However, since some terms are incomplete and many 
higher ones are missing, the provisionally designated 
terms have been omitted in the term combinations 


459 





5d 
Sd 


Config 
S) fis 
tis 

S)fis 
C7) he 


—_ . 9 fu } } 
Pasi 4 Ene gy levels and spectral lerms of Re uy 


I ler J La Interval Ot ( D> f 
M 4 if 2 { fit ity 
M 1377 0 . ms 
Ms | 1 14824 sng 1. 50 7 
M 50 2 ‘ <se I 1 
9 rin ve 2535p l 
M x ; 14930. 5 ! 
] 34 is 
a itt MI 
‘1 { 14882. ¢ l 
M 2 17223. 5 l 1 6 
M SO » ISS45.8 ) 
9 ’ 0. 64 
eis NM 
M 50 19139. 7 1.12 ‘ 
1323. 5 di ¢ I 
ss ae ; eine a ) ID M 
CS 4 M976. 4 1. OF ‘ 
254 N 
C'S 54 t 22 l 1. 24 a 
CS 54 ! 23722. 4 1. 90 ( 
M OO 2 23340. 8 1. 67 
1490 ; : Ld ; 
1429 dt Gs L))é ? 
M 7 : 24763. 1 1. 41 : , oa 
M 4 22544. 7 1. 2) : 
a neal 15. ¢ M 
24 
ad se | ; : (C'S 
2 | 
C'S ! 5 ss] ] ) 
‘ ; ’ 1 
M 2 0 oO 
cs 
papas " ) . . Cs 
sie 2 M 
CS M4 is 4 Cs 
- ( > 
Cs l ttt x 
1 4 
M 
M 50 M 
1 
- ; I 
MI MM t 1.49 ( 
M MM a l CS 
(‘Ss 1 | ‘ 
cs 
rT, 
M , 
cs 
140 ¢ 
M & 
M ) 
( 
( } . tos vp 
Ve" = } 
( 
M 0294.8 “i ‘s 
( 
Cs { ) ] 
M ' " aia ‘s 
( 
CS M4 1728 
j ( 
( 
CS 54 . CS 
RS ( 
c's | ] 
( 
( s 5 ’ 
{ 
cs 
CS 4 ; ‘ 
( 
( 4 } % 
CS M4 x 
( 
CS i ! ; 
‘ ( 
( f ies 
; A { 
( 
Cs t 
, | 
( 
we : ( 
af ' 
Nl 4 : - 
I 
CS 54 $ 


460 


tery 
25. 1 
1904. 1 
tit 
2191 
29.1 
1 ¢ 
1991 ¢ 


of 
foul 
A 
Re 
W 
hon 
it a 
be 
dp 
elec 
inv‘ 
the 
to t 
and 
beel 
sho’ 
the 
so V 
S 
ava 
ion 
pols 
cori 


1. 18 





of table 1 where they are represented by the first 
four digits of the energy levels. 

According to the above, the low even terms of 
Re 1 resemble those of the isoelectronic spectrum 
W 1 [10, p. 156] more closely than those of the 
homologous spectrum Mnu. For the same reason. 
it appears that the lower odd terms of Re 1 should 
he attributed to d* sp, as in W i, rather than to 
@p, asin Mn This arrangement avoids double 
electron jumps to produce the strong Re ut lines 
involving d‘s?a °D. It also explains why d°(*P)p 
y PF; could not be found [10] despite diligent search ; 
the two other levels, 99059.2 and 55150.8, assigned 
to this term have been reassigned to d‘s(*D)p z 7D, 
and the doubtful designations, z °G° and z *P° have 
heen abandoned. Beyond the first five odd terms 
shown In table 2, no attempt to group and designate 
the remaining levels has been made, because to do 
so would be highly arbitrary. 

Since no spectral series could be found among the 
available lines of Re ul, we tentatively adopt the 
ionization potential of 16.6+0.5 electron v extra- 
polated by Finkelnburg and Humbach [11]. This 
corresponds to a series limit at about 133900 K. 
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Density Formula for Alkali Silicate Glasses From 


Annealing to Glass-Processing Temperatures 
Leroy W. Tilton 


Vol. 61, No. 6, December 1958 Research Paper 2915 


An additive formula is given for the computation of specific volumes of molten and 
compacted alkali silicate binary glasses with modifiers up to 50-mole percent and from anneal- 


ing temperatures (400° or 500° C) to 1,400° C. 


silica are postulated as 


The effective partial volumes, vg, for the 


(1/vgs) = 2.198 +7,,Cs(1723 —2), 


where r, is the mole fraction of nonsilica and Cs is a constant to be evaluated from glass- 


density data. 


This is based on the idea that silica networks can contract in volume as 


temperatures are lowered provided, and in proportion as, modifier ions are present in the 


glass. 


The effective partial volumes of R,O, the nonsilica, are assumed as linear functions 


of temperature but not of the fraction of silica present. 


Computed and observed densities agree within approximately 1 percent. 


The effective 


volumes of nonsilica are somewhat smaller than published estimates for the oxides themselves, 
as should be expected because of interpenetrations. 


Previous formulas for the computation of densities 
of glasses have related principally to variations in 
composition. The formula proposed by Stevels [1] ' 
is based on ideas of structure and has very few con- 
stants. It assumes for restricted ranges in compo- 
sition that volume is independent of the nature of 
the cation. As modified by Stanworth [2] it shows 
that volume is dependent on the size of the cations, 
the glass volumes being larger as cation radii are 
larger. 

The additive formula proposed by Huggins [3] has 
many constants, and it is notably successful. That 
work is of considerable interest to anyone interested 
in structure because Huggins finds that the addition 
of every cation oxide molecule linearly increases the 
volume of a glass per oxygen atom, regardless of 
previous additions. However, at least for R,O 
silicate glasses, it can be shown that the linearity is 
not dependent on the referral to oxygen, and his 
corresponding partial volumes of nonsilica (not 
referred to oxygen) also are found to be linear if 
plotted against mole fraction of the added oxide.’ 

The Huggins formula is of further interest because 
the partial volumes used for silica are almost identical 
in the presence of various kinds of nonsilicas and are 
s0 treated except for the presence of a few elements, 
Most interesting of all in connection 
with the present paper, his partial molar volumes 
(P.M. V.) attributable to silica * decrease as larger 
amounts of modifier are added. The Huggins 
stepwise procedure in dealing with the silica compo- 
nent is roughly equivalent to the procedure men- 
tioned by LeChatelier [13], namely an arbitrary dual 





' Figures in brackets indicate the literature references at the end of this paper. 

? Although the Huggins partial volumes for silica (whether or not referred to 
orygen) are more nearly representable by a smooth curve when plotted against 
mole fraction of silica than against the Si to O ratio (used by Huggins), their 
better representation by segments of straight lines can still be shown. The 
Huggins special points of intersection correspond to compositions where there 
are integral numbers of cations per cage according to vitron theory [5, 6]. 
be able from Huggins’ table 1 as P. M. V.=2bgi+csi. See White [8], 
». 24 to 28. 





calculation of the volumes attributable to silica, 


partly from the density of silica glass and partly from 
that of the crystalline silica phases. Biltz and 


Weibke [14] proceeded in a slightly less arbitrary 


manner and calculated volumes of “average mols’”’ 


of glass by the additive rule after taking into account 


the specific volumes of certain crystal phases that 
could be derived from the equilibrium diagrams and 
therefore might be presumed to exist in the liquid 
silicates. 

; , : , 

Ihe procedures mentioned are certainly not incon- 
sistent with a conclusion from data on annealing [4] 
that the nonsilica particles, chiefly in proportion to 
their mere number, seem to govern the extent of 
partial collapse of the silica framework. A decrease 
in the bond angle Si-O-Si, where the tetrahedra are 
connected, was suggested. 

The object of this paper is toseeif the additivity of 
volumes in alkali silicate glasses at high temperatures 
is compatible with the idea of a sili¢a that is con- 
tinuously variable in density as a function of both 
temperature and modifier content. The following 
considerations also are relevant to the development 
of the proposed formula: 

(1) A bond-angle theory of annealing [4]; 

(2) Known linear variations of densities of glasses at 
equilibrium in the annealing’ ranges [4]; 

(3) Shermer’s [7] curves for densities of molten glasses 
as functions of composition, showing (qualitatively ‘) the 
marked increases in expansivities as cation is added to 
silica; 

(4) White’s [8] finding that densities of molten alkali 
silicates vary nearly linearly with temperature over a 400° 
C range; , 

(5) White’s report that the silicon-oxygen structure in 
molten silicates is not only similar to that in vitreous 
silica but is independent of the kind of cations present; and 

(6) White’s report that the partial molar volumes of 
silica are identical for the alkali silicates and decrease as 
cation is added. 


* The slopes of Shermer’s density curves from room temperature te 400° C are 


about 44 as steep as is indicated by the coefficients in his table 1. 
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The basis for development is the well-known 


additivity of volumes, so that the specific volume of 


a glass may be written: 
(1) 


where vs and rg are effective partial specific volumes 
attributable to the silica network and the nonsilica 
modifier, respectively, and Sy and rr, are the respec- 
tive weight fractions. 

For each modifier oxide, the effective density ts 
assumed as a linear function of temperature 


(1/vp) = Ag+ Be(1 400-1 


where 1,400° C is used because it will be desired 
later to compare directly the empirically determined 
A constants with White’s estimates at 1,400° C. 

The problem, then, is to find a suitable expression 
for vs, the effective volume of the silica. As a first 
approximation one might constant 
independent of temperature and modifier concentra- 
tion. If this be done then it is found necessary to 
express each of three values of rg not only as a 
function of temperature but of modifier content. It 
is simpler to assume that only silica volume, ¢s, is 
a function of modifier content. 

It is desired that 1/rg have the value 2.198 g/em 
at 1,723° C, the melting point, and at all tempera- 
tures in the absence of modifier oxides; also that it 
decrease monotonically with temperature for a given 
degree of modification, and increase with modifica- 
tion at constant temperatures other than 1,723° C 


use Ys as a 


Hence it is postulated that in the presence of any 
alkali modifier 
(1/¢s) =2.198+7r,.Cs(1,723—0), (2) 


where 7, is the mole fraction of nonsilica and Cg is 
a constant to be evaluated from the glass-density 
data. 


Other procedures ean, of course, be followed in 


seeking to express as a continuous function the 
changes with temperature and modification that 


seem necessary for successfully computing glass 
densities. It seems preferable to follow precedent 
in attributing the changes mainly to the silica 
component [3, 13, 14], as mentioned above. 

The assumption of eq (2) leads to much lower and 
more reasonable temperature coefficients, B, for 
UR than would otherwise be found necessary in eq (3 
Incidentally, the form of eq (2) accords with ideas 
previously expressed [4] concerning the nature of 
the structural rearrangements and volume changes 
in silicate glasses in their annealing ranges, namely 
that the silicon-oxygen bond is temperature sensitive 
(undergoes more electron cloud distortion, and be- 
comes more angularly directed) in 
nonsilica particles (metal ions) 
operation of formula (1) 


proportion as 
are present. The 
together with the auxiliary 
eq (2) and (3), depends importantly on the use of 
the factor 7,, as included in eq (2) for evaluating the 


effective partial specific volume, rs of the silica. 


In eq (2) the constant Cs is the slope or the tem 
perature-contraction rate, Ad/At, for a etl 
of vitreous silica collapsed as fully as can be aceall 
plished by the presence of alkali ions. Ip Othe 
words, the slopes of a family of curves of the effectiy, 
partial densities of silica itself (but permeated yi 
alkali modifiers) versus temperature are expressed br | 

rm Cs and the partial derivative of such slope 
respect to 7, 18 Cs. 

Krom experiments [4] on volume changes of cop, 
plex glasses in their annealing ranges only, and thy 
for an average of numerous types of modifiers, , 
value of Cs 0.00057 might be expected.®  Aety. 
ally, for use in eq (2) for computing densities of 
alkali silicates in the annealing range and up 
1,400° C, a still lower value of +0.00040 is nog 
empu ically recommended (after several trials With 
larger slopes °). 

The constants A and B were successively deter. 
mined from the glass densities after assuming various 
values of Cs. A value Cy=0.00040 was found to giy, 
the best fitting curves when residuals were averaged 
for all three modifiers at four temperatures, nameh 
annealing, 900°, 1,150°, and 1,400° C 

The “observed’’ values for molten 
computed ? from those that White |S] tabulated a 
1,400° C with coefficients for other temperature 
from near 1,000° to 1,500° C. White shows that his 
results at 1,400° C are in fair agreement with those 
by Heidtkamp and Endell {9} and those by Shartsis, 
Spinner, and Capps {10}. The “observed” values for 
compacted glasses at room temperature and at their 
critical temperatures are on the densities 
published by Young, Glaze, Faick, and Finn {1 
and on thermal coefficients tabulated by Shermer 
(7! and by Schmid, Finn, and Young {12}. 

The best effective densities, 1/vg, for use in eq (I 
are 





S With 


glasses are 


based 





l/r, 1.700 + 0.000332 (1,400—1?), 
l/¢yg=2.117+-0.000416 (1,400—2), | 
1 /rg = 2.066 +0.000428 (1,400—2). 


This implies that the effective densities of the | 
oxides in the silicates at 1,400° C are 1.70, 2.12, and 
2.07 which are 10, 12, and 14 percent higher than the 
densities computable for the oxides themselves from 
volume estimates as given by White [8] for Li,0, 
Na,O, and K,O, respectively. These percentages, 
can be considered to some extent indicative of inter! 





From eq (1) of the reference cited (An/A/ 0.00010r , and from the Newtom 
Laplace relation between refractive index and density, namely 1)=dC, one 
finds by differentiation that (Ad/A 1/0.176. Henee (Ad/Af 0.00057 Tm for | 
i silica network stabilized by the presence of various metallic ions, and G 

1/107 » +0. OST This value of Cs extrapolates to a molar volume ¢ 
18.7 cm! for stabilized silica network at absolute zero, as compared with 27.3 @ 
for metastable silica glass This would be lica with 39-percent voids, even if ont 
may consider oxygen and silicon as solid spheres with radii, 1.31 und 08 BI 
respectively Such a condition is then not unre onable because it is intermediat 
between the known conditions for densest packing of sphere NH} voids) and f@ 
open packing (48 voids 

It can be shown that the higher value of Cs—0.00057 can be used if s threshold | 
value of r»--0.12 is assumed for the beginning of configurational thermal expal 
sion See White IS p 106 

The computations are in several cas (tray itions from Whit s dats 0 
temperatures below the liquidii, partic rly 1! 1 potash glasses * 
vou) ( 
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Table 1. 


Comparison of computed and observed densities 


Mole Annealing tempe 1 m0? ¢ 1,150° C 1,400° C 
per 
m0 Computed Obsery c— 1 Computed Observed c—o) X10 Computed Observed e—o) X10" Computed | Observed | (e—o) X108 
Li.O-S;O. 
10 2. 232 2. 239 7 2. 213 2. 183 2. 180 
0 | 2.263 | 2.275 12 | 2.215 | 2. 206 +9 | 2181 | 2 187 6 | 2.152 | 2. 169 17 
30 2. 283 2. 298 15 2. 208 2. 216 8 2. 164 2. 177 13 2. 120 2. 138 18 
40 2. 209 2. 313 14 2. 197 2. 194 3 2. 138 2. 146 s 2. O83 2. O98 15 
50) 2 208 2. 172 2. 167 5 2. 104 2. 11) 7 2. 040 2. 055 15 
60 2. 137 2. 108 29 2. 062 2. 044 18 1. 988 1. 980 +8 
NagO-S\O2 
10 2. 272 2. 276 1 2. 242 2. 223 2. 202 2. 210 S 
20 2. 342 2. 35 12 2. 278 2. 271 7 2. 241 2. 247 6 2. 202 2. 222 20 
30 2. 404 2. 420 16 2. 306 2. 308 2 2. 252 2. 259 7 2. 200 2. 210 10 
10 2. 457 2. 468 11 2. 331 2. 322 q 2. 263 2. 261 +2 2. 194 2. 201 7 
50 2. 496 2. 514 18 2. 344 2. 336 8 2. 267 2. 264 +3 2. 187 2. 192 5 
60 2. 524 2. 352 2. 323 29 2. 264 2. 243 21 2. 176 2. 163 +13 
K.O-S;O 
10 2. 283 2. 282 | 2. 242 2. 216 + 26 2. 217 2. 216 + J 2. 191 2. 216 25 
20 2. 352 2.355 3 2. 271 2. 262 7) 2. 225 2. 223 +2 2. 179 2. 184 5 
30 2. 413 2. 405 S 2. 293 2. 284 + 2. 229 2. 220 + 2. 166 2. 156 10 
10) 2. 457 2. 307 2. 290 17 2. 230 2. 213 17 2. 151 2. 136 5 
50 2. 491 2. 313 2. 226 2. 140 


penetration of the silica network cavities with these 
modifiers. Correspondingly, in eq (3) the tempera- 
ture coefficients of density are 6, 8, and 16 percent 
higher than those computable from differencing 
White’s estimates of oxide volumes at room tempera- 
ture and at 1,400° C 

The complete formula yields densities of which 
two-thirds agree with observed densities within 
one-half percent over the range 1,400° to 400° or 
500° C for both molten glasses and those in the 
annealing ranges, as shown in table 1. 

From the residuals, c-o, in table 1 it is evident 
that the proposed density formula fits the data 
fairly well over the long range in temperature. This 
is illustrated in figures 1, 2, and 3 where, for the three 
systems, rough comparisons are also made with data 
by observers other than White. The slopes in 
figures 1, 2, and 3 show how regularly the expansivi- 
ties Increase as modifiers are added. The approx- 
imate constancy expansivity from annealing 
temperatures to of molten glasses accords 
with the idea that glass structure does not change 
fundamentally within these temperature ranges. 

When more numerous and more precise data 
become available it may be found, for glasses nearly 
saturated with modifiers, that observed densities 
in the annealing ranges differ systematically from 
those computed by this formula as adjusted to fit 
densities not only for annealing ranges but for 
molten glasses also This expectation is based on 
the postulated characteristics of vitrons. Saturated 
5, 6] cannot be expected to respond 


of 


those 


cages of vitrons |: 
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as freely to bond-angle contraction during cooling 
as is assumed for unsaturated cages considered in 
this simple formula. The resulting stresses that 
develop on further cooling may eventually cause 
some breaking of bonds so that the delayed bond- 
angle contraction can then proceed. Some indication 
of such systematic tendencies is shown by the 
predominantly positive residuals in the 900° column 
of table 1, and the almost exclusively negative 
residuals in the annealing-temperature column. 

As modifiers are added, there is in these residuals 
another and more definite systematic trend. The 
computed densities must be regarded as averages 
over the composition ranges because the formula 
has no terms relating to degree of saturation. The 
computed densities are too low at low modifier 
content and become too high at higher percentages 
of nonsilica particularly for molten glasses at highest 
temperatures listed. This is clearly indicated in 
figure 4. 

It will be noticed that maximum excess of ob- 
served densities tends to shift toward higher per- 
centage modifications in the order K,O, Na,O, Li,0. 
Indeed the maxima for K.O and Na.O occur at 
lower percentages than 16.7 and 23.1, respectively, 
at which it has been suggested that complete in- 
closure of oxygen in the dodecahedral cages should 
cease and the breaking of oxygen bridges should 
begin. It happens also, that the excesses change to 
deficiencies not far from the percentages 28.6, 37.5, 
and 50 that have been identified with complete 
saturation of the cages. 
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In addition to the results listed in table 1, data 
for two ternary glasses investigated by White were 


used. Each contained 60-mole percent of SiO, and 
20 percent of Li,O. For one with 20 percent of 
Na,O the computed and observed densities at 


1,400° C were 2.147 and 2.155; similarly for one 
with 20 percent of K,O, 2.123, and 2.126, respec- 
tively. 


In resumé the following points are emphasized: 


(1) A simple additive formula [1], if written in accord 
with an effective partial value of silica as postulated by 
eq (2), represents the data on alkali silicates approximately 
as well as they are known. The systematic trends in the 
residuals seem qualitatively in accord with expectations 
from vitron theory. 

(2) The use of a factor rm in eq (2) is analogous to 
Huggins’ [3] use for silica of a parameter variable with 
modifier concentration. The need for such a procedure 
was indirectly implied by White [8] and by Callow [15]. 
(See considerations as above listed.) 

(3) The empirically determined constant, Cg, for silica 
in eq (2) is reasonable in that it is interpretable as a 
temperature collapse rate for silica network that extrap- 
olates to reasonable volumes at 0° K. 

(4) The empirically determined constants for nonsilica 
in eq (3) are reasonable in that they approximate the 
densities of the oxides and their temperature coefficients 
within 10 or 15 percent of estimated values. 


It would appear from the analysis in this paper 
that to a first approximation the sensitive-bond-angle 
concept for silica with alkali modifiers is semi- 
quantitatively compatible with our knowledge of the 
densities of alklali silicate glasses. 


References 

[1] J. M. Stevels, J. Soc. Glass Technol. 30, 31 (1946). 

[2] J. EE. Stanworth, J. Soc. Glass Technol. 30, 54 (1946). 

[3] M. L. Huggins, J. Opt. Soc. Am. 30, 420 (1940). 

[4] L. W. Tilton, J. Soc. Glass Technol. 40, 338T (1956). 

[5] L. W. Tilton, J. Research NBS 89, 139 (1957) RP2782. 

[6] L. W. Tilton, J. Research NBS 60, 351 (1958) RP2854. 

(7| H. F. Shermer, J. Research NBS 57, 97 (1956) RP2698. 

[8] J. L. White (Ph. D. Thesis, Univ. of California 1955). 

[9] G. Heidtkamp and K. Endell, Glastech. Ber. 14, 89 
(1936). 

[10] L. Shartsis, 8S. Spinner, and W. Capps, J. Am. Ceram. 
Soc. 35, 155 (1952). 

fli] J. C. Young, F. W. Glaze, C. A. Faick, and A. N. Finn, 
J. Research NBS 22, 453 (1939) RP1197. 

[12] B. C. Sehmid, A. N. Finn, J. C. Young, J. Research 
NBS 12, 421 (1934) RP667. 

[13] H. LeChatelier, Kieselsaiire und Silikate, p. 240 (Leipzig, 
1920), 

[14] W. Biltz and F. Weibke, Z. anorg. allgem. Chem. 203, 
345 (1932). 

[15] R. W. Callow, J. Soc. Glass Technol. 36, 137 (1952). 


467 








E 


Thermodynamic Preperties of Gases at High 
Temperature: 


fournal of Research of the National Bureau of Standards Vol. 61, No. 6, December 1958 Research Paper 2916 


| J. Chemical Equilibrium Among Molecules, Atoms, 


and Atomic lons Considered as Clusters 
Harold W. Woolley 


The equilibrium thermodynamic properties of gaseous mixtures at high temperatures 
are treated by an extension of the cluster theory of Ursell, omitting the assumption of addi- 
tivity of pair energies. An effective partition function is introduced, which is a convenient 
function in expressing the increase of the partition function for the entire gas due to the 
joining together of the parts of the cluster. The law of mass action and the dependence 
of second and third virial coefficients upon the cluster integrals for pairs and triples of 
molecules in a mixture have been obtained therefrom. Extension to a partially ionized 
gas is made by incorporating Mayer’s ciuster-based extension of the Debye-Hiickel theory. 


1. Introduction 





The calculation of the thermodynamic properties of gases at high temperatures involves 
a number of special problems which may well be disregarded at more moderate temperatures. 
The thermodynamic functions for a gaseous constituent in a mixture in chemical equilibrium 
are usually calculated from a partition function definable as a sum of Boltzmann factors for 
all states of the elementary unit of the constituent. The states of a single isolated molecule, 
atom, or atomic ion may be arranged in some kind of joint or compound progression of quantum 
numbers up through a range of moderate energies. For still greater energies, however, the 
number of states associated with the physical system comprising the constituent is enormously 
greater, with the states forming a continuum of values nominally of infinite extent, as the 
constituent dissociates, rather than occurring discretely. At some elevated temperature and 
beyond, the usual criterion for neglect of high-energy states based on the smallness of the 
Boltzmann factor is not fulfilled below the dissociation limit. The infinity of states immedi- 
ately encountered at this energy does not permit the continuation of the same type of sum of 


states to indefinitely higher energies. A common procedure for avoiding the divergence which 


threatens to occur involves the termination of the sum of states at the dissociation energy or 
the inclusion of only that limited group of states above dissociation which still imply some of 


the typical structure of the undissociated constituent. 


A different procedure for avoiding this divergence difficulty is to be discussed, involving 
the use of a method thought to have the merit of formal correctness, although it is admittedly 
not the only type of analysis with this appropriate feature. For this purpose, a special kind 


of partition function will be used, here to be termed an “effective partition function”’. 


effective partition functions are related to the increase in the occupancy of states of a combined 
unit due to the joining together of parts of the unit, putting it together from any of all of the 
ways in which the unit may be considered as separated into component parts. The present 
| procedure comprises an analysis using clusters of the Ursell type [1]! to specify the net 


thermodynamic effects of particle groupings in the combinations involved in forming mole- 
| cules, atoms, and ions. As the populations of these clusters are measures of contribution to 


concentration due to interactions among several particles, it is admitted that they can give 


contributions of reversed sign if exclusion effects exceed the effects of combination. 


The calculation procedure presented is a treatment of a multicomponent system including 
the principal quantum effects remaining at ordinary or elevated temperatures. As such it is 


in some degree a generalization of the treatments for pure substances given by Kahnand 
| Uhlenbeck [2], Kahn [3], and Kilpatrick [4], and for the nonquantum multicomponent treat- 
: ments of McMillan and Mayer [5], Fuchs [6], and Muto [7]. A restricting assumption of 
| additivity of pair energies, made in some of these treatments, has not been retained in the 





Figures in brackets indicate the literature references at the end of this paper. 
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present one. The assumption is considered to be a rather good approximation for mixtures of 
ordinary molecules but not universally for the many types of interaction occurring among 
the chemically unsaturated fragments of molecules existing at very high temperatures. 

Though the cluster integral procedure provides a way of avoiding divergence difficulties 
associated with ordinary dissociation, there is a type of dissociation for which the separating 
parts of the molecule or cluster are oppositely charged. The distribution of rotational states 
for which such fragments occur as unexcited parts is similar to that existing when a single 
electron is nearly separated from a positive ion. In the latter case, just below the energy at 
which the continuum due to ionization begins, there is an infinity of discrete states, and each 
level of the outer electron specified by the principal quantum number n is actually a group of 
states whose number has 2n? as a factor. In the method of Urey [8] and Fermi [9], discussed 
also by Fowler [10], the summation was halted at such electron orbits that at the given density 
an orbit based on one atomic ion would just “touch” an orbit of a neighboring atomic ion, 
Their method, which appears capable of extensive development and application, forms part of 
the basis for a treatment of high temperature thermodynamic properties of a gas given by 
Megreblian [11]. The method of Planck [12] involved the evaluation of a Gibbs phase integral 
with simplifying assumptions. The present treatment is somewhat related to that of Planck 
but takes some account of a known property of an ionic assembly in involving an exponential 
Debye shielding function multiplying the long range Coulombic ionic interaction, such as oc- 
curs in the Debye-Hiickel theory of electrolytes. In their theory of a plasma, Bohm and Pines 
[13] made a canonical transformation to new coordinates of two kinds, one associated with 
pure translational motions, and one associated with a short-range interaction showing at least 
some resemblance to the Debye screened Coulomb potential. Of more immediate application, 
Maver [14] derived thermodynamic quantities for ionic solutions including all ring-type graphs 
of interionic interaction, using an exponential shielding function multiplying the Coulombic 
ionic interaction as a tool in the mathematical analysis. In effect, as shown by Meeron [15], 
Mayer converted the calculated effect of interionic interactions as indicated by virial coefficients 
from a nonconvergent form for the unscreened Coulomb potential to a convergent form with the 
Debye screened potential and obtained also the well-known limiting law Debye-Hiickel term. 
The avoiding of divergence in the integrals related to virial coefficients by using Debye screening 
corresponds in some degree with the avoiding of divergence by Ecker and Weizel [16] who noted 
that n, the principal quantum number for atomic states, must have an upper bound if the 
electron moves in a Debye screened potential field. The correspondence is not that of identity, 
however, since the virial-coefficient-type quantities receive part of their contribution from 
regions of phase space above the dissociation, or ionization energy. 

As to the arrangement of the present paper, section 2, which follows, presents a number of 
items of cluster theory leading to the method of calculation used to treat chemical equilibrium 
and the virial coefficients related to further clustering of the aggregations considered as the 
chemical constituents. Many of the relations are so similar to those in earlier derivations of 
cluster theory as clearly to form obvious extensions from them. A somewhat extensive dis- 
cussion pertaining to an analysis in rather elementary algebraic form has been included, how- 
ever, which may serve to show the very simple physical meaning of the extended Ursell-type 
formulas and thereby permit an easy inference of the essential correctness of the treatment. 
The effect of quantization on the cluster integrals or effective partition functions for molecules 
is also treated. 

In section 3, formulas are given for the Helmholtz free energy for the gaseous equilibrium 
mixture including molecules, atoms, and atomic ions, based in part on the ionic-solution theory 
of Mayer and in part on the equilibrium theory of clusters of section 2. Formulas for the 
equilibrium constants for chemical reactions which can occur among constituents are also given. 
Expressions for other themodynamic functions for the mixture in terms of functions for its 
constituents are given in the appendix. 

Further discussion of the status of cluster theory, particularly in regard to its present 
application, will be found in section 4. 
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2. Application of Cluster Theory 


Apart from effects associated with ionization, an accurate theory for the thermodynamic 
properties of equilibrium mixtures of chemically reacting substances at high temperatures may 
be directly formulated involving equilibrium constants identical with the quantum equivalent 
of the cluster integrals of gas theory. In the nonquantum treatment, the Gibbs phase integral 


; pier’ N ‘ “ . 
@ IN. J — - fexp pp h “SB pdr, (2.1) 


is analyzed according to the method of Ursell [1]. The extension to the exact quantum mechan- 
ical treatment involves the replacement of the Boltzmann factor by a Slater sum so as to 


provide an exact equivalent to 


Q=di exp (2.2) 


The treatment should give results similar to those obtained by Kahn and Uhlenbeck [2], and 


Kahn [3], for a pure gaseous substance. 
The theoretical analysis is made in terms of cluster integrals 6,,, defined by 


° ad t & 
Vb Hl! | fumes, (2.3) 


with U’, to be defined in terms of products of quantities W, to be described. The number of 
particles of type 7 within a given cluster is /;, with the different types of particles extending to 
i=t. This general cluster will be said to have composition (/,). The cluster is to be thought 
of as subdivided arbitrarily into fragments, including every possible way of subdividing the 
cluster. The number of particles of type 7 within a particular type j of fragment of the cluster 
is to be indicated as (/,); and the number of fragments of type j by k;. Conservation of matter 
then indicates that 

ald j=l (2.4) 


4 


for each type of particle in the cluster. The total number of fragments according to a particular 


way of subdividing the cluster is 
k= 2 (2.5) 


The quantity (7, may now be defined as the sum 


, l Wig 
‘(a \k-1 (J 
=, l k—1 Ne I 


II / , 2.6) 
with the factor g, arbitrarily separated to indicate an explicit counting of spin multiplicities 
for the fragments, if desired. The summation in eq (2.6) includes products formed according 
to every possible way of subdividing the cluster into fragments or groups of particles. Each 
product term involves different W,’s which together involve once and only once each particie 
of the cluster. W, represents the probability of occurrence of an isolated fragment j of the 
cluster in a given geometric configuration as a relative probability compared with the occurrence 
of the fragment when there is no internal energy of interaction. Thus, W, will represent a 
Boltzmann factor in a nonquantum treatment or a suitable Slater sum in a quantum treatment, 
such as is indicated by 


W =I (h?/2am,kT)* > y .Weexp[—e,/kT]. (2.7) 
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The integral of the Slater sum for each fragment or group of particles provides the sum over its 
quantized states, as far as quantized states exist, with the appropriate integral in the continuum 
for the remaining contribution. A number, q;, of states of like spin may have essentially 
identical energy and then can naturally be counted together. The separation of this factor 
from eq (2.7) as written in eq (2.6) is clearly arbitrary. When completely separated, the spin 


factor IIq,*/ is the same as II (2s,+-1)‘, based on the individual spins of the constituent particles, 


0 


j 
The denominators II(/,);! under the W,’s in eq (2.6) occur as shown if the eigenfunctions in 
eq (2.7) are regarded as not restricted by the exclusion principle but instead as having the full 
number of states corresponding to the Boltzmann statistics as a limiting condition. If the 
Slater sum were obtained exactly according to the exclusion principle, then a compensatory 
multiplication by this factor would be required, as indicated in the theory for a pure gaseous 
substance as given by Kahn and Uhlenbeck and by Kahn. 
With the preceding formulation, the cluster integrals are the equilibrium constants for 


cluster formation according to the equation 


N N 
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where V is the volume accessible to the gas, where ,,;,, is the number of clusters of composition 
(/,), and where N,, is the number of separately occurring individual particles of type 7, i. e., 
as ‘‘clusters’’ of one particle. This result can be inferred in accord with the analogy between 
the present case and the treatment for clusters for a single pure substance [2, 3]. Alternately, 
the result can similarly be easily inferred from the nonquantum formulation of Fuchs [6], on 
the grounds that the present formulation is equivalent except for its inclusion of quantum 
effects. Further indication of the correctness of these equations will be presented later in this 
section. 

It is to be noted that the equilibrium constants of eq (2.8) are the theoretically exact 
ones for a treatment of the real gas as a mixture of ideal gas components. Hill [17] has 
pointed out that any modification of definition that would give clusters other than these 
(he terms the present ones “mathematical clusters’’), would eliminate the ideal gas law de- 
pendence of each component cluster gas; this would not allow the simple mass action law of 
eq (2.8). Excluded volume effects are automatically included in our theoretical formulation 
at this stage and would not be added separately unless a compensating readjustment were 
introduced to remove an equivalent quantity from the cluster representation. It may be 
recognized that a cluster of a given composition can be called a molecule or free radical or ion, 
or a group of such entities, according to the normal custom for nomenclature. If a molecule 
happens to have various isomeric forms, the requirements of the particular application could 
determine whether the different forms would be treated together or separately according to a 
reasonable apportionment of accessible phase space. 

The Helmholtz free energy for a mole of like clusters of type (/,) in the absence of 


ionization may be represented as A=k71nQ with 


Q= [Qi 1Qr]*/N! 
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The first bracket indicates Q,‘, the effective internal partition function, and the second in- 
dicates Q,,, the translational partition function for a single cluster. 
From eq (2.9) with eq (2.3) and (2.6), one may infer that the effective partition function 


for the cluster may be represented formally as 
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with dV; taken as 
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in accord with eq (2.1) and (2.7), since each W, in eq (2.6) is a function of coordinates internal 
to the fragment /. 

Even though the equations given follow plausibly from earlier investigations, it may be 
advantageous to examine some of the relations in greater detail. By such a procedure we may 
hope (1) to confirm definitely that the Ursell type formula is actually being extended cor- 
rectly to mixtures with “non-additive”’ potentials, (2) to show that the structure of the cor- 
rect Ursell-type formulas is simply that which provides the incremental contribution (per 
cluster) to the state sum, due to the particles being in the cluster rather than present in other, 
less combined forms, and further, (3) to obtain the result that the formation of a combined 
cluster from two component clusters is strictly analogous to the combination of two particles 
into a cluster of two, so that repeated application of a process analogous to calculating the 
second virial coefficient gives an analysis concerning clusters of any degree of complexity. 

Thus it is informative to compare the effective partition functions for two cases differing 
in that in one case one additional particle is present within the cluster while in the other case 
it is counted as present in the assembly without interaction with other constituents of the 
cluster. The difference in magnitude between the state sum for the entire gas when the 
particle is combined and its value when the particle is uncombined includes the effect of in- 
serting the particle into every subcluster or fragment of the cluster involved according to eq 
(2.6) or (2.10) for the effective partition function for the initial cluster. The numbers of sub- 
clusters after insertion of a single particle into one of them include k,, of the type of the 
receptor subcluster (initially there were k,,+-1 of this type), &;, of the type of the product 
subcluster (initially there were k,;,—1 of this type). The number of single particle subclusters 
for particles of type g after addition of the added particle of this type is to be indicated as 
k,,. The increment in the total effective partition function for the system including the 
initial cluster and the added particle due to addition of the particle to the cluster will be indicated 
by 


Plaka .. >be « +) =OOGle...t—l,.sd 2.12) 


where @ is an operator operating on Q*(/,, /. . . . l,—1, . . .), the effective partition function 
for the initial cluster. The total contribution to the effective partition function of the com- 
bined cluster formed by adding the new particle includes a contribution from each factor 
of each term of the partition function for the initial cluster. From this it can be seen that the 
numerical factors arising in use of the operator @ are identical with those occurring in a 
differentiation process. @ can be represented as Q,*—Q,, in which Q,* is an operator which 
produces the ordinary partition function of a combined system when it operates on the 
ordinary partition function of the partial system, as Q:*Q,=@,@>, while Q; is an ordinary 





partition function of the added particle. 

Making use of eq (2.10) for the effective partition function of the initial cluster, the 
effective partition function of the combined cluster obtained with this operator @ is indicated 
in accord with eq (2.6) and (2.11) using the corresponding (7 / !in terms of the W’s, now 


to be written omitting the spin factor IT(2s;+-1): 
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(2.13) 
On the right-hand side of the equation, the expression before the bracket is the regular ex- 


pression for the combined cluster. The expression within the bracket includes contributions 
of various kinds arising from the addition of the particle to various possible fragments of the 
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initial cluster. For the first term in the bracket, coming from the insertion of the added par- 
ticle to a subcluster, the factor £,,+-1 in the numerator is due to there having been this number 
of subclusters to which the particle might be added to produce one particular type of product 
cluster. This modifies a particular product of the type preceding the bracket since there had 
been a factor &,,+-1 in the old factorial in the denominator and this factor now appears in the 
denominator of the term in the bracket. Similarly the new factor outside has a new factor 
k;, in a factorial in its denominator (included in one of the (4,)! factors) which is balanced to 
the original by the £,, in the numerator of the term in the bracket. Also, a new factor (/,), is 
in the denominator in front of the bracket (included as a net increase of one factor in the 
general (/,);!) because of the &,; indicating /,, instead of a former /;,—1, and this is balanced 
in the term in the bracket by the (/,),; in its numerator. 

The other two terms in the bracket arise from the addition of the new particle as a separate 
subcluster arising in the —Q, term of the @. For the initial cluster, the value of k was less by 
one than its value in such corresponding terms of the combined cluster, as are obtained using 
the Q, factor thus giving the opposite sign by (—1)*, which, when combined with the minus 
sign on the—Q), gives positive contributions within the bracket for both terms. The occurrence 
of k—1 in their denominators is a compensation concerning the (k—1)! factor and is also due tok 
for the terms for the initial cluster being one less than for the combined cluster. The first of 
these last two terms in the bracket has an origin somewhat analogous to that of its predecessor 
but with Q, involved rather than the combining Qf. The prime on the summation sign is to 
indicate that the summation is to include only those 7’s for which the k’s do not represent k,,, 
so as to obtain the count of the number of terms which had been involved for the Qf, except 
for those due to adding a single particle to a like particle. The effect of adding the separate 
single particle subcluster is to increase the exponent in front of the bracket for single particles 
q to k,, from k;,—1, its value in the initial cluster. To compensate for the matching extra 
factor in the factorial expression in the denominator preceding the bracket, the value k;, occurs 
in the numerator of the present term in the bracket. The last term in the bracket is likewise 
analogous to the first term, with Q; involved rather than the combining Qf. This term counts 
the number of separate occurrences in which the added particle enters into the formation of 
clusters of two like particles. This term appearing as a balancing quantity to which the effect 
of QF is compared, has the factor & 1 for the number of initial subclusters in connection with 
which this final comparison process arises. In addition, the numerator contains the factor k; ,, 
as in the term preceding, to balance the augmented factorial in the denominator of the ex- 
pression preceding the bracket. 

The quantities in the bracket may now be combined. For the first term, /,;, occurs only 
for the subclusters produced by adding the new particle to a previously existing subcluster. 
Thus, it does not include 78 for single particle subclusters, so that 


Sr (Ue) p= Da (le) yee = Le — Boje 2.14 


—— j ommed 7 


For the second term, >4’k,4;,, the summation is over j’s for which the /,’s do not represent k 

As the full sum S*%,; would have been equal tok, S\’k; is | k.. The three terms are now 
_-— —s 

readily combined as 
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which reduces to he This is cancelled by the l, in the denominator after the bracket, represent- 
ing an additional factor in the general factorial of the denominator needed to convert from 
phase integral to partition function as indicated in eq (2.1). From this demonstration eq (2.6), 
or its equivalent, eq (2.10), may be seen to be an appropriate extension of the ordinary Ursell- 
tvpe formula to mixtures. Thus, an effective partition function agreeing with the Ursell-type 
formula may be generated relevant to the combined cluster, starting with an Ursell-type 
formula for the cluster of one less particle by simply evaluating the increase in the net state 
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sum due to counting an additional particle of the entire system as belonging to the cluster. On 
the one hand, one may see this as a logically correct process to build up the effective partition 
function, so that on the basis of the obvious correctness of the Ursell-type formula for small 
clusters of one or two particles, mathematical induction permits one to regard the present 
Ursell-type formula as clearly holding for clusters of arbitrary size and composition, but with 
due regard for details of particle statistics as required. Alternatively, for one who finds the 
present Ursell-type formula to be an obviously correct extension of the Ursell formulas previ- 
ously known, the foregoing demonstration can show that the addition of a particle to a previous 
cluster can be treated formally in the simple form of a second virial coefficient as indicated, 
using the operator 6=QT—Q, with eq (2.12). It may also be noted that the joining of one 
cluster with another to form a combined cluster may be treated similarly using an operator @ 
built from effective partition functions 


0— Qe*— Qi, (2.16) 
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= O06 = 1105 — Gi (2.17) 
has a form identical with the pattern for obtaining the second virial coefficient for the effects 
of interaction between the two clusters. Thus, the calculation of the effective partition function 
for a diatomic molecule is directly related to and essentially identical with the calculation of the 
second virial coefficient for the interaction of its constituent atoms, no matter whether one 
regards the atoms as interacting particles or clusters of nuclei and electrons. Similarly, the 
treatment for a triatomic molecule is directly related to that for a cluster of three particles such 
as is involved in the cluster theory for the third virial coefficient. In all cases, states of higher 
electronic excitation should be allowed for if significantly present and should be understood as 
included in the present notation. 

We now come to the particular advantage obtained by use of the effective partition func- 
tions as defined in the foregoing discussion. The partition function for the complete gas as- 
sembly can be written out readily when they are used in the analysis. Using total effective 
partition functions Q5=@Q5(Q) trans) for a cluster for which the number of particles of kind 
iis (/;),, the partition funetion for an assembly involving .V; particles of each kind 7 is 


Qs S[ mia m4 ms], (2.18) 
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with the sum including every kind of term consistent with the condition 


> m,(/ 
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N,. (2.19) 


The most probable state may readily be found in the sense of determining the greatest term in 
the sum. This may be done on the basis that, in the immediate vicinity of the maximum, 
neighboring terms are almost exactly equal. If for a neighboring term, the number of clusters 
of tvpe 7 is taken as m,+-n,; (where each n, is a small integer) rather than the m, of eq (2.18), 
then eq (2.19) takes the form >3(m;+n,)(/);—N;, and the difference between the two is the 
general conservation equation 

Sin, (U;) ;=0. (2.20) 
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The requirement that neighboring terms near the maximum be approximately equal gives 
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since m;>1 and n; is a small integer. By substituting QQ; ans for @§ and dividing by 
IV", the relation 
l 


” n 


ui( ud ) j 11(Q5)" @ =) , (2.24) 


is obtained. Since Q contains V as a factor as shown in eq (2.9), or according to 


trans 


CV erene VA * ’ (2.25) 


with 
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the right-hand side of eq (2.24) is independent of V, and by further reduction using items 
relevant to eq (2.9), is equal to II(4,,),;" and forms an equilibrium constant to which the 


product of concentrations on the left-hand side of the equation is equated. The reaction 


equation to which this law of mass action is relevant is 
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where n,; represents both nj’ and —n‘{, with ); representing the chemical symbol for the type 
of cluster (or molecule) referred to. For the reaction of formation of a single cluster (or 
molecule) from elementary particles or clusters, the mass action equation can be reduced to 
eq (2.8), with V.,, corresponding to the m, of the equations just concluded 

If the representation of gas properties is carried out directly in terms of particular kinds 
of clusters, then these clusters may be considered afresh as the constituents of a gas mixture 
whose interactions lead to the virial coefficients for this mixture gas. These virial coefficients 
for the gas are of course related to the cluster integrals—or the effective partition functions 
for the clusters formed among its constituents. An adequate theory representing their effects 
formally can be given using the grand partition function. Exactly the same results may be 
obtained in the present case, however, by noting that the gas acts as a mixture of ideal gases 
subject to chemical equilibrium, as shown by eq (2.24). Then, the pressure of the gas is the 
sum of the partial pressures of all kinds of clusters of the constituents. Similarly, the net 
concentration of any kind of constituent particle is proportional to an analogous sum in which 
each term of the pressure sum is multiplied by the number of that kind of constituent particle 
in the corresponding cluster. While the various terms of these two expressions are expressible 
using cluster-integral-type or effective-partition-function-tvpe equilibrium constants, multiplied 
by products of powers of density-type fugacities of the various constituents, the ideal-gas- 
type relations which apply to each constituent permit the density-type relations which apply 
to each constituent permit the density-ty pe fugacity to be equated to the quotient of the partial 
pressure of the constituent, divided by &7, so that 
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The mole fraction of constituent ¢ is z;=N,/N. The process of eliminating the fugacities 
(pJkT) from the expression for pV/NkT is more involved than that for a pure substance 
[18]. Using 5,, and b,,, as the cluster integrals 6,,,. for two particles 7 and 7 and for three 


particles i, j, and k, respectively, the first two virial coefficients are obtained as 


B= D02:2;B,; (2.31) 
and 
C du eutitjstrCur (2.32) 
with 
B,, Nb,,. (2.33) 
and with 
N? 
Cua —— [405 bj 4-40 by 4- 4p ,b, 6b; 4). (2.34) 


Higher virials may similarly be obtained, but much less readily due to the proliferation of 
algebriac expressions in the transformations. Even for the third virial, the problem of non- 
additivity of pair potentials arises. For pure substances, Rosen [19] found nonadditive effects 
in the three helium atom interaction, Axilrod [20] has found such effects for the dispersion 
energies, and Jansen [21] has examined the entire problem at some length. Kihara [22] also 
has recognized the problem and has made use of Axilrod’s results in estimating the effect on 
third virials for the pure noble gases. For gas mixtures, the problem is of added complexity 
as it involves the estimate of interaction between two unlike constituents as well as the non- 
additivity among all the differing kinds of constituents considered three at a time, alike and 
unlike. In applications at elevated temperatures where the constituents may be other than 
ordinary stable substances, nonadditivity can be expected to have large effects. Where the 
prediction of properties of a gas mixture using rigorous theory is particularly impracticable, 
the use of empirical mixture rules becomes a natural expedient. 


If one attempts to obtain the effective molecular partition functions from second virial 
or cluster integral calculations, a classical type calculation may be simplest and thus might 
most naturally be used, with potential functions based on empirical spectroscopic constants 
and other known facts. Corrections for quantization effects in such cases would also be needed. 
A study of the effects of the change of spacing of levels when the quantum of action, h, is 
considered as approaching zero, as discussed in Appendix 1 (section 5.1) shows that the quantum 
effective partition function (referred to the lowest quantized level) may be obtained approxi- 
mately at high temperature by multiplying the classical effective partition function (referred 


to the classical zero energy) by the correction factor, 
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where for nonlinear polyatomic molecules 
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while for linear molecules 
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The constants in these expressions follow the usual spectroscopic notation with vibrational 


energy given by 
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with rotational constants 
: / 
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for the nonlinear molecule, and 


d 
B B > ( r 9 ) 
for the linear molecule. The diatomic molecule may be regarded as a special case of the 


linear molecule. The terms purely in w, in the bracket in eq (2.55) amount to a simple exten- 
sion of a result which Herzfeld and Teller [23] obtained by the method of Wigner and Kirkwood, 
and which others have used in other forms [24, 25]. Except for form, the terms linear in 
he/kT agree exactly with the result of Mayer and Mayer [26] for diatomic molecules. It will 
be noticed that if the two kinds of partition functions, quantum and classical, are to be 
referred to the same zero of energy, Say the classical zero or the zero of the potential energy 
function, then the additional required factor exp [| he/kT) Go|, where Gy is the zero point 
energy, removes the factor of exp [(he/k7)S(w'd,/2)] of eq (2.35) and replaces it with an expo- 
nential factor based on the much smaller anharmonicity constants. Thus, for such a quan- 
tization correction, the zero point energy might be considered as not being involved. The 
detailed procedure of taking a calculation with the zero of energy with the cluster particles 
dispersed and adjusting to a zero of energy at the lowest quantized level for the cluster to 
conform with custom in the tabulation of molecular thermal functions involves, in addition 
to eq (2.35), the use of the factor exp [(—e/k 7) DJ =exp [(he/kT) (Dy + G)|. This gives as the 
combined correction the expression in the bracket in eq (2.35), multiplied by exp [(—he/kT) Dy 
exp [(he/kT)(X\wd,/2}—G)|. In the adoption of eq (2.35) as the approximate quantization 
correction, not only has the existence of other than Boltzmann statistics been ignored as of 
no importance in calculations at high temperature and up to moderate densities, but also the 
quantization correction contributions associated with fragments of the cluster have also been 
disregarded, since at moderate temperatures the net quantization effects due to the fragments 
would be very small because of reduction due to powerful Boltzmann factors. 

The molar thermal functions for molecular-like clusters according to this formulation are 


to be consistent with 


A k7T1nQ (2.36) 


as in eq (2.9), but with each effective internal partition function Q,* including its appropriate 
quantization effect given by eq (2.55 and also multiplied by eXPp | he/kT)D, so as to be based 
on the usual molecular zero of energy. 

In eq (2.17) it is indicated that the effects of interactions between two different clusters 
or between a cluster and a particle may be indicated formally as a cluster integral for the 
joining of the two parts into a combined cluster. In making use of this principle as a guide 
in obtaining approximate thermodynamic functions for a gas mixture, it seems appropriate 
to emphasize its indication of the permissible arbitrariness of assignment of regions of phase 
space either completely to the formation of the combined cluster (in the style of the so-called 
mathematical cluster), or partly to the combined cluster (now possibly called a physical 
cluster), and partly to a second virial type interaction between the two parts of the cluster. 

One application of this principle is in regard to the effective internal partition function 
for the interaction of two atoms forming a diatomic molecule. One alternative is to use the 
entire effect in this partition function and the corresponding thermodynamic functions. On 
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the other hand, any part of this partition function may be used in obtaining the thermody- 
namic functions for the constituent, provided that the remainder of the effective partition 
function is retained as contributing in its appropriate way to the second virial type interaction 
between the corresponding atomic constituents. 

Somewhat similarly, the combination of an added outer electron with a positive atomic 
ion to form an atom or a positive atomic ion with one unit less of charge may also be treated 
in good approximation using an arbitrary subdivision of phase space. In a gas in which 
considerable excitation of levels near ionization occurs, there is also considerable formation of 
free ions, for which the effective interaction potential in the assembly as a whole may be taken 
as the Debye screened Coulomb form. Meeron [15] has shown that the nonquantum- 
mechanical treatment of Mayer [14] for ionic solutions can be given this simple interpretation, 
using quantities corresponding to second virial coefficients, even though expansion of higher 
order terms would not be the equivalent of a higher virial evaluation using the Debye screened 
Coulomb potential. Application to a high-temperature gas involves details of quantized 
states of the combined atom with ionic screening to which the approximate treatment by 
Ecker and Weizel [16] is relevant. Further discussion of this subject may be presented in 
a later paper. It would appear that a conventional sum of states for the atom up to some 
moderate principal quantum number of the outer electron should be acceptable provided that 
the remainder of the effective partition function be used in obtaining corrections for equation 
of state quantities. A very different division of phase space has been used in the treatment 
which will be presented in the following section and in the next paper of the present series. 

As the virial coefficients based on additive pair interactions as used by Mayer [14] in 
deriving quantities for his theory of ionic solutions are consistent with the cluster integrals 
for the Coulombic interaction, since this is additive, it is possible to treat the effect of ioni- 
zation using Maver’s classical ion-solution theory, which he arranged primarily for an assembly 
of ions with each pair separated by a distance greater than its particular sphere radius. The 
interior of each ion sphere is a region actually accessible to the electrons of the system, as they 
can enter and become energetically attached so as to form ions of various numbers of electrons. 
As each such group is a cluster of particles, one would wish to evaluate contributions to the 
Ursell integrals for electrons within the sphere but with such quantization corrections as 
would be needed to approximate the effects due to the discrete states of the ions. Unfortu- 
nately for purposes of evaluation, the integrand of the general integral for the many-particle 
case is of such complexity and the integrals are so highly multiple that precise direct evalua- 
tion according to this procedure appears impossible. It is therefore appropriate to make use 
of convenient approximations. The effective partition function may be estimated approxi- 
mately using the observed and estimated low lying quantized states with classically estimated 
contributions from higher energy regions of phase space within the assigned radii. These 
may depend on the representation of effective potential-energy functions, with simplifying 
approximations made possible by partial compensation for cognate fragment terms in the 
cluster integrals. 

Because the thermodynamic functions obtained are specifically arranged for use with an 
ionic-gas treatment based on Mayer’s development of ion-solution theory with specific values 
for the ion radii, it seems appropriate now to summarize the results of his approximate 
treatment when applied to such an ionic-gas mixture in chemical equilibrium. 


3. Ionized Gas as an Ionic Solution 


The nonquantum-mechanical theory of ionic solutions given by Mayer [14] retains terms 
applving to ions free to move to any accessible pair separation greater than their collision 
radius. By also including the contribution of Mayer’s eq (46) due to the departure of the 
non-Coulombic pair potential from rigid-sphere form, an approximate expression consistent 
with the terms retained in his original derivation may be obtained. On this basis, the 
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Helmholtz free energy for the gaseous assembly of molecules, atoms, atomic ions, and electrons 
with the accompanying thermal radiation may be taken approximately as 
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where C’, is the C, representing electron concentration and where terms are given for the 
first-order quantum corrections for an electron gas including Fermi-Dirac statistics [27, 28] for 
electrons. 
The svmbols used in the present ion-solution equations have the following meanings: 
(,=concentration relative to that of a pure ideal gas at standard conditions, 1. e., C 
N,Vo/NoV, 
Z,=ionic charge of each ion of type s, in protonic units, 
e=electronic charge, 
D=effective dielectric constant, ~ unity for low density, 
k= Boltzmann constant, 
as a,,.DkT €, 
a,,=collision radius for ion s and ion r. 
b.(¢), a(¢), and h,(¢) are functions defined and tabulated bv Poirier [29]. The argument of 
these functions is defined by 
29a, (D> IZ2C,) 
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where x tre DZ DkT. as known from the Debyve-Hiickel theory. 

The omission of 6;, g;, and h, (A, appears in appendix 2) on the basis of net electrical 
neutrality for the complete assembly Is acceptable in case there is but one value of a,,, the case 
for which Poirier’s tables were prepared. If there are several values of a,,, then cancellation 
does not occur in general for ion separations between the least and greatest a,,. Evaluation 
of the net contribution may then be based on the explicit expressions 
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A(—A/RT) has been included to account for any further corrections omitted from the equation 
of state. As pointed out in the preceding section, nonionic equation-of-state effects at low and 
moderate density may be treated by cluster theory and could be accounted for by including all 
clusters as separate new constituents in the first term in the brackets in eq (3.1 Alternatively, 
the effects may be indicated in accord with the usual procedure emploved to obtain thermo- 
dynamic corrections for nonideality. Thus, one might use 


y VT(B—B Cnt D—D 
at —_ 5° 24 ‘ge My >°C)4 - oe ef 3.2 
4(—pr) Al y, (= ay ‘aut gye (Sut ] 


where B— By, C—Cy, and D— VD, represent additional parts of the virial coefficients required 
over and above the zeroth order values implied by the remainder of the ion-solution-cluster 
treatment being used. Corrections to a nonvirial type representation could similarly be 
provided. 
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The quantity /,, in eq (3.1) has the value 


» » aa 2. an 
exp {[ Y +p (1—< ier}. 


for O R a and the value 
, £,b,e ' -_ LL,€ 
exp4[ Usetap (l—e**) [kT > —exp| —p 5 ape" 
LL Di DR,kT 
for a R In the present method of describing effects between ion cores and an outer 
electron, h,, will be taken as zero in the range 0< R,,<a,,, since this is a region which con- 
tributes to the internal sum of states determining the function (—A°/RT7), for the corresponding 


constituent of ion core with outer electron, and the effect is taken into account by that means. 
If a,, is well chosen for the approach of two atomic ions, A,, can be small in almost all of both 
regions, as pointed out by Maver for his closely related &,,. Thus for R<a,,, h,, is small 
because the pair potential (’,, is very large; while for R>a,,, it is small because U’,, is almost 
identical with Z,Z,e/DR, 

If the contribution within the braces of eq (3.1) due to A,, for the interaction of ions of 
like charge is to be based on the use of a Debye screened Coulomb repulsion for R,,<a,, as 


well as for PR a,,, it may be given as 
2 Ne een 1 , 
2 ay, |” 22a! ( 4, Ler) Der; 
with 


I Z,4,¢7/Da,,kT=Z,Z,| >. 


Formulas for qg,(2) are given in Appendix 3 (section 5.3). 

Equation (3.1) has been written in rather general form which allows it to be used in more 
than one way. An approximate treatment for an assembly including atomic ions under condi- 
tions of high-ion density may be obtained by choosing a single a,, to apply throughout or a very 
small number of different a,,’s, with suitably limited sums to represent the internal partition 
functions. For conditions of relatively low density, however, it may be desirable to use 
internal partition functions for the atomic ions based on an extensive set of levels extending to 
moderately large electron separations. The so-called “collision radius’? between an electron 
and an atomic ion could in principle be taken correspondingly large in this case, but with the 
added detail that the v2, or initial index in the summation in eq (3.1) and in similar expressions 
in Appendix 2 (section 5.3) would be greater than zero, signifving free entry into the ion sphere 
for this one case rather than having a zero value signifying the exclusion from the ion sphere 
such as is used in regard to the close approach of two atomic ions. 

Calculations concerning the concentration of the various constituents of the gas mixture 
under conditions of equilibrium can be made based on this ion-solution treatment according 
to well known principles. Thus, we may consider a stoichiometric chemical reaction among 
substances }, with stoichiometric coefficients ny. and ny.» in the reaction equation (c.f., 


eq (2.27) 


Done Vena Vy 
which may be written as 
> nes -=0 (3.3) 
with n,, representing both nj,» and —nyzy. The equation for chemical equilibrium for this 
kth reaction equation is 
> Nest, =0 (3.4) 


in terms of the chemical potentials u,, with 
us=(OA/ON,)7 y. (3.5) 
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This gives 


Ki, =T1Cy#=( 7) exp| Sm ( — rr) a7 | Hy, = Kt 3.6 

where AH 4 is >)n;,, Hy for the reaction with 
: us SOZzC,)* (SoZ?) EN ] = 
y exp V: DkT) Od 


according to the Debye-Hiickel theory, and with 
- l = 2 Wy "Jyl Zi ] y } 
Y2=CXP 9 TSA Ga) Va Qa da(— lV)" Lia, Nislby (Ger) + ole 


(222M) ame a as . 

~(S3Z?C,)? heh hana nea ' ICC, ZZ ia Ju\Psr) 0-9) 
as the principal correction based on the tabulated functions representing Mayer’s extension 
of the theory. A correction for the effect of the A,, terms may be expressed by 


’ ya 4 - 
a XD oN — hi No F ] Dar f?? IP . os ) 1S \/ ( N, Oh Ie PR? IR 
Ya OND A cals rT en a SSAC, Ae 8 Bl Dies 
. >. 9) 
which for Debye screened Coulomb repulsion between ions takes the form 
2 
¥3—— exp 2 é 7 2252249 re oon 
1 We No 1/ it \*N BQ 
2 TR asap) 7 Cote tH] seta Nt et | Zim | f- 3-9) 
12a mke*T? V, = 32\amk7T Vo 2(>Z:C,) “a 


An estimate for effect of dependence of dielectric constant J) on composition may be made 
using 


T (S°772C € N, | “e . 
% c — — = S=, ie 
. x {[ DET) \V,) ~2(S3Z3C) By COZ 


- 0 


d ) 
[vb, Yy y—2)grl¢ | | = _ }}: 3.10 


although a more detailed microscopic analysis of the kinds of screening affecting both y, and 
vy, should be made. In fact, the rapid increase of atomic polarizability with principal quantum 
number poses another convergence problem unless each atomic state be taken as defining an 


independent constituent. Additional effects related to the equation of state may be estimated 
from. 
. exp >on (0/0¢ B I > 56 V 

C—O) (Se 2V D—D,)(>5C,)*/3V .ac 3.1] 


It is also necessary to take into account the dependence of energies of states and the heats of 
formation and reaction on the ion density. Some indication of the dependence for atomic 
ions mea) be obtained from i study of atomic levels for a Debve screened atomic ion. 
Discussion of these effects will be deferred until a later publication. 

The equations for “mass’” and charge balance are well known and obvious but may be 
represented here for completeness Thus if the number of atoms of element -Y in a molecule 
of constituent J)’, is /,,, then 
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is the local average concentration of this element, usable either directly or in ratio to other 
similar quantities in the expression of conservation of each chemical element. Similarly, if 
the electric charge on each molecule of constituent Y, is Z, in protonie units, then 


20 2:=0 (3.13) 


states that, on the average the gas is electrically neutral. 

The concentrations of constituents as given by these equations may be termed nominal 
concentrations of clusters considered as molecules. To obtain the actual molecular concentra- 
tions C;” from the corresponding nominal cluster concentrations C,, one may make use of 


the ratio 


CH 1C.= QS” |G = 359: exp e,/kT)/Q3", (3.14) 
with the partition function 
Qs” =>049; exp é/kT), (3.15) 


I 


formed by a sum over whatever group of states is deemed to comprise the specification of the 
molecular form (or similarly for atoms and atomic ions). The concentration of molecules in 
their ground states is given by 


CJo EXp (—e/kT)/Q, (3.16) 


where go is the @ priort weight of the ground-state level and 6 is its energy relative to the 
energy zero used in obtaining Q,°". 

It is to be emphasized that nominal concentrations of clusters are in general not identical 
with concentrations of molecules or ions bearing the equivalent name. If accurate estimates 


for concentrations of these species are required, then (3.14) and (3.15) must be used. 


4. Discussion 


The calculation of thermodynamic properties of gas mixtures in chemical equilibrinm 
at high temperature and moderate density, with each molecule, atom, or ion considered as 
a cluster of its constituent particles, has been discussed. The method of calculation may be 
used for very high temperatures at which the number density of ions is high, using ion-solution 
theory in Mavyer’s formulation. 

The detailed discussion of the method of calculation of thermodynamic functions for the 
atoms and atomic ions that form the principal constituents of the gas mixture at high tem- 
perature is reserved for paper II of this series. Under conditions of considerably elevated 
density, the proportion of small molecules may be enough to make accurately descriptive 
functions desirable, even at rather high temperatures. Methods of calculating such functions 
for diatomic molecules may be discussed further in paper IIT of this series, including both con- 
ventional partition functions and the effective partition functions provided by cluster in- 
tegrals. Molecular cluster integrals automatically include distortion and stretching effects 
due to centrifugal forces and so are especially suited to calculation at high temperature. The 
effect of quantization on the molecular cluster integrals has been included in considerable de- 
tail, vielding estimated “effective partition functions.” Ratios between these functions give 
proper equilibrium constants for the treatment of a gas mixture as a mixture of ideal gases 
in the absence of ionization for the purpose of representing the over-all thermodynamic 
properties of the mixture. 

Special emphasis can be given to the concept of the “effective partition function” con- 
sidered as a convenient function for expressing the increase of the partition function of the 
entire gas due to the joining together of the parts of the cluster. It is thought that this con- 
cept is useful in clarifying the physical meaning of the kind of cluster which has sometimes been 


called a ‘mathematical cluster.” 
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The analysis of the gas-mixture properties by means of equilibrium constants for cluster 
formation, or the cluster integrals for clusters of all compositions for small numbers of par- 
ticles is formally no less complete than the representation of the gas properties according 
to a virial series expansion carried through virials of an equivalent order, or involving inter- 
actions of an identical number of particles as in the highest virial coefficient included. Indeed, 
in the low density region the cluster integral representation as a treatment of chemical 
equilibria is considerably the more complete when applied to effects of chemical bond for- 
mation under conditions that favor the formation of small clusters (molecules) only. This 
detail of good representation may be touched on again in connection with paper III of this 
series, dealing with diatomic molecules. If triatomic molecules were the largest molecules 
formed, then a cluster treatment including all clusters of three atoms would include all 
chemical bonding effects varying directly exactly as the third power of the free atom con- 
centration and would include higher order virial type effects due to chemical bonds in pairs 
and triples of atoms. In the high-density region, however, the higher virial coefficients have 
large excluded volume effects which are not included in a treatment covering only pairs and 
triples of atom particles as clusters. For these, the cluster treatment would need to be 
carried through for clusters of the same number of particles as the numerical order of the virial 
coefficient in order to give a complete representation for the given coefficient. Application In 
the high-density region is complicated further by the detail that for an ordinary chemically 
inert gas, the virial expansion itself, to which the cluster treatment is formally equivalent, 
is found to be progressively poorer as a representation when extended to higher and higher 
density. This might suggest that the virial expansion itself may be correct only asymptotically 
in the limit of low density. Nevertheless, the virial expansion is to be regarded as a useful 
approximation capable in principle of giving good estimates up into regions of considerable 
density. As the equivalent cluster treatment gives a still less suitable extension to high density, 
the problem of practical estimation of properties of chemically reacting gas mixtures at high 
density and temperature calls for further estimates of volume effeets bevond accessible cluster 
evaluation, so as to be either in accord with empirical and relevant theoretical knowledge of 
higher virial coefficients or in accord with gas theory of still broader application. 


5. Appendices 
5.1. Appendix 1 


For high temperature calculations of thermodynamic functions for molecules according 
to the present cluster approach, one might first obtain the ordinary nonquantum cluster in- 
tegrals based on potential functions agreeing with empirical spectroscopic constants and 
other known facts, following which, one would wish to combine quantum corrections based 
on known molecular spectroscopic constants. By taking the ratio between the ordinary 
partition function for the molecule and a limiting partition function consistent with it, ob- 
tained by reducing the spacing between levels uniformly to zero with a proportionate reduc- 
tion in weight per level, one may readily obtain a principal correcting factor for the effect of 
quantization. Since the fragment submolecules implied in the cluster integral occur at much 
higher energy than the combined molecule, only a small error is introduced by ignoring the 
quantization corrections for these fragments. 

In a well-known system of notation, the logarithm of the ordinary partition function, 
with energy referred to the ground-state level of a molecule, may be written as 


4, he ‘ 
In@ InQrt pr 2a4 ln l—e “i —( in) d.(d,+1)e72*i(1—e7*s) -2 
th 2,,d,d,e7%e- "(1 —e-*) “(1 —e-4) 84-25% 9,,e7%«(1—e * ft aad — 
Al.1) 


Here Qe is the rotational partition function for the ground vibrational state in classical 
approximation, @; represents the Stripp-Kirkwood [30] correction for nonlinear molecules or 
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the Mulholland correction for linear molecules due to the finite spacing of the rotational levels. 
The variable u,; may be expressed in terms of the usual vibrational constants of 


G= Vw, (v,+d,/2) + S524 (0; +d,/2)?+ > 35> 52:4; (0; +d,/2) (v,+d,/2) + SSgiul? (A1.2) 


according to 


ul * w, + (d,;+1)2;,4 32d Zi |} (A1.3) 


In view of the principal vibrational dependence of the rotational partition function, repre- 
sentable by a factor exp >3a,v;, In@e may be replaced by InQ@g.+ 34> 5d,a;, where Qe, refers 


i 
to the equilibrium configuration. The derivation of the main part of the vibrational effect 
proceeds as already stated with a factor \ for the level spacing appearing wherever a quantum 
number factor (7,+-d,/2) appears in the fundamental expressions: 


he . 
InQ@ir > 5d, ind S 3d, ln | 1—eo** ]- -( ip)™* > rid, (dit le~™™™ (1—e-"*™) -2 
> tf se IM e 4 (1 —e-™) "(1 —e -9™) — + 25> ge “™ (L—e i ™) -2 
| ASS¢ d le ur (] ¢ use ) sl 5] (A1.4) 
with 
he . . 
u(r ay (Nw, +A? (d+ Daye t+ 3D d 24). (A1.5) 
ke] ji 
The quantity InQ(A=1)—InQ(A=0) gives the vibrational quantization correction, to be 


put in suitable form by expansion of the various functions, either directly or giving the form 


l he 


VY(A= 1) /QGA=0 {1 T19 pp laud (d,+l)2r;,4 D2 d jty(3 ww; '—wjw 
1 (/he\? l he\' ;, 
Be ret. 2 = S“‘d.w S*(5d,+2 » 4 
Sart 23 ul a4 ( gp) Bade? + 5760 (jer) (a Cait 2d 


‘ . » 9 he ral, 
ODM aia] bexp Pe (ALG) 
v i 


The rotational quantization correction is already known, being given by the factor 1+-6,/T 
as indicated earlier above. 

It may be noted that the expression is in such form as to change from a classical treatment 
with energy measured from the classical 7,=0 or potential minimum value to a quantum 


treatment with energy measured from the quantum v,;=0 or ground-state energy. 
5.2. Appendix 2. 


Additional thermodynamic functions and relations for the equilibrium gas mixture as 


derived from eq (3.1) are: 
PV Vf 1 2(S5Z20,) 6 ( No\' 
i . ae ¥ —_ 
RTO V3 DRT) (7) cml 
N ‘ Kk Oh l he N 1 h? m N, 
: ee leet 2 24 ? - — — 7 — . Cz 
cae! sf VY | (/ 2 OK ) an lti,dt E mke’T? V, 35(=neT) | 
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The dependence of the dielectric constant P) on the gas density using the Clausius- 
Mosotti relation provides as a further contribution? to PV/RT, the quantity 


D+-2)(D—1) (6E 
- D ( RT) sons 


? 


where (6//?7),.,,, includes the ordinary ionic terms in eq (A2.2), 1. e., the limiting law Debve-. 
toni { ve 


Hiickel term and the Mayer-Poirier terms for which the dependence on PD is of the same form 
as for 7. 
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For such ion interactions as are taken to have identical radii a@,,, the factors >)>5C,C,Z;Z? 
& Tr 
in (3.1) and (A2.1) to (A2.6) can be grouped together as (20.2; )( p> Lr), reducing further 
s Tr 


to the form used by Poirier if all are identical. 


5.3. Appendix 3 


Functions relevant to the Coulombic interaction in the close approach region for ions of 


like sign include a contribution to —A/RT due to h,, given by 
V 2 N, . ‘ \ ‘ 
: - S 2 Uae SO T lor, (A3.1) 
Vo 3 Vo Os 1=0 
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The related contribution to PV/RT is 


A similar contribution to E/RT is given by 
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The contribution to H/RT is given by 
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with the ¢, as given above, and with rde,/dr obtainable by differentiation, giving relations 


such as 


rile 


rel ¢ 





da © qt(3: a 5a) 
da 3 u-(G? at if a )e 


343 l 123 523 _..,. J 
da war rd +( r in r’ = r® ye 


6. References 


fl] H. D. Ursell, Proce. Cambridge Phil. Soc. 23, 685 (1927 


2} B. Kahn a 


13] B. Kahn, " 


4) J. Kilp: 
5] W. G. Me 
I6) K. Fuchs, 


17) ¥ Muto, 


nd G. Kk. Uhlenbeck, Physica §, 399 (1938 

Phesis, Utrecht (1938 

atrick, J. Chem. Phys. 21, 274 (1953 

Millan and J. Ek. Maver, J. Chem. Phys. 13, 276 (1945 
Proc. Rov. Sor London) [A] 179, 408 (1942 


J. Phys. Soc. J pan 4, 23% 238, 242 (1949 


[8] H. C. Urey, Astrophys. J. 59, 1 (1924 


19] Kk. Fermi, 

[10] R. H. Fow 
L936 

[11] R. Vio Me 


Propulsi 


[12] M. Planck, 


[13] D. Bohn a 
[14] J. bk. Mave 
115] Kk. Meeron 
16] G keke a 
(17) T. L. Hill 
(IS) H. W. Woe 
119) P Rosen, 
[20] B. M. Axil 


[21] L, Janse ni 


Z. Physik 26, 54 (1924 


le r, Statistical NI chanics 2d | litic Nn The \I icmillan Co London, engl . 


grebliar Dissertation, (Calif. Inst. Technol 1953 or Guggenheim Jet 
on Center Tech. Report No. 7, (Calif. Inst. Technol., Nov. 1952 
Ann. Phys. [4 75. 673 1024 
nd D. Pines, Phys. Rev. 88, 338 (1952): 92, 609 (1953 
r, J. Chem. Phys. 18, 1426 (1950 
J. Chem. Phys. 26, 804 (1957 
nd W. Weizel, Ann. Phys. [6] 17, 126 (1956 
J. Chem. Phys 23. 617 1055 
Hlev, J. Cher Phys , # 236 (1953 


J. Chem. Phys. 21, 1007 (1953 


rod, J. Chem, Phys. 19, 719 (1951 
Dissert iio} Leider 1055 see also. L Janse! ind Z l Slawsky, J (Chem, 


Phys. 22, 1701 (1954 


[22] T Kihara 
Japan 1 
11, 1050 

\. F. Herz 
Bigeleise 
sigeleise 
| \I ive 
York, N 


Rev. Mod. Phys. 27, 412 (1955): Y. Midzuno and T. Kihara, J. Phys. Soe. 

1, 1045 (1956); S. Koba, S. Kaneko, and T. Kihara, J. Phys. Soc. Japan 
LO56 

feld and Ek. Teller, Phys. Rev. 54, 912 (1938 

nand M. G. Maver, J. Chem. Phys. 15, 261 (1947 

n, J. Chem. Phys. 23, 2264 (1955 

r al ad \I (; \I iver, St itistical Mechat cs Jd hn Wile v & Sons Ine . New 
Y., 1940 


[27 R. J. Plock, Dissertatio1 Yale L957 

[28] Kk. C. Stoner, Phil. Mag. [7] 21, 145 (1936): [7] 28, 257 (1939 

124 J. Poiri - 4 nem P} vs 21, 65 Q72 1953 

150} Ix | Stripp ind J (; IK rkwood, J ( hem Phvs 19, 1151 1951 


WASHINGT: 


wN, July 18, 1958 


490 





pi 
th 
tit 
pe 


m 





tions 





nal of Research of the National Bureau of Standards Vol. 61, No. 6, December 1958 Research Paper 2917 





Parallel Testing Interferometer 


James B. Saunders 


The conventional methods of testing the parallelism of opaque bodies, such as gage 
blocks, by interferometry require wringing of the body to an optical flat. This operation 

disturbs the temperature equilibrium, necessitating long periods between tests, especially 
for long blocks. It often injures the surfaces of both the optical flat and the test body. 
\lso, if the body is a standard gage block, repeated wringings during use ultimately change 
the dimension. This paper describes an interferometer for measuring the parallelism of 
gage blocks and other bodies of any reasonable length without the necessity of the wringing 
operation. Two forms of this instrument are used—one for testing long blocks and another 
for testing short blocks. Either form can be constructed for testing blocks of any length, 
but two forms are found to be more practical. 


1. Introduction 





The conventional procedures for measuring the 
parallelism of gage blocks * ? require the wringing of 
the blocks onto an optical flat. The wringing opera- 
tion often injures the contacted surfaces and re- 
peated wringings necessitate frequent refinishing of 
the optical flat that is used as a base. A method for 
measuring parallelism, without the wringing opera- 

| tion, has significant advantages because the danger 
of injury the contacted surfaces is eliminated. 
Accordingly, two instruments that utilize this 
method are described: one is for testing very long 
blocks and similar whereas the other is 
designed for short Both instruments 
low orders of interference and neither requires the 
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Késters double-image prism. 


Figure 1. 


use of a standard. 
A small angle is formed by the two 30° edges shown exaggerated in B. 


2. Optics of the Interferometer 
becomes hard by cooling after adjustments are 


complete. Furthermore, the tilting of any plane 
surface outside the prism that affects the two com- 
ponent beams between division and recombination, 


A description of the optics for either form of this 
instrument covers a large portion of that for the 





other. For distinction we will designate them as id recom! 
“the long-block interferometer” and ‘the short- | will produce equal effects in this direction and, 
block interferometer’. The double-image prism consequently, will not affect the fringe width. 


However, the rotation of plane surfaces about an 
axis normal to the plane of figure 1A will produce 
equal effects on the two component beams but in 


used im these instruments is adjusted during con- 
struction * so that a ray of light, shown in the plane 
of figure LA, after division into two component rays, 


1 and 2, at P, will, on reflection at P,; and P», lie in 
planes that are parallel to the semireflecting plane of 
the prism but deviate equally toward or from oppo- 
site sides of the plane of figure 1A. The projection 
of the light ravs on the dividing plane is shown in 
figure 1B, which is perpendicular to the plane of 
figure 1A. This deviation is effected by rotating 
one component of the prism relative to the other 
about an axis normal to the dividing plane of the 
prism. If this deviation is held constant, the width 
ofjthe interference fringes in the direction normal to 
the plane of figure 1A is fixed. This component of 
fringe width is, therefore, frozen into the system * 
when the cement between the component prisms 
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opposite directions, thus producing a proportionate 
effect on the fringe width in the direction parallel to 
this plane. 

Because of the above-described properties of this 
prism, adjustments of the instrument in which it is 
used affect the fringe width in one direction only. 
Consequently, when measuring the parallelism of 
gage blocks the test can be applied to parallelism in 
only one direction at atime. To test for parallelism 
in other directions the block must be rotated. 


2.1. Long-Block Interferometer 


Figure 2A is a horizontal section through the 
optical elements of this instrument. The light from 
a source at 8,(S.) is collimated by lens L,(L.) and 
divided into two equal components by the beam- 
dividing plane B,(B,). Each component suffers total 
internal reflection in the prism and emerges in planes 
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parallel to B,(B,), but at a small angle to the 
ere C.(C,). T he latter condition is obtained by 

i slight rotation of P,(P,) about an axis normal to 
the div iding plane B,(B,). This serves to permit the 
elimination of light reflected from C,(C, 

The two prisms, P,; and P,, are separated by a dis- 
tance exceeding the length of the longest block to be 
tested. If desired, this distance may be made 
adjustable. The two dividing planes, B, and Bz, are 
adjusted to be coplanar and the two base faces, C;, 
and C,, are adjusted parallel to each other. <A line 
joming the centers of P, and P, is adjusted to form a 
small angle with the normal to faces C, and C, 

The light emerges from prism P,(P,) as two sepa- 
rated components, one on each side of the dividing 
plane. It enters prism P,(P,), again suffers total 
internal reflection, and each pair of component rays 
recombines in the plane of B,(B,). One-half of each 
beam proceeds to the neighborhood of source S,(S 
and the other half to S.’(S,’).. An observer at 8,’(S 
sees a set of interference fringes that cover the entire 
aperture. 

If a gage block, G, is inserted in the position shown 
with its end face, G,, adjusted normal to the light 
beams, it will reflect equal and corresponding parts 
of the two component light beams from 8S, back 
through P; to S,’. Accordingly, the observer at 5,’ 
sees a background set of fringes, produced by light 
from S, and another set on the face of G,, that is 
produced by light from S,. Since G, is normal 
the light beams these two fringes will be 
parallel to each other and to the plane of figure 2A 


, 


sets of 


bs ase 


If the other end, G., of the gage block is parallel to 
G,, the fringes seen at 5,’ will likewise be parallel to 
each other and to the plane of figure 2A. If, how- 
ever, G, is not parallel to G,, in the plane of figure 
2A, it will not be normal to the light beams, and the 
sets of frmges seen at 5,’ will not be parallel to each 
other The angle between these two sets of fringes 
is a measure of the angle between G, and Gz, in the 
plane of figure 2A. The component of the angle 
between G, and G, that is perpendicular to the plane 
of figure 2A (or horizontally in figs. 2B to 2G) does 
not affect the fringes because it affects all pairs of 
component beams equally. If the component of 
the angle between the gage-block surfaces that is 
normal to the plane of figure 2A (or vertically in 
figs. 2B to 2G) ts desired, the block must be rotated 
90° and the operation repeated. 

Since each prism is adjusted for complete com- 
pensation in the plane of figure 2A, white light can 
be used. A measure of the vertical width of the 
fringes (perpendicular to fig. 2A) for a known mono- 
chromatic light, with a micrometer eyepiece at 
S.’, gives a calibration of the micrometer scale in 
units (microns, millionths of an inch, etc.) of length 
for measuring the displacement of white-light fringes 
from a reference port on the cvage-block 
surface. 

The procedure for adjusting a block 1s ex- 
plained with the aid of inserts in figure 2. In general, 
when the block is placed on its supports, the light 
reflected from its end will not the observer 
because of excessive angular d from the 
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evepiece. The block will appear in silhouette, as 
indicated im figure 2B. When the surface G, is 
adjusted approximately normal to the light, fine 
fringes will usually be visible in the area covered by 
both images of it, as shown in figure 2C. The 
images of the two parts of G, are made to coincide by 
rotating G about the center line of figure 2A, causing 
the image to change from that of figure 2C to figure 
9D. A lateral motion, without rotation, will then 
change the image from that shown in figure 2D to 
foure 2E. <A further small rotation of G about an 
axis normal to figure 2A brings G, normal to the light 
and the fringes on G, will appear horizor tal and 
parallel to the background fringes as shown in figure 
oF, The observer then moves to position S,’ and 
observes the set of fringes shown tn figure 2G. The 
angle between these two sets of fringes corresponds 
to the angle between surfaces G, and Gy). 

A photograph of the long-block interferometer is 
shown in figure 3. The base of the instrument is 
designed for rigidity so as to avoid flexure. Rigidity 
is quite important. The two prisms are mounted on 
rigid tables at each end of the base. Three screws 
D, and two others not shown) permit raising and 
lowering of prism P,. They also permit rotation 
or tilting of this prism about any chosen horizontal 
axis. The prism housings are fastened to the table 
tops by means of large-headed screws, in oversize 
holes, that permit lateral adjustments of the prisms 
relative to each other. Two screws, D, and another 
concealed by the housing of prism P,, permit small 
rotations of this prism about a vertical axis by apply- 
ing lateral torques to the legs of the table. Similarly, 
screws D, and D, permit rotary adjustments of the 
other prism housing 

The adjustments described above permit the 
alinement of the two prisms. This adjustment is 
critical, rather difficult to attain, but when once 
obtained is very stable. The final adjustments are 
executed while observing interference fringes pro- 
duced by lens L’,, in one of the many images of the 


source S,. This author uses a pinhole source of 
approximately 1-mm diam and a short-focus lens to 
observe its image at S,’. There are two sets of 
fringes, superimposed upon each other, in the proper 
image to be used. When these two sets of fringes are 
horizontal and very broad, fringes can be seen with 
the two eyepieces when focused on the pinhole. 
These fringes, in white light, are horizontal with the 
zero order in the center of the field. 

The supports on which the long blocks rest are 
located at the Airy points so as to reduce changes in 
the angle between G, and G, due to gravitational 
distortion.’ These supports rest on an adjustable 
plate which, in turn, is supported at one end by two 
points and at the other by one point. This plate is 
adjustable at one end, laterally with screw D, 
and vertically with another screw D, which is con- 


cealed in figure 3. 
2.2. Short-Block Interferometer 


The optics of the short-block interferometer are 
shown in figure 4. The double-image prism, lenses 
light source, and viewing position in figure 4B are 
identical to that of either end of the long-block 
instrument described above. A_ reflecting prism, 
P, in figure 4, replaces one of the prism assemblies 
of the long-block instrument. Also, the optical 
axis of the instrument is vertical instead of hori- 
zontal. Figure 5 is a photograph of the short-block 
interferometer. Figure 4A may be considered a 
section through the center of 4B, coincident with the 
dividing plane of prism P. The indicated rays, 1 
and 2 in figure 4A, however, do not lie in this plane. 
Their positions relative to it are indicated in 4E, 
which is a vertical view through 4B. The two 
surfaces, G, and Go, of the gage block (figs. 4A and 
4D) appear as G and G’ in figure 4E. G’, in figure 
4A is an image of G, as seen by light reflected from 
the right-angle prism, PR. 


F. H. Rolt, Gauges and fine measurements II, 340 (Macmillan and Co., 1929). 





Figure 3. Photograph of the long-block interferometer. 
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Figures B, C, and D are sections through A Figure 
fringes tor a parallel and a nonparallel block, respectively. 
A plane optical wedge, W,, and its image, W’,, 


(in figure 4A) are shown between G and its image G’. 
Figures 4B, 4C, and 4D are sections through 4A, 
along the lines indicated. The rage block is not 
located in the center as was the case in the instru- 
ment described above. 

The two component rays of light, 1 and 1’ (figs. 
4A and 4D), are caused to reflect normally from G, 
by adjusting G with a leveling screw. They return 
into P where they recombine to produce the inter- 
ference fringes, F;, shown in figure 4F. The two 
component rays, 2 and 2’ (figs. 4A and 4C) are 
transmitted downward through the optical wedge, 
W;, suffer two internal reflections in Pp, and if Pp, 
is properly adjusted, return upward and _ parallel 
to their directions of incidence through W, to Go. 

The wedge, W,, is adjusted initially by rotation 
so that its thickness is constant at all points In either 
of the planes B, C, and D. When in this neutral 
position it does not affect the interference fringes 
because of compensation in each pair of component 
beams that pass through it. The function of W, 
will be explained later. 
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In order to measure the angle between G, and 
G,, the deviation of the light by Pg toward or from 
the dividing plane of P must either be reduced to 
zero or its effect eliminated by measuring the ob- 
served angle for two orientations of G, which are 
180° apart. The light that is not intercepted by G 
forms an interference pattern of uniform tint (or 
color) that fills the background about and between 
the two images of the gage block, shown in figure 


{Ff or 4F’. Figures 4F and 4F’ represent the con- 
ditions observed when the ends of the block are 
parallel and nonparallel, respectively. A typical 


pair of component rays, which form this interference 
pattern, is indicated by 3 and 3’. They travel 
downward in figure 4C and upward in figure 4D. 
If the right-angle edve of Pe is normal to the divid- 
ing plane of P, the pair of ravs 3 and 3’ in figure 4C 
can be made to return in planes that are parallel 
to the dividing plane, by rotating Pg about an axis 
parallel to the plane of 4A and normal to the inel- 
dent light. This condition is attained when the 
background fringe is infinitely broad. The direction 
of the background fringes, when not infinitely broad, 
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| Fiat RE 5 Photograph ot the short-bloac A inte rferometer. 


remain parallel to the dividing plane, because for 
each pair of component rays, such as 3 and 3’, there 
‘isa corresponding pair, 4 and 4’, that travel identical 
paths but in opposite directions. The optical path 
‘differences are, therefore, equal to each other and 
| also equal to that for any other pair of component 
i in the plane of 4C and 4D. The order of mnter- 
ference along the dividing plane corresponds to the 
| optical path difference that was introduced into the 
double-image prism by the built-in wedges at the 
pomt that corresponds to the point of intersection 
of the right-angle edge of Pg and the dividing plane 
fP. This point is located in the center of figure 4E. 
In general, due to the inherent error of judging 
when the background frmges are mfinitely broad 
and to imperfections in the optical elements, a more 
precise method of evaluating the wedge between G, 
and G, is to measure the wedge for two positions 
that differ by 180 If the background fringes are 
maltered, the instrumental errors will be equal for 
jthe two positions and the value of the wedge un- 
changed except in sign. Consequently, the albebraic 
difference yields twice the value of the wedge. 


There are three ways that one might evaluate 
the wedge between G and G The first is to 
rotate Pp until the order of interference at points 


Cand E (fig. 6A) ar equal; then rotate the gage 
block until the orders at A and B are equal; and 
finally, observe the difference in order of interfer- 
at pomts F and H. The second method is: 
\fter performing the above operations, mstead of 
rading the order difference between F and H, re- 


ence 


| 
| 
| 
| 
} 







































AL | OIRECTION 4 
J OF WEDGE \ 
/ —> — 
| We-W2 Fy 
Cr) | Os 
| i 
a NG | WA 
A 8 c 


Ficure 6. Figures A and C represent the positions of reference 
points relative to the different sets of interference fringes. 


Figure B shows the wedge from which the two smaller wedges, W? and W”, 
were cut and the difference in optical thickness between them, 


duce this order difference to zero by rotating W, 
and read the resultant change on a scale attached 
to W,, figure 4E. This scale may be calibrated 
with monochromatic light and the units may be 
radians, degrees, or the corresponding variation in 
height of the block. A third method is to leave 
the wedge in its neutral position, adjust G so that 
the orders of interference of A and B are equal, and 
change the order at H to equal that at F by rotating 
P, about an axis normal to the incident light and 
parallel to the plane of figure 4A. The order of 
interference between two points such as C and D 
(fig. 6A) after rotating P, will be equal to one-half 
of that between F and H before this rotation was 
performed. By choosing a point, such as E im fig- 
ure 6, such that CE equals K times CD, the order 
difference between C and E will be K times that 
between C and D. 

When using this last method for testing gages 
that are almost parallel, the angle between G, and 
G, will be small and the background fringes will be 
too broad for reading fractions of fringes. To elim- 
inate this difficulty, an optical wedge, illustrated in 
figure 6B, is constructed and from it two sections 
W, and W,’ are cut and placed on I’, as shown in 
figure 7. The wedges W, and W,’ are equal, but 
when placed in the position shown, the effect is to 
narrow the background fringes seen through them. 
The results are illustrated in figure 6C. The differ- 
ence in thickness of W, and W,’ at a selected refer- 
ence point E (fig. 6C) is determined by the choice 
of the corresponding positions on the plate from 
which they were cut. This difference in thickness 
is chosen so as to cause the zero order of interference 
to pass through the chosen point when the back- 
eround fringes about W.-W,’ are infinitely broad. 

If the angles of wedges W, and W,’ are properly 
chosen, the width of the fringes seen through them 
will be most favorable for measuring the fractional 
parts of fringes. Also, the position of the zero- 
order fringe, relative to point E, may be calibrated 
to read directly the angle between the ends of the 
gage blocks. 

The recommended procedure for measuring a block 
is: (1) Adjust the two sets of fringes seen on the ends 
of the block so that they are perpendicular to the 
dividing plane, ‘as in figure 4F; (2) note the position 
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hig. 6C); (3) rotate the block 180° about a vertical | that is perforated so as to transmit the required part not 








axis through its center; (4) readjust the fringes to | of the light beams used for making mez al ewes 
restore the condition of (1) above; (5) again note the | See also figure 8 for a vertical view of M, al its| fron 
position of the zero-order fringe (or absolute order at | supporting parts. The plate, M,, rotates in an anny vy, 
E) ; (6) the difference in the two observed orders at E, | lus ring, Mz, which in turn is rotatable, from outsid 
or positions of the zero-order fringe, multiplied by the | the tastrumeat, in another annulus, Ms. The aper- 
constant K, described above, is a measure of the | ture in M, has its center displaced from the center of 
angle between the ends of the block. the instrument in a direction parallel to the dividing 

Figure 7 shows two vertical sections through the | plane of P and by an amount equal to one- half the } 
center of the short-block interferometer that are | horizontal separation between the centers of they, { and 


ali 
( ple un 


mutually perpendicular to and through the centers of | gage block, G,, and its image, G’. The gage block} « = 
each other. A pinhole, S, illuminated either with | rests on the center of M, (see fig. 6B). Accurate lh 
monochromatic or polychromatic light from outside | Placement of the block is facilitated by stops. ser 
the tube H, serves as source. The position of the | The eccentric annulus, M : has an arm that pro-| ;, , 
pinhole is adjustable in the focal plane of the collima- | jects through the wall, Q, of the instrument and is 
tor lens, L, and the collimator tube is adjustable in supported by this arm at one of its three supporting ie 
length. The prism, P, rests on a thick plate, I, to | pomts by an adjusting screw, R,. The other two i 
which is fastened the lens-holding plates, J and J’. | supports for M; are steel balls, T; and T,, (fig. 8B),| | 
The plates, K and K’, which cover the ends of P, | Which are held in conical holes by means of UW! has 
are not fastened to J or J’. Consequently, small | screws, U, and U,, respectively. The ends of these oh 
stresses applied to tube H’ while adjusting the evye- | Screws have eccentric conic ‘al depressions that permil f 
piece or manipulating a micrometer in it, are not | ® limited amount of rotation of NI about an 8) 1, 
ened to the dividing plane of P. The screw Ry)" 


a ; > 
appa . - :; : permits fine adjustment of M., and conseque ntly the 
_ The material of the instrument, except for the col- gage block, which it supports, about a horizontal (ep 
limator, the eye-piece tubes and optical elements, is | jyjs parallel to the dividing plane of P. A similar} opt 
made of steel. Steel was chosen because its expal- | pair of balls, screws, and the adjustab le screw Ry} Aft 
Sivity approximates that of the glass elements more permits rotation of P, about two axes parallel 10} alo 
ne arly than other usable materials those used for adjusting Msg. wit 
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FIGURE & Several views showing individual parts of the short-block 


ynterferometler. 


If the refracting edge of the optical wedge, W,, is 

made perpendicular to the dividing plane of P (i. e., 

‘parallel to the plane of fig. 7B), each pair of compo- 

nent rays will traverse this plate at points of equal 

plate, M thickness. ‘onsequently , W, in this orientation, does 
ulred parts not affect the fringes of interference. It does, how- 
surements} over, serve as a window, protecting the prism P, 
Mi and it} fom the accumulation of dust. Other functions of 


1 an anny W, will be discussed later. 
om. outsid 
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The prism, P, is centered over the aperture in plate 
I (figs. 7A and SA) with its ends parallel to plates K 
and K’ (fig. 7B). Using the base surface of P as a 
(plane mirror, the pinhole aperture 5 is located in the 
focal plane of lens L. by varying the length of tube H 
and at a point in this plane where the light will form 
m image of the pinhole upon itself. This places the 
light beams, after division at the dividing plane, 
normal to the prism base and, consequently, parallel 
toeach other. 





A gage block, whose end faces are parallel to each 
other two “ther, is placed on the center of plate M, with its 
(fig, 8B),| wer surface parallel to the top surface of M,. When 
, the surface of the block is adjusted parallel to the 
base of P by means of screws Ry, U,, and U, the 
at permit light from S is reflected normally from the top surface 
t an axis of the block and observed at the eveplece aS inter- 
screw h ference fringes 

iently the' The parts of the light beams that are not inter- 
iorizontal| epted by the block and its support M,, traverse the 
A similar) optical wedge W,, and enter the right-angle prism Px. 
screw Ry} After two internal reflections in Px, the light returns 
arallel ad ilong a path that is symmetrical to its incident path 
‘with respect to the 90 The 90° edge 


is of two| 
s of these! 


edge of Pe. 


of Py is made normal to the dividing plane of P 
by means of screws in its support similar to U,; and 
U, and the fine-adjusting screw R, (fig. 7A). This 
light forms the background fringes used in the test, 
but no interference is observed until this prism edge 
is nearly normal to the dividing plane. When it is 
not normal the two images of this edge, formed by 
the two component light beams, intersect in the 
extension of the beam dividing plane. The prism 
P, is rotated until the two images coincide—perfec- 
tion being attained when the background fringes 
have maximum contrast. 

The above-mentioned contrast in the background 
fringes is not affected by screw, Ro, ‘since it has no 
vertical rotational effect on the prism. Its effect 
is to change the width of the background fringes only. 
When the background fringes are made infinitely 
broad, the collimated beam of light returns toward P 
parallel to the incident beam—all rays having 
suffered a horizontal shift in Pp, as illustrated in 
figure 4A. 

The holes in the gage-block support, M, (fig. 8C), 
are so-spaced that when M, is rotated to one of the 
four positions for which the rectangular sides of the 
block are either parallel or perpendicular to the 
dividing plane, all light that goes through M, 
returns again through it. That is, the apertures in 
M, are symmetrical both with respect to the dividing 
plane and to the 90° edge of Pg. Parts of the beam 
(rays 3, 3’, 4, and 4’ in fig. 4) will pass downward 
through M, to Px, shift horizontally in Pg, and pass 
upward through other apertures in M, to P. Other 
parts of the beam (rays 2 and 2’ in fig. 4) will pass 
downward threugh M,, shift horizontally in Pz, pass 
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upward through M, to the lower surface of the gage 
block, return through M, to Py and again upward 
through M, to P. If the two end faces, G, and Gag, 
are parallel, the light will be incident on G, and G, at 
equal angles. When G, is adjusted normal to the 
light, G, will also be normal to it. 

We have three sets of fringes, shown in figure 4F’, 
to consider: (1) The set F, is formed by light re- 
flected from the top surface, Gi), 


») 


(2 


of the race block: 
) the background fringe between and about F, 
and F;, formed by light reflected from P, but not 
incident on the gage block: and (3) the set F. 
formed by light reflected from Py to G, and back 
through Py. The direction orientation of F, 
determines the angle between G, and the incident 
wavefront; the width of the background fringes 
determine the direction between the incident and 
reflected beams to and from Pr; and the orientation 
of F, determines the angle between G. and the wave- 
front that is reflected from it. 


498 







The set of fringes, F, is adjusted by means of R 
normal to the dividing plane of P, for which eon 
tion G, is normal to the incident light. The back 
ground fringe is made infinitely broad, for whl! 
condition the light beams returning from ] “ 
parallel to the incident beams. If the 
of G are parallel, the set of fringes a 
G, will be parallel to those on G, an indicated , 
the set F, in figure 4F, If G, and Gi. are not paralld 
the fringes seen on them and indicated as F, and P| 
in figure 4F’ will not be parallel to each Other | 
The angle between these two sets of frin 
measure of the angle 
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Correction for Instrumental Drift in Flame Photometry 
Barry W. Mulligan and Alan F. Haught 


The 


mental measurements, specifically in flame photometry, 


application of the principles of experimental design to correct for drift in instru- 


is described. Based on the procedure 


proposed by Youden, samples and reference standards may be measured in combinations of 


pairs in 
of drift can be determined. 
measured 
the procedure 
severe drift, 


i range 


to flame 


of 0.6 to 1.2 


3.5 TO : 


1. Introduction 


In any series of measurements, drift in the response 
of the measuring weg og may introduce appre- 
ciable errors Alth iough the errors ms L\ be reduced 
by improved adjustment or design of the instrument, 
some drift may be inevitable. 

An unusual example of drift was encountered 
measurements with a flame photometer. Replace- 
ment of the detector, a photomultiplier tube, was 
later found to overcome most of the drift; however, 
the possibility of compensating for drift by appro- 
priate design of the experiment in this and similar 
cases Was explored 

Youden' proposed a design and statistical treat- 
ment of data using a method of mcomplete blocks to 
minimize errors caused by mstrumental drift, and 
at the same time, to measure the magnitude of this 
effect. The applicability of his procedure to this 
drift problem was examined. 


2. Experimental Procedures 


During the development of a method for deter- 
mining the alkali metals with a Beckman model 
DU* flame photometer with photomultiplier attach- 
ment, it was noted that the readings of the instru- 
ment were shifting with time. The determinations 
were bemg made on a long series of solutions, so that 
a shift in readings would seriously affect the results. 

In Youden’s method to compensate for drift,’ the 
objects to be measured (here the sample and standard 
solutions) are arranged in an appropriate sequence. 
Specifically, the objects are grouped in pairs, so that 
each object appears 1a a pe uir with every other object. 
For example, 5 objects are grouped in 10 pairs, as 
shown in table 1, the sequence of the pairs being 
immaterial. The objects are then measured on the 
instrument in a manner such that the time between 
measurements on the objects in a pair is small in 
comparison to the time between pairs. 

The rigorous method Youden for computing 
the drift involves long arithmetic calculations. A 
simplified, yet accurate, alternate method, suggested 
later by Youden, will be described here. No dif- 
ference was found in the values computed by the 


of 
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i manner such that errors caused by drift can be 
Where the drift is dependent on the magnitude of the quantity 
grouping by magnitude prior to final measurement is necessary. 
photometric determination of alkali elements, 
resulted in improving coefficients of variation from an original range of 2.0 to 





minimized and the actual amount 
Application of 
under conditions of 


TABLE 1. Design of experiment 


Pair__. 
Object AB 
Pair 6 
Object 


CD 
7 s 9 10 
DA 


within the accuracy of the data. The 
method involves, first, the determination of cal- 
culated, corrected, values for the measurements 
on the objects, B, C, ete., and second, the com- 
putation of drift using the observed values and cal- 
culated values of the objects. 

calculated value for each object is determined 
as is shown by the steps in table 2. In the table, 
A, represents the measured value of object A in 
pair 1, C; the measured value of object C in pair 7, 
etc. The differences between A and the object 
with which it appears are summed, +A and —A 
are added to the left side of the equation, the equa- 
tion is divided by five, and everything but A is 
transposed to the right side. Since the expression 

B,+C;+ D+ E,)/5 will be approximately equal 
to the average of the 20 measurements, (that is, 
four each on A, B, C, D, and E, divided by 20), the 
last equation is obtained. A similar computation 
is performed to obtain the calculated values for the 
other objects, B, C, D, and E. 

In addition to correction for drift, the method 
also provides a measure of the drift. Table 3 shows 
the method of computing the average drift for each 
pair. The differences between the observed values 
and the calculated values of the two objects in the 
pair are averaged, and the result represents the 
average drift for the pair. The drift for the pairs 
is then plotted as the deviation from the mean drift. 
It should be noted that the validity of the statistical 
method rests on the assumptions that the objects do 
not change in the course of the measurements and 
that the drift is independent of the magnitude of the 
measurements made. 


two methods, 


TABLE yA Calculated value for obje ct A 


A; —Bi=X 
Ar—C xX 
A,o—Dye=X 
Ay— Ey=Xy4 


{1A —(Bi+Cr+ Diet Ex 


A = Average of all measurements +— 
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TABLE 3 Computation of drift in pair 1 


Value Value 
observed calculated 


Difference 


rage drift= 


Figure 1 shows an instrumental-drift curve ob- 
tained by this method for a series of readings on 
sodium. This is a normal type of curve with narrow 
limits. The abscissa is the number rather than the 
time of the reading, since the only time requirement 
is that the 2 objects in any 1 pair are immediately 
consecutive. Thus, for example, the third point 
on the graph is plotted midway between readings 
5 and 6, since these 2 readings were used to compute 
the point. The ordinate is the deviation from mean 
drift in terms of instrument readings and in this 
particular case one instrument unit represents ap- 
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pre sent 


ite unit equa pproximately 2 percent ¢ f sodiun 


The value of this method became immediately 
evident upon examination the curve shown in 
figure 2. This represents a series of readings on 
potassium in which, unexpectedly, the range of the 
instrument had be changed from low high 
during the measurement of the series. The instru- 
ment response was evidently slow 

The original method emploved in these deter- 
minations consisted of making measurements on a 
series of unknown solutions with standard solutions 
measured periodically during the series. A large 
change in readings on an object during any series 
of readings with other objects is, of course, evident 
to the observer. Instead of having to discard data 
and begin over, as was done in the original method, 
the information obtained when using a drift correc- 
tion is perfectly valid. Similarly, any small in- 
strumental drift, unnoticed by the observer, is 
corrected for in the calculations. Another advantage 
of the method is that the experimenter may in- 
terrupt a readings at any point between 
pairs, and upon return may resume at the point of 
interruption, rather than start over again 

The coefficients of variation obtained when using 
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Dh nge of the i s 1S and 19 
this refinement show a marked improvement over 
those obtained with the previously mentioned origina] 
method, in which no instrumental drift correction 
was applied. Coefficients (obtained from 10 deter- 
minations on each of 2 samples of lithium solutions 
and 2 of potassium) which were in the range of 2.0 
to 3.5 were improved to fall in a range of 0.6 to 12 
by the new method. 

In an effort to determine the nature and cause of 
the drift, the statistical method of incomplete blocks 
was again employed. A sodium determination was 
made on 10 solutions of widely different concentra- 
tions under standard-operating conditions in a sta- 
tistically designed sequence during a 2% hr time 
interval Krom the observed readings, the drift was 
computed and was plotted as a function of time. 
The curve obtained (fig. 3) indicated an initially 
rapid but continually decreasing drift with an 
irregular drift’ superimposed on it. Because the 
“reality” of the superimposed drift was doubted, 
additional measurements were made on a solution 
containing 3 ppm of lithium and on a solution con- 
taining 20 ppm of potassium. Experiment had 
shown that determinations of these elements are less 
subject to random fluctuations than those for sodium 
and they are less liable to error because of contami- 
nation by contact with dust in the air. The runs 
in each case were made on a single solution and the 
drift was computed as the deviation of the observed 
readings from their mean. With both the potassium 
and the lithium solutions only the initially rapid but 
continually decreasing drift was obtained, and none 
of the superimposed, irregular drift was present. 
This would indicate that the drift may be dependent 
on the magnitude of the measurement being made, 
contrary to one of the assumptions of the statistical 
method. 
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FIGURE 3 Instrument drift obtained when measuring a series 


of 10 sodium solutions having widely different concentrations 

To test the effect of a wide range of concentrations 
on the drift, runs were made on 2 lithium solutions, 
| twice the concentration of the other. The drift in 
each case was computed as the difference between 
the observed readings and their mean and was 
plotted as a function of time as shown in figure 4. 
The ordinates of the two curves are in a constant 
ratio at all points along the abscissa, within the 
accuracy of the readings, showing the drift is de- 
pendent upon the intensity of the light being 
measured. 

The variation of the drift with the incident light 
intensity explains the irregular character of the 
drift obtained previously (fig. 3) on the group of 
sodium solutions having different sodium concen- 
trations. The larger deviations were caused by the 
large drift of the readings for the highly concentrated 
solutions and the points closer to the average drift 
eurve by the small drift for the less concentrated 
solutions. The particular sequence of solutions of 
high and low concentrations resulted in the irregular 
drift. 

The cause of the drift was finally ascertained to 
be due to fatigue of the photomultiplier tube, a 
phenomenon usually associated with these tubes 
but only at much higher phototube current levels 
than those encountered here. This effect may be in 
either a positive or negative direction, although it is 
most commonly negative. 

The phenomenon of fatigue is not a strictly re- 
versible process, that is, the recovery curve is much 
less steep than the fatigue curve. However, the 
statistical method described here is based in part on 
the assumption that the drift is independent of the 
measurements made; accordingly, a modification 
was made in the procedure to ensure accurate results. 
When a particular element in a group of unknown 
solutions is to be determined, a quick initial run is 
made without the use of any statistical sequence or 
repetition of readings This gives a good approxi- 
mation of the concentrations of all the unknown 
solutions. The solutions are then grouped accord- 
ing to these preliminary results so that the concen- 
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lithium solutions. 


Instrument sensitivity is approximately 0.04 ppm per unit. 


trations of all the solutions within one group are 
very near one another. Standard solutions for each 
group are selected to approximate the concentrations 
of the group. The group of samples and standard 
solutions is then measured using the method of 
incomplete blocks, the normal drift curves are com- 
puted, and the results are corrected for drift. There 
is no irregular trend of the curves, since all the 
readings in a group are close together. This group- 
ing is also advantageous from the standpoint that 
smaller errors in interpolation are introduced when 
readings on standards and unknowns are close 
together. 


3. Conclusion 


This paper has discussed a problem of instru- 
mental drift in the laboratory to which a statistical 
solution involving incomplete blocks has been 
applied. While the magnitude of the drift en- 
countered was not characteristic of the instrument, 
it served as a good example of a complicated drift 
problem. The statistical method improved the 
precision and provided a measure of the drift, except 
when quantities of widely different magnitudes were 
intercompared. Investigation showed that the drift 
resulted from fatigue of the photoelectric tube 
involved in the measurement, and that the drift was 
dependent upon the intensity of the light to which 
the tube was exposed. When drift is dependent on 
the magnitude of the measurement, it appears 
necessary to group the quantities according to order 
of magnitude before final measurements are made 
and the drift correction is applied. 


Suggestions received from W. J. Youden during 
the course of this work are gratefully acknowledged. 
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; Mechanism of Stress-Corrosion Cracking in the AZ31B 
Magnesium Alloy 
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The mechanism of stress-corrosion cracking of annealed AZ31B magnesium alloy in an 
aqueous NaCl keCrQO, solution was investigated. Cracking was predominantly an electro- 
chemical process and was shown to be dependent on the rate of strain in the specimen follow- 

f ing loading in tension. It is postulated that cracks develop if the protective film on the 


metal surface is ruptured 
repaired in the corroding medium 
filmed metal 


1. Introduction 

) An investigation of the mechanism of stress- 
corrosion ¢r acking of metals, sponsored jointly by 
hhe Corrosion Research Council and the National 
Bureau of Standards, is in progress. Certain of the 
magnesium alloys are admirably suited for such an 
While reports of failure in service 
laboratory techniques that will 
soduce stress-corrosion cracking in minutes are 
vel known. As a result, extensive series of experi- 
ments can be performed in a conveniently short time. 
E frsmed work on the AZ31B magnesium alloy in 
s laboratory [1] ' and elsewhere [2, 3], it was 
es that in certain corrodents a threshold stress 
ould be obtained below which cracking would not 
geeur in 1,000 hr. Above this stress level, specimens 
lwould generally fail in less than 5 min. Furthermore, 
\. was found, in this laboratory and elsewhere, 
that the threshold stress Was procedure sensitive [4]. 
It was appreciably higher if the specimens were 
loaded in tension before the corrodent was added 
than if the specimens were already immersed in the 
erodent when the stress was et While a 
threshold stress was readily obtained in laboratory 
lnxperiments, such threshold stress ws not been 
lobserved in weather-exposure tests, most probably 
ecause the ‘se specimens were also subject to gener: al 
Furthermore, there was little difference 
naverage times to failure, at a given stress, whether 
he specimens were exposed in a marine atmosphere 
it Hampton Roads, Va., or in an urban atmosphere 
the Bureau. It was noted that, in both environ- 
most failures during or following 


investigation. 
are extremely rare, 





both 


wrrosion. 


/ 


iments, occurred 


(5, 6] showing the 
cracking of magnesium alloys have 
Both showed plastic deformation 
preceding fractures, and Harwood [7] has suggeste “ 
plastic extension must precede 
wrosion cracking in all materials 
Priest [8] reported that cracking of the AZ61 alloy, 
ee started, could be stopped by the application of 
ithodie protection. It was considered important 
determine whether or not this was also 
ith the AZ31B alloy and if so, to see what Giaaee- 


Two excellent 
stress-COrrosion 


eel) prepared. 


motion pict ures [5, 


dat some stress- 


possible 





—_— ndicate t terature references at the ¢ this paper 


over narrow segments of the specimen) 
This would expose film-free metal that is anodic to the 


at a greater rate than it is 


tion the application of cathodic protection might 
give as to the mechanism of stress-corrosion cracking. 
The results of these investigations are reported in 
this paper. 


2. Materials and Experimental Procedure 


Four lots of the annealed AZ31B alloy material, in 
two were used in the investigation. The 
chemical compositions of these materials, as reported 
by the manufacturers, did not differ appreciably 
from the nominal composition of aluminum 3.0 
percent, zinc 1.0 percent, manganese 0.3 percent, 
with the impurities iron, nickel, and copper be low 
the maxima permitted in ASTM Specification B90 
57T. Both 0.064- and 0.125-in.-gage materials were 
used; however, most of the work was done with the 
thicker The materials were machined into 
,-in. reduced-section tensile specimens. In order to 
remove any impurities that may have been present 
on the surface of the sheet as a result of rolling, most 
specimens were immersed in an aqueous solution 
containing 30 ml of HNO, and 10 ml of H,SO, per 
liter until approximately 0.001 in. of material had 
been removed from each face. A few specimens 
machined from sheet that had been chrome pickled 
were used with the original surface, intact. The 
specimens were mounted into cells similar to those 
used in the earlier work [1] and were stressed in 
tension using lever systems. The corrodent [9] 
was an aqueous solution containing 3.5 percent of 
NaCl plus 2.0 percent of KsCrO, used at 30° - +2° C. 

The threshold stress was determined at 30° C for 
the 0.125-in.-thick materials, in terms of the percent 
of vield strength of the individual lots of material. 
These yield strengths (0.2% offset) ranged from 
21,900 to 23,100 Ib/in?. 

in stressing the specimens, they were first subjected 
to the load of the lever alone, the equivalent stress 
being 2,000 to 4,000 Ib/in.?, depending on the specimen 
thickness. The full load was then added as a unit 
without impact in 1 to 3 sec, thus increasing the stress 
to the desired value? of 21,500 to 23,000 lb/in.2 Wire 
strain gages attached to opposite faces of specimens, 
read individually, were used to determine the strain 
produced in specimens by the addition of both the 


gages, 


gage. 


? The lever system used ir stressing the specimens had aratioof25tol. It was 
easily possible for one person to apply the full load to the lever without overload 
or impact. 
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minor and major stresses. As was to be expected, 
the major strain accompanied the application of the 
major load. It was noted, however, that specimens 
exposed in air at 30° C continued to extend after 
loading as long as the stress was applied. For the 
first few minutes, at any rate, the logarithm of the 
strain rate decreased linearly with the logarithm of 
the time. 

Wire strain gages, of course, could not be used to 
determine the strain rate of specimens immersed in 
the corrodent. It was found, however, that in the 
strain-rate range that was to prove most important 
in the investigation (namely less than 500 10o-° 
min) a relationship could be obtained between the 
strain rate determined using wire strain gages and 
that computed from dial-gage readings, the dial gage 
being actuated by the lever system. ‘Thus, it was 
possible to follow the strain rate of a specimen stressed 
either in the corrodent or in air (beginning approxi- 
mately 30 sec after stressing) from periodic readings 
of a dial It was also possible to follow the 
strain rate of a specimen to some preselected value, 
using wire-strain and then rapidly strip off 
the gages and add the corrodent. The extension at 
the moment the corrodent was added could be 
estimated very closely from the extrapolated strain 
rate-time curve. 


gage. 


r ron 
yages, 


3. Results and Discussion 


3.1. Specimens Immersed in Corrodent Before 
Application of Stress 

The strain rate increased markedly as the stress 
was increased above the threshold value (approxi- 
mately 954° of the vield strength). After tabu- 
lating the data from many specimens it was evident 
that if the corrodent was added before the major 
stress was applied and the strain rate exceeded 500 

-10010-*/min 1 min after the stress was applied, 
early failure by stress-corrosion cracking was to be 
expected. If the strain rate had fallen below this 
value, specimens generally did not fail in 1,000 
min and were then removed from the corrodent; 
the earlier work had shown that if a specimen did not 
fail in this period it probably would not fail in 10,000 
min. The times for failure for the various specimens 
are plotted against the percents of their respective 
vield strengths in figure 1. 

In an earlier investigation [10] i 
straining of a magnesium alloy ptt its 
strength ruptured the protective surface film and 
exposed film-free metal that was approximately 240 
mv anodic to the filmed surface. If the stress was 
held constant so that the strain rate fell off rapidly, 
the protective film re-formed on the surface where 

had been ruptured. 

It is postulated that 


was shown that 
vield 


threshold stress 


above the 


and at very high strain rates the protective film is 
ruptured on most of the surface of the specimen 
and corrosive attack is of a general nature. As the 


rate decreases, certain crystals work harden rapidly, 
the strain rate on them decreases to the point W here, 
with the aid of the K.CrO, inhibitor, a stable film 
is formed \.\ith the average strain rate above 
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500 10-°/min, however, the strain rate will bel 


sufficiently high on other crystals so that the rup-|* 
tured film over narrow areas perpendicular to the 
direction of strain will not be repaired. They 
film-free areas will be anodic to the filmed areas jy 
the NaCl-K.CrO, solution and = most probably | 
will be very small compared to the cathodic areas 
Electrochemical corrosion will result and unles 
the strain rate decreases to such an extent that the| 
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protective film can re-form, stress cracks will he 
initiated and rapid failure of the specimen may ya 
follow. ase. 
If no cracks are formed before the average strain) Psy 
rate falls below 500 10-°/min it is. believed that expel 
film repair generally is as rapid as film rupture 99) i)” p 
that stress-corrosion crac king does not occur. pe 
3.2. Specimens Stressed Before Immersion in the a 
Corrodent a 


that the specimen {fH,SC 


It would be expected, then, 
in air and the|the ¢ 


would not fail if it had been loaded 
corrodent was not added until the strain rate had {mitt 
fallen below 50010-*/min. Data indicated, how-| subse 
ever, that if the specimen had been given the nittie sult 
acid-sulfuric acid pickle prior to it would fail , of th 
unless the strain rate had decreased to less than Th 
100 10-*/min before the addition of the corrodent. /fom 
A typical strain rate-time curve for one of the spect {at a 
mens is shown in figure 2. It is seen that the¢subse 
strain rate increases very rapidly with the acldition | reas 
of the corrodent. The maximum strain rate reached | ens 
was of course the average over the whole specimen strat 
The subsequent behavior of the specimen indicated {at t 
that the measured strain was most probably con- mites 


use, 


centrated in short segments. ir ollowing this rapid | strail 
increase in the strain rate, it decreased over most )itom 
of the specimen length. However, the extension} lpr 


rate was still so high over several short lengths of the} Fr 
specimen that the protective film was not able to} made 
re-form and stress-corrosion cracks deve loped, prob | ipo: 


ably between the 6th and 7th min. They had )sigg 
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probably 
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nd unless 
t that the 
cS wil] by 
men ma) 


‘stress been removed 


penetrated so deeply by the 8th min that the average | 
train rate for the specimen was increasing rapidly 
and early failure would have occurred had not the 
Subsequent examination of 
the specimen showed that the cracks were segregated 
in several short lengths of the specimen and not 
more or less uniformly distributed as was the usual 


ase. An anomalous change in strain rate has been 


ieported for polycrystalline iron [11] and for single 
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crystals of cadmium [12] under somewhat similar 
experimental conditions. No report of such behavior 
in magnesium alloys has come the author’s 
attention. If explanations given for other materials 


to 


‘are also applicable to this alloy, it is suggested that 


| the 


from the metal was 
blocked by the film formed by pickling in the 
HSO,-HNO, solution. Contact of specimens with 
the corrodent momentarily dissolved this film per- 
mitting the escape of trapped dislocations and the 
subsequent slip of individual crystals with the re- 
silting marked change noted in the rate of extension 
f the specimen 

The corrodent was added to specimens, machined 
fom a chrome pickled sheet, that were extending 
iia strain rate as high 280><10~-°/min without 
subsequent failure. There was no measurable in- 
rease in the average strain rate with these speci- 
mens when the corrodent was added. At higher 
train rates there was apparent dislocation escape 
it the machined with high local strain 
utes but not appreciable over-all increase in the 
‘train rate. Failures resulted in these instances 
om cracks penetrating into the metal from the 
uprotected edges 

From these data, a tentative explanation can be 
made of the fact that most failures of weather- 
posure specimens occurred following rain. It 


escape of dislocations 


as 


edges 


Is 


They had\siggested that solutions formed by water collecting 


on the specimens would dissolve the oxide film over 
some areas on the metal sufficiently to permit the 
escape of block dislocations. This would allow slip 
on certain microscopic areas and microscopic exten- 
sion at a sufficiently high rate to initiate stress- 
corrosion cracking in the rainwater solution. Prop- 
agation of the cracks was sometimes stopped by 
evaporation of the solution before fracture was com- 
plete, as is shown in figure 3 taken from the earlier 





NBS. 


stress-corrosion crack extending about \4 through specimen. 


FIGURE 3 Specimen exposed to weather at Note 


Cracks of this type, found after rain, are believed to have been initiated by 
microscopic strain after dissolution of the oxide film by the rainwater-corrosion 
product solution which permitted the escape of blocked dislocations. Propaga- 
tion of the cracks continned in this solution until failure occurred or evaporation 
of the solution stopped the reaction 
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paper [1]. In many cases, however, specimens were | the time just before and immediately follow; 
. . 0 





found to have failed completely after a rain. removal of the cathodic protection. It is seep th 
the specimen began to extend rapidly within > iin 
3.3. Effects of Cathodic Protection after the protection was removed. It would hak 
failed within another minute had not the stam 
The effects of the application of cathodic protec- | been removed. There was no lasting polarizatigg | 


tion to a specimen in which stress-corrosion cracks | of the specimen as the result of cathodic Protection 
are developing offers much information as to the | With its removal, the potential of the specimes | 
mechanism of the corrosive process. It was deter- | immediately became approximately 400 my mop 
mined that 0.125-in.-thick specimens stressed to | positive than it had been before protection wy 
98.6+0.4 percent of their yield strength would fail | applied. The potential soon changed to the nog. 





























in an average period of 152 sec (standard deviation | inal value of 1.5 v and as the extension of th 
39 sec). Cathodic protection * was applied to other | specimen opened up  stress-corrosion cracks and ( 
specimens subjected to the same stress after expo- | exposed unfilmed metal, became still more negatiyp 
sure periods of 90 and 105 sec. A typical strain 
rate-time curve for these specimens is shown in fig- . 
ure 4. Examinations, using a binocular microscope, | 4.48x10 , . 
of specimens that had been removed from the cor- 
' 1 fter 35 10) ; ' QR 5 ——o—— EXTENSION 
rot —_ © tel 5 T5) to 40 see exposure (stressed to YS.o ——e—— POTENTIAL 
to 99% of the vield strength) indicated that crack- 
ing began during this interval. Hence, it is logical a rr 4 dees 
. . . . t 
to assume that stress-corrosion cracking is an impor- _ 
. . ° . c 
tant factor in the increase in the extension rate as = 1 | 
shown in figure 4 after 1 min of exposure. The 3 | 3 
. . . . . (7) 
marked change in extension rate with the applica- 2 _ | 
tion of cathodic protection (at a current density of r 440 F l ‘ 4—i5' Ficu 
0.5 ma/em* over the exposed area) indicated that ” 7 PROTECT r 
" ‘6 a CA I 
cracking was almost if not completely stopped. a Ss RS ong 
—— ' . 5 REMOVED 
This is substantiated by the fact that a specimen re : l Le 
subjected to the conditions described above had | =| a | | 
not failed in 20 hr although it had extended about F “UT in. fe 
0.04 in. Figure 5 shows the total extension aad rood lv. i. 
potential curves, for this specimen, plotted against a mens 
A 
lindrical anode of sg | ow 
2A cylindrical anode of stainless stecl was placed around the inside of the , . 
cell. The ends of the specimen were masked off with *‘ Duco Cement” so that 4.32 —— L__ ! —L failu 
only 15 cm? of the reduced area of the specimen was exposed to the corrodent 1260 1265 1270 1275 1280 resul 
7 TOTAL EXPOSURE PERIOD, min the : 
tenia emia atintipsittiadinmivareats oe ; 7 ' 0 
4000 x10 TT TiTrTtt T FiGuRE 5 Total extension-lime and potential-time curves for + 
- ] specimen after immersion for 1,260 min in the corrodent ihe 
O - CATHODIC poin 
i“ ) PROTECTION Note the change in potential with the removal of cathodic protection. Nomina = 
foe <A APPLIED potential of material with respect to a saturated calomel electrode was -Li irom 
i", \ Che increase in extension rate after the removal of cathodic protection ind: abov 
1000 ; cated that early failure of the specimen would have occurred had the stress not ; 
= ] been removed 
Aad b 4 
5 ’ { 3.4. 
> * J 
= _ ; 
= _ 1 The fact that stress-corrosion cracking can be 
w - ] stopped by the application of cathodic protection If 
ps L and that it will start again after the removal of this the 
= . . ° . . ° °, 8 1¢ ¢ 
= protection indicates that, in this alloy at least, it is| and 
Z wot primarily electrochemical. Specimens removed from) poy 
? = the solution after the application of cathodic protec-| the 
ee —_ . . . . . . n { 
o \ tion, contained cracks as is shown in figure 6. The pene 
x depths to which cracks penetrated into the metall ,. | 
‘“ . . . . 
s : were determined on longitudinal metallographic see p46 
a | tions taken well away from the edges. Maximum] 4) 94 
crack depths were determined in 0.125-in.-thiek} 5 9 
- Lope ra ee aoe specime ns to whic h protection was applic d afte US a 
01 LO 10 so | and in 0.064-in. specimens to which the protection was| 5) Ds 
TIME, min added after 60 to 65sec. Cracking was assumed to 
SF «a » . . In 
Ficure 4. Extension rate-time curve for specimen stressed to | have started in each case after 35 sec exposure. Pen- 
approximately 99 percent of its yield strength etration rates for some of the deeper cracks amounted as h 
+: : : : ; ; ° as [1 
. LT a 5 rs , ’ , ’ ods 
Deviation of the curve from i straight line (after 1 min) was duc to development to 4 - 10 In. /Sec. In the ave rage exposure pe n othe 
of stress-corrosion cracks First cracks develop in 0.6 to 0.67 min The applica- | to failure for the two thicknesses of specimens, the 
tion of cathodic protection stops the crack propagatior fter 3 min the extension | : ; ‘ 7 
rate 2 " sqgneainntety that s he cupeeted had Deanna “ racer ; penet rat 1ons would be approximately ).02 and 0.04 } 
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Edge of specimen given cathodic protection SO sec 


FIGURE 6. 4 
afte) essa of major load. 


Cracks extending into specimen from edges are believed to have developed 
in period between 35th and S0th sec afte r loading. Specimen was removed from 
corrodent after approximately 18 hr exposure x Ww 


in. for the 0.064- and 0.125- in. specimens, respective- 
ly,i. e., 30to 37 percent of the thickness of the speci- 
mens. Satisfactory measurements of the depth of 
penetration of the cracks at the point of 
failure are difficult to make, in that failures frequently 
result from a crack or cracks that penetrated into 
the specimen simultaneously from an edge and the 
opposite faces. It is believed that the weakening of 
the specimen by stress-corrosion cracks, to the 
point where it will fail by tensile fracture, resulted 
from penetration of the cracks at the rate indicated 
above, and was therefore primarily electrochemical. 


stress 


3.4. Current Density Associated With Stress-Cor- 
rosion Cracking 


If it can be assumed that Faraday’s law applies to 
the dissolution of a metal as well as to its deposition, 
and if the electrochemical equivalent of the alloy is 
then taken to be the same as that for magnesium, 
the current density necessary to produce a crack 
penetration of 4.0 10~* in./see may be easily shown 
to be approximately 14 amp/em*®. It should be 
noted that this calculation assumes only that the 
anodic reaction is 100 percent efficient; it requires 
no assumptions as to the width of the anodic area at 
the tip of the crack. If the anodic reaction is even 
in part, 


Meg Meg ré 


’ 


as has been suggested by Hoey and Cohen [13] and 
others instead of the conve ntional 





_ slower rate, in 


Mg=Meg** 


19 
g 2é, 


the current are materially 


than 


density requirements 


reduced. They are, however, much higher 
any that have come to the author’s attention for 
other materials. Hoar and Hines [14] estimated 


that a current density of 1.5 amp/em?* would account 
for the rate of crack penetration (2.2 10~° in./sec) 


that they found in austenitic stainless steel in a 
boiling 42 percent MgCl, solution. They further 
report that under certain conditions, iron can be 


anodically dissolved at five times the rate of propa- 
gation of the cracks in the austenitic stainless steel 
reported above. 

It is the author’s experience that the 
corrosion cracking of austenitic stainless steels in a 
boiling 42 percent MgCl, solution is a very much 
slower process than that of the AZ31B alloy in the 
NaCl-K.CrO, solution. No failures of 0.064-in.- 
thick stock of type 304 stainless steel have occurred 
in this laboratory in less than 1 hr, although cracking 
is believed to have been initiated in less than 30 min. 

The NaCl-K,CrO, corrodent had a specific resis- 
tance of 11 ohm cm; with a difference in potential 
between the cathodes and anodes of less than \ v, 
it is seen that the cathode areas must be largely 
within the cracks if current of the magnitude 
indicated above is to flow. The motion pictures 
mentioned above showed copious quantities of gas 
(probably hydrogen) escaping from the stress- 
corrosion cracks which would also indicate that the 
cathodes were located primarily on the sides of the 
cracks. 


stress- 


4. Summary 


1. A study has been made of the mechanism of 
stress-corrosion cracking of the AZ31B magnesium 
alloy in an aqueous solution containing 3.5 percent 
of NaCl plus 2.0 percent of KyCrQ,. 

2. Specimens dead loaded in tension by means of 
lever systems continued to extend for an appreciable 
time following stressing. Specimens stressed in the 
corrodent and extending at an average strain rate in 
excess of 500+10010-°/min, 1 min after the stress 
was applied, generally failed in less than 5 min by 
stress-corrosion cracking. Specimens extending at a 
most instances, did not fail in 1,000 
min. 

The strain rate for specimens stressed above 
the threshold stress prior to immersion in the corro- 
dent must fall appreciably below 500X10~°/min or 
the specimens will fail. The strain rate on addition 
of the corrodent depended on the surface condition 
of the specimen. It was suggested that the addition 
of the corrodent to specimens, pickled in the 
H.SO,-HNO, solution, permitted the escape of 
blocked dislocations and an appreciable increase in 
the average strain rate. Escape of dislocations from 
abraded or chrome pickled surfaces was prevented 
(except in some highly localized areas) and no in- 
crease in the average strain rate was noted. 
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The advancement of stress-corrosion cracks was 


halted for at least 20 hr by the application of 
cathodic current of 0.5 ma/em?* of exposed area. 
With the removal of cathodic protection, stress- 


corrosion cracking continued so that the total time 
to failure (without protection) was little, if any more, 
than would have been required had the process not 


been interrupted by the application of cathodic 
protection. 
The evidence indicates that stress-corrosion 


cracking of the AZ31B alloy is primarily an electro- 
chemical process. Segments of a stressed specimen 
are strained at a higher-than-the-average rate 
that the protective film on these crystal surfaces 
ruptured by the first strain is not repaired. These 
film-free surfaces are anodic surfaces on which 
film repair has taken place. An _ electrochemical 
reaction occurs with very rapid attack at the anodes. 
This forms a groove or notch at the anode and acts 
as a region of stress concentration so that the ma- 
terial at the tip of the notch can continue to extend 
at stresses lower than the initial stress. The protec- 
tive film does not re-form on these areas and the 
crack propagates primarily by the electrochemical 
action until failure occurs. 

6. If it is assumed that the alloy has the same 

electrochemical weight as magnesium, it is calculated 
from crack-penetration rates that there is an effective 
current density of 14 amp/cm? at the crack tips. 
7. Itissuggested that the effect of rain on weather- 
exposure specimens was to form a solution with the 
corrosion products that would permit the escape of 
dislocations and hence the rapid extension of micro- 
segments of the specimen as described in 3 above. 
This would result in partial or complete failure in the 
rainwater solution by the process postulated in 5 
above. 


SO 


to 
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Evaluation of Lens Distortion by Visual and 
Photographic Methods’ 


Francis E. Washer, William P. 


Tayman, and Walter R. Darling 


The evaluation of lens distortion by photographic and visual methods is discussed. 
Measurements made on a single lens using the two methods are reported. The precision of 
measurement of each method is determined which shows that the observed differences must 
be attributed to systematic error. Various sources of systematic error are considered. Un- 
compensated differential plate tipping is identified as the most probable cause of the observed 
differences. A method of correction is developed. It is concluded that when work is done 
with extreme care and with due account taken of various insidious sources of error, it is 


possible to achieve comparable results with « 


1. Introduction 


The measurement of radial distortion in the focal 
plane of photographic objectives has been the subject 
of intensive study since the advent of aerial mapping 
from photographs. This particular aberration is of 
prime interest as its magnit ude determines the accu- 
racy with which the final negative maintains the 
correct relationships among the array of point 
images making up the photograph of the corre- 
sponding array of points in the area photographed. 
It was early realized that the reliability of the quanti- 
tative information obtained from a photograph in- 
ereased as the distortion in the taking lens decreased. 
Consequently, the development of improved lenses 
was encouraged with the result that succeeding series 
of new lenses were characterized by ever lower values 
of distortion. 

During this period of change, diverse methods of 
evaluating the distortion of lenses came into use at 
various laboratories. The reason for such diversity 
was primarily the availability of given types of 
measuring instruments in various laboratories. There 
are now three principal methods of measuring dis- 
tortion plus numerous additional methods that are 
either the inverse of one of the principal methods or 
a variation thereof. The nodal slide bench is one of 
the oldest methods; this is a visual method capable 
of high accuracy and is perhaps the most widely 
sed. The photographic method is more recent and 
arose Out of the desire to make measurements under 
conditions approximating that of use. The third 
principal method is the goniometric method which is 
ised to considerable extent in Europe. 

Because of the diversity of methods being used in 
the evaluation of distortion, it seems worthwhile to 
nvestigate the results of measurement made in a 
single laboratory on a number of lenses by a variety 
of methods to determine whether or not the values so 
obtained varied appreciably with method. This is 
being done for several different methods and some 
preliminary results are reported herein for two 
methods. 

The two methods are as follows: 

A. Photographic, precision lens testing camera, 
B. Visual, nodal slide optical bench. 


ieiencintemenieeeiemene 
This work w iS perlormed ip connection with the research project sponsored by 
eU.s Air Force. 


‘ither method. 
2. Methods of Measurement 


Methods A and B have been described at some 
length in the literature. However, for purposes of 
clarity, in the following section, a brief description 
of these methods is given. 


2.1. Precision Lens Testing Camera. Method A 


The precision lens testing camera [1]? shown in 
figure 1 was developed at the Bureau by I. C. Gard- 
ner and F. A. Case. It is one of the earliest successful 
devices developed to measure the performance of 
lenses by photographic means. It consists of a bank 
of collimators spaced at 5-deg intervals having reso- 
lution test charts as reticles. The lens under test is 
mounted at the center of convergence of the collima- 
tor fan and can be aimed at any one of them by ro- 
tation of a carriage which carries the lens holder and 
camera back whereon the photographic recording 
plate is mounted. As presently constituted, the lens 
testing camera has 10 collimators covering a total 
angle of 45 deg. When used to test wide angle lenses, 
the camera is aimed at one of the extreme collimators 
(position 1) and the test made. In order to cover a 
complete diameter, a second test is made with the 
camera aimed at the collimator at the opposite ex- 


Figures in brackets indicate the literature references at the end of this paper. 





Figure 1. Precision lens testing camera (Method A). 


* : 

This photograph shows the camera back with a light metering device in place 
on the plate holder. The values of illumination in the focal plane are read 
from the large meter above the collimator. Use of this device enables the 
operator to adjust for uniform exposure in making the test negatives from 
which the values of distortion, Da, are obtained 
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treme (position Il). It has been found by experience 
that the results obtained from two negatives made 
in this manner are quite reliable. 

To determine the value of the distortion, the best 
row of images is selected on the negative made in 
position |. The separation of all images on this row 


from the image formed by light from the 0° collimator 


is measured. The equivalent focal length (/F 7) is 
then determined from the relations 
f d cot B-, | 
fo=dy Cot Bro, 2) 
: f f 
EFL ——. 2) 
Where d, and dy are the measured distances on the 


negative separating the images formed by the lens 
for the target in the O° collimator from the corres- 
ponding images formed by the lens for the targets in 
the collimators located at B=5° and 10° respectively 
from the 0° collimator. One then determines the 
distances from the 0° image to the other images using 
the relation 
d3=EFL tan 8 

The difference between the distances found with the 
aid of eq (4) and the corresponding measured distance 
for a given value of 8 is the value of the distortion 
at that point. It is positive if the measured location 
of the image ts farther from the central image than 
the computed location. Similar computations are 
made for the negative obtained for position II and 
the results averaged and accepted as final 

It has always been maintained that the probable 
error of values of distortion obtained by this method 
did not exceed +0.020 mm 
arising from error in angle in plate measurement are 
discussed in an earlier publication {2} where it was 
established that errors from this source should not 
exceed +0.013 mm 

In making the negatives, the collimator targets are 
illuminated by light from a tungsten source after 
passing through a K-3 filter. Eastman Kodak 
Spectroscopic plates, emulsion Type V—F, are used to 
the image formed by the lens under test. 
The exposed plates are processed in trays containing 
Eastman Developer D—-19 and 68° F $ min with 
continuous agitation. 


The sources ol error 


record 


2.2. Visual Optical Bench, Direct Nodal Slide. 
Method B 


The visual optical bench has long been the basic 
tool for evaluating the optical constants of lenses. 
The one used at the Bureau has been in existence for 
approximately 30 yr and can still be regarded as a 
precision instrument [3]. For measuring distortion, 
one uses a collimator, nodal slide lens holder shown 
in figure 2, and micrometer microscope. The lens is 
carefully alined in the holder and the axial image 


formed by the lens under test of the illuminated 














nodal slide 


FIGURE 2 Visual 


Vethod B 


reticle of the target is brought into coincidence with 
the object plane of the viewing microscope. Th 
reticle is illuminated by filtered light from a tungsten 
The effective wavelength is approximately { 
575 mu By a series of successive adjustments, a | 
condition is found for which a small rotation of the 
lens about il vertical aXis does not produce a displace- 
ment of the axial image viewed. The rear nodal 
point of the lens is then considered to coincide wit! 
the center of vertical rotation of the nodal slide 
Assuming the equivalent focal length, f, to b 
known, the nodal slide is rotated by amount 8 about 
the vertical axis using the calibrated circle of th 
nodal slide to position it exactly. The entire saddl 
carrying the nodal slide and lens is then moved away 
from the microscope toward the collimator by an 
amount, /(see B | The viewing 
shifted laterally to the new position of the —_— 
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and its dial read and recorded as reading R. The 
nodal slide is then rotated to position 8 and a 
ing of the microscope, L, is made. The 
_ is then obtained from the relation 


dD, (A = L) sec B. (5) 





second sett 
> distortion, Da, 


3. Results of Measurement 


At the time this study was initiated, it was planned 

io make measurements of distortion on a series of 
seven or more lenses by different methods and to 
compare the results. Measurements have been made 
on seven lenses by two different methods and com- 
parative tables of results prepared. All values of 
distortion are referred to the calibrated focal length 
4] for ease of comparison. In general the results are 
comparable and all values fall within the range of 
0.02 mm which has been established here as the 
accepted maximum value of the probable error for 
distortion measurements. However, the agreement 
fell somewhat short of that expected in that system- 
atic differences were observed. Because of these 
systematic differences, it was deemed wise to con- 
centrate on the analysis of the results of measure- 
ment made on a single lens by the two methods, in 
the belief that more information on the causes of lack 
of agreement could be gained thereby. This has 
proved to be the case, and in the following pages 
some of the steps in this investigation and analysis 





are reported. 


3.1. Values of the Distortion by Methods A and B for 
Wide Angle Lens No. 3 


The values of the distortion referred to the eali- 
brated focal length obtained photographically with 
the precision lens testing camera (1,4) and visually 
with the nodal slide optical bench (),) are shown 
in table 1. These values are based on a single run 
by each method. For method A, a single run 
} provides two negatives, one each for positions I and 
| Il; each negative is measured 5 times and_ the 
) results are averaged. For method B, a single run 

consists of two independent sets of measurements; 

the accepted values of distortion are based on the 
It is clear from table 1, that 





average of two sets. 


) TABLE l. Veas ired alues of the distortion versus anqular 
separation B fron the axis for wide anale lens No. 3 
‘ The values obtained with the precision lens testing camera are designated 


D,« while those obtained on the nodal slide optical bench are designated Dg, 
All values are referred to the calibrated focal length and are given in microns 
8 D Dr Dy Da 
le 
0 0 i 
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the values of Dy, and D, do not depart from their 
average by more than 12 « at most and accordingly 
either set of values can be regarded as accurate 
within the usual tolerance of +20 wu. However, the 
differences in the values obtained by the two 
methods D,—D, which are also shown in table 1 do 
not appear to be random but rather appear to 
indicate the presence of systematic error. This is 
shown more strikingly in figure 3 where D,—D, is 
plotted as a function of 8. 
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FIGURE 3 Variation of distortion difference, Dg— Da, with 


angular separation from the axis B. 


3.2. Precision of Measurement for Visual Nodal Slide 


One’s first thought when confronted with dis- 
crepancies in the results of measurement obtained 
by two different methods to make additional 
measurements using the method most likely to be 
suspected. The visual nodal slide optical bench in 
use at the time of these measurements was somewhat 
antiquated so it seemed proper to check its per- 
formance first. Five additional independent sets of 
measurements were made. A pinhole reticle was 
used initially, but a transparent crossline was sub- 
stituted for it during two sets of measurements. 
Some time elapsed between the first two sets of 
measurements and the 3d, 4th, 5th, and 6th cali- 


Is 


bration. During this time, the lens cells were 
removed from the test barrel and replaced. A 


similar operation occurred between calibrations 6 and 
7. Additional care was taken in checking alinement 
in the course of calibration 3 and 7. Each set of 
data was carefully processed and the values of 
distortion adjusted to a calibrated focal length with 
maximum plus at 35° and maximum minus at 45°. 
The results are shown in table 2. Cursory examina- 
tion of these values indicates quite good agreement, so 
an average for all seven including the first two 
calibrations was made. The table also shows the 
probable error of a single determination, and the 
value of PE,. The values of PE, are so low that it 
seems there is no cause to question the precision of 
method B. While there may be a systematic error 
present that impairs the accuracy, none comes to 
mind at the present writing. 
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TABLE 2. Values of the distortion, Dp, versus 8 for seven sets 


of measurement using method B 


These values of De were obtained for wide angle lens No. 3 with the visual 
nodal slide bench All values are referred to the calibrated focal length and are 
given in microns The values of PF, are also shown 
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3.3. Precision of Measurement for Lens Testing 
Camera 


a. First Trial 


Comparison of the average values of distortion 
based on the seven separate determinations shown in 
table 2 with the values of J), for the first two deter- 
minations shown in table 1 indicates that the 
difference between dD, and dD, can not be explained 
by lack of precision in method B. The average 
value of D, for the seven calibrations differ at most 
by 2 u from the average of the first two calibrations. 
Accordingly further tests by method A were made. 
Five new negatives were made in position l on the 
precision lens testing camera and five new negatives 
were made in position Il. Test negatives | to 4 
were made under identical conditions in position | 
were test negatives 5 to S&S in position I] 
The lens was removed, replaced, and realined for 
test before making negatives 9 and 10. Each nega 
tive was measured and the distortion evaluated 
separately and the results are shown in table 3. 


so) also 


Part (a) of table 3 shows the results for position I, 
and part (b) shows the results for position Il. In 
part (c) the results shown in parts (a and (b) are 


combined by pairs and results of this combination 
are shown. The values the index, 
Pk., are somewhat higher in parts (a) and (b 
than those shown in part (c) where the results are 
averaged in the usual manner. The differences in 
the average value of 1), for positions I and II are 
not 


of precision 


regarded as serious as it is presumed to arise 
from a small amount of plate tipping ol a nature 
that negligible on averaging. The final 
values average, D),, shown in part (c) of 


becomes 


of the 


table 3, differ by 1 few microns from the values of 
1), shown in table 1, but have not changed in a 
manner that appreciably reduces the values of 


D,—D,. The values of the index PE 
shown in part (c) compare favorably with the similar 
values shown in table 2 for method B. However, 
the magnitudes of PE. for both methods are 
small to justify the existence of differences in tbe 
values of the distortion as great as are found. 


precision 


too 
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b. Second Trial 


While the first values of the distortion obtained 
in 3.3a compare favorably with the values 
shown in table 1 for method A, it is disquieting lo 
see the large discrepancies that eXist between values 
obtained in positions I and II. While it is probable 
that these discrepancies arise from actual differ- 
ences in performance along the opposing radii, none- 
theless it seemed worthwhile to determine whether 
or not any maladjustment of the lens testing cam- 
era could produce this effect. Careful analysts 
showed that when the camera was properly alined 
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<o as to point properly at the center of the collimator 
reticle for position I that it did not so point when 
swung into position II. Instead there was an error 
of 36 sec in the pointing for position II. Compu- 
tation indicated that this defect in alinement in 
position II could produce small errors in the distor- 
tion although the indicated magnitude was too small 
to change the values of D4 by more than a few 
microns. ‘This alinement error was corrected, the 
angles were recalibrated, 10 additional negatives 
were made, and the measurements described in sec- 
tion 3.3a were repeated. The new results are shown 
in table 4. 


TABLE 4 Values of the distortion, Da, versus B for the second 
trial using method A 
These values of the distortion were obtained for Wide Angle Lens No. 3 with 
recision lens testing camer lowing its recalibration All values are referred 
the ealibr ited focal length and are given in microns In part (a values are 
riven for position I for five negatives toge ther with the average value of D4 and 
the value of PE,. In part (b) similar results are given for position II, and in 
rt yveraged values for paired negatives for positions I and II are given 
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When the final results obtained from negatives 
11 to 20 are compared with those obtained from 
negatives 1 to 10, it is clear that while there is 
slightly better agreement at the large angles than 
before, the agreement is somewhat worsened at the 


small angles. Finally, the average values for 
the second trial, shown in part (c) of table 4, are 


substantially the same as those obtained from the 
first trial, shown in table 3. 


C. Summation 


The final accepted values obtained by methods 
A and B are brought together for comparison in 
table 5. The values of D4 and Dg, differ slightly 
from the values shown in table 1 but are believed 
to be more reliable as each value is the average 
of many more determinations. It is noteworthy 
that the magnitude of AD has increased slightly. 
This table also shows the probable error of the mean 
for the values of Dy, Dg, and AD. Consideration 
of these various values indicates that the systematic 
error still exists and has not been reduced by the 
multiplication of measurements. 


TABLE 5. Comparison of the average values of distortion, Da, 
derived from 20 negatives using method A and the average 
values of distortion, Dg, derived from sevenc alibrations using 
method B. 


Ds 


ilso shown 


Values of the probable error of 
All values are expressed 


is shown as AD= Dz 
of determinations are 


Che difference 
the mean for each set 
in microns 


PF ., for 
8B Da D AD 
Da—Da 
Ds Dr AD 
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30 oo 87 7 1.0 0.4 1.1 
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1) R3 79 4 6 15 1.6 
45 07 106 v 3 0.8 0.9 
3.4. Effect of Plate Curvature 


In view of the known effects of plate curvature 
on the values of distortion, [2] it seemed worthwhile 
to examine the emulsion surfaces of the negatives 
used in method A. This was done and a small 
amount of plate curvature was found. However, 
in no instance was curvature present in sufficient 
amount to produce more than one-fourth of the 
measured differences in distortion values. The 
average departure from flatness for the 20 negatives 
could not produce differences in distortion in excess 
of one-tenth of that found. It may therefore be 
stated that plate curvature is not a prime cause of 
the differences in values of distortion found by the 
two methods. 


3.5. Effect of Plate Tipping 


The plate holder in the precision lens testing 
camera is so constructed that the emulsion surfaces 
are coplanar for positions Land II. For the coplanar 
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or parallel plane condition, small departures of the 
plane from true normality to the optical axis of the 
system under test would not produce variations in 
the average value of distortion even though the 
measured distortion on either side of center would 
be different. However, if a slight warpage of the 
holder has occurred, the plate in position IT would 
not be coplanar with that in position I and a small 
amount of asymmetric distortion would persist and 
adversely affect the final average. Accordingly a 
procedure was developed for checking this possibility 
based on the plate tipping analysis described in an 
earlier paper [5]. The final product of this analysis 
is shown in table 6. When the values of AY) in 
table 5 are divided by the appropriate multiplier ./, 
then the magnitude of fe, can be obtained for each 
value of AD), where e, is the angle between the plane 
of the plate in position I and that in position IT. 
This was done and the results are shown in table 7. 


TABLE 6. Values of tan B and the m iltipl er, WV, for a value 


ol tes 1.00 micron 
rhe multiplier, .M, is so calculated that when it is multiplied t given value 
of fe,, it vields the distortion correction referred to the calibrated ! 
Angular 
eparation tan-g Vl 
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TABLE 7 ( omputation of fe, to determine effective amount of 
encom pe nsated plate { pping that may be present in me thad A 
Values of fe: are obtained for each value of AD giver 
je: for 8 ranging from 10° to 30° is 120 microns; tl value used to deter? 
the probable contribution to distortion Mfe I'he ist colun 
nitude of difference that remains. All values are given 
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The values of fe, that are shown in table 7 indicate 
that some effect of uncompensated plate tipping 
may be present in the values of distortion obtained 
by method A. All values of fe, are of the same sign 
and are of the same order of magnitude except that 
computed for B=40°. The value of AD for this 
point is so small (4u) that little weight can be given 
to the values of fe,, derived from it. While there is 
appreciable variation in the remaining values of 
fe,, it must be remembered that the warping that 
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introduces this effect probably bends the Dlaet 
slightly as well as tipping it. An average fe,, for th. 
range 10° to 30° was determined and found to * 
fe, 129 microns which indicates that the effective 
angle between the plates in positions I and J] _ 
approximately 3.0 min of are. The quantity We 
was then evaluated for each value of 8. The degres 
of compensation achieved by assumption of this 
amount of uncompensated plate tipping is showy 
graphically in figure 4. The solid curve shows thp 
variation of the quantity Mfe, with B while the 
circles show the corresponding values of AD. The 
quantity, D+Mfe,, would equal zero if fe, wag 
invariant for all values of 8. However, the do. 
parture from zero is sufficiently small that it jg 
evident that uncompensated differential plate tipping 
makes an important contribution to the systematic 
error whose cause is being sought. 
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Following this work careful measurements on the 
plate holder were made to determine if this warpage 
actually existed. While measurable values of war- 
page were found that were somewhat smaller than 
might be inferred from the results shown in table 7, 
they appeared large enough to warrant correction. 
Accordingly a new heavier plate holder has been 
constructed and has replaced the one in use during 
this experiment. Recent measurements show that 
the values of AD=D,—D, now obtained are less 
than one-fourth the values of AD obtained prior to 
the change. 

In table 8, the correction Mie, Is applied directly 
to the values of distortion D,. The corrected value 
D,’=Q), Mte, is given in the table and agrees 
closely with the corresponding values of D,. The 
last column shows the differences still remaining 
between values of D,’ and D,. It is clear that the 
remaining discrepancies are so small in comparison 
with the values of D,’ and D, that they can be 


neglected. 
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4. Conclusion 


It is evident from the foregoing study that highly 
accurate values of distortion can obtained by 
either the photographic method using the precision 
lens testing camera or the visual method using the 
nodal slide bench. For each method, values can be 
obtained that are to within +4 microns. 
When proper care is taken, the values obtained by 
either of the methods do not depart from the common 
average by amounts exceeding +5 microns. 

It is however clear that to obtain and maintain 
such high accurarcy it Is necessary be constantly 
on the alert for various insidious sources of 
such as plate curvature, differential plate tipping, 
and incorrect prime calibration of angles used 1 
determinations. 
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error 
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